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Abstract 

This  thesis  examines  the  utility  of  commercial 
satellite-acquired  imagery  for  the  surveillance  of  the 
Canadian  North.  Analytical  performance  models  are  developed 
for  visible  and  thermal  wavelength  sensors.  These  models 
form  the  basis  for  evaluation  of  an  individual  sensor’s 
potential  contribution  to  surveillance.  The  mission  of 
surveillance  is  sectioned  into  five  separate  missions.  For 
each  mission,  sensor  system  evaluation  algorithms,  which 
combine  individual  sensor’s  probabilities  of  detection  and 
tracking,  are  proposed  and  optimization  techniques 
identified.  Sample  algorithms,  using  a  representative 
target  set,  are  provided  for  each  mission.  Analysis  shows 
that  the  selection  of  specific  sensors  is  mission  and 
situation  specific.---  , 
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USE  OF  COnnERCIAL  SATELLITE  IHAGERY 


FOR  SURVEILLANCE  OF  THE  CANADIAN  NORTH  BY  THE 
CANADIAN  ARMED  FORCES 

l_.  Introduction 

Background 

The  first  civilian  satellite  built  to  observe  the  earth 
was  the  National  Aeronautics  and  Space  Administration’s 
CNASA)  Earth  Resources  Technology  Satellite  launched  in 
1972.  That  satellite  mas  later  renamed  the  Land  Satellite 
(Landsat-l)  and  uias  the  Forerunner  of  a  series  of  remote 
sensing  satellites  built  and  flown  by  NASA.  In  1905  the 
Landsat  system  was  transferred  from  NASA  to  the  Earth 
Observation  Satellite  Company  CEOSAT):  a  private  company 
formed  to  market  Landsat  to  the  world.  The  remote  sensing 
of  the  Earth  had  become  a  commercial  enterprise. 

Interest  in  providing  satellite  sensed  data  for  profit 
is  growing.  Systems  Probatoire  d ’ Observation  de  la  Terre, 
or  SPOT,  Joined  E05AT  in  the  remote  sensing  marketplace  in 
1906  with  the  launch  of  its  first  satellite,  SPOT-1.  The 
following  year,  in  1907,  the  Soviet  Union  began  selling 
photographs  taken  from  its  satellites  through  a  commercial 
organization  known  as  Sojuzkarta.  In  the  near  Future  other 
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nations  plan  to  build,  Fly,  and  market  their  own  systems. 
Among  them  are  Japan,  Canada,  India,  and  a  consortium  of 
European  nations  which  participate  in  the  European  Space 
Agency  CESA) . 

Not  only  will  there  he  more  systems  available  in  the 
ISSOs,  but  the  resolution  of  the  data  the  systems  will 
provide  will  be  improved.  For  example,  Landsat  5  carries  a 
sensor  that  achieves  a  ground  resolution  of  30  meters.  A 
similar  sensor  on  Landsat  6  will  achieve  a  ground  resolution 
of  15  meters  C7S:360D. 

As  the  resolution  of  the  imaging  sensors  increases,  so 
does  the  potential  for  use  of  the  data  by  the  military  for 
surveillance  and  reconnaissance.  At  present  SPOT  can 
provide  images  having  a  ground  resolution  of  10  meters. 

This  resolution  applies  irrespective  of  whether  SPOT  is 
looking  at  wheat  fielos  or  military  installations.  While 
military  intelligence  gathering  satellites  provide  better 
resolution  than  does  SPOT,  there  is  sufficient  detail 
available  on  SPOT  images  to  reveal  information  of  military 
value  such  as  the  existence  of  facilities  and  tracks  made  by 
the  gathering  of  vehicles  C21:5).  "SPOT’s  implications  are 
profound  because  it  blurs  the  distinction  between  civilian 
and  military  observation  from  space  in  direct  proportion  to 
the  clarity  of  its  imagery"  (10;3EB). 
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Although  the  trends  toward  more  and  better  data  are 
clear,  the  specific  military  value  of  data  from  commercial 
systems  is  difficult  to  ascertain  —  difficult  since  the 
value  of  available  information  is  largely  determined  by 
comparing  it  to  the  information  required  for  a  specific 
mission  which  can  be  hard  to  identify,  and  difficult  since 
the  information  requirements  for  the  diversity  of  modern 
military  missions  are  equally  diverse. 

Although  there  are  many  responsibilities  entrusted  to 
the  military,  the  primary  mission  of  any  nation’s  armed 
forces  is  the  establishment  and  enforcement  of  control  over 
national  territory,  airspace,  and  territorial  waters.  This 
requires  information  on  the  current  situation  within 
territorial  boundaries,  and  on  potential  threats  to  the 
nation . 

Notwithstanding  the  fact  that  the  enforcement  of 
sovereignty  is  a  common  mission  of  all  militaries  and  that 
this  requires  similar  types  of  information,  the  magnitude  of 
the  task  is  significantly  different  for  different 
militaries.  Perhaps  the  military  most  challenged  by  this 
mission  is  the  Canadian  Armed  Forces  CCAF) .  Canada  is  the 
second  largest  country  in  the  world,  covering  3,920,330 
square  kilometers.  floreover,  approximately  30^  of  its 
territory  is  isolated  and  bears  no  permanent  settlements 
C13:3).  Yet  the  vast  stretches  of  the  Canadian  North  offer 
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economic  and  scientific  wealth,  and  occupy  a  region  of 
significant  strategic  importance  between  the  United  States 
and  the  Soviet  Union. 

The  Canadian  North  is  large  and  it  must  be  controlled. 
To  attempt  to  exercise  control  through  permanent  stationing 
of  forces  throughout  the  area  would  be  impractical.  In 
contrast,  it  is  reasonable  to  maintain  surveillance  of  the 
Canadian  North  so  that  forces  may  be  assigned  as  required. 

The  task  of  maintaining  surveillance  in  the  arctic  is 
itself  a  formidable  one.  However,  it  is  possible  that 
commercial  satellite  imaging  systems  may  be  able  to  provide 
data  of  sufficient  quality  to  support  this  mission. 
Therefore,  it  is  appropriate  to  assess  the  value  of  such 
data  in  terms  of  its  ability  to  provide  useful  surveillance 
information  about  the  Canadian  North  to  the  CAF .  Moreover , 
the  imminent  proliferation  of  commercial  satellite  imaging 
systems  makes  this  a  good  time  to  conduct  this  assessment. 

Problem  Statement 

In  that  the  Canadian  Armed  Forces  CCAF3  requires 
information  concerning  current  situations  and  potential 
threats  in  the  Canadian  North,  and  in  that  surveillance 
information  derived  from  satellite  imaging  sensors  will  be 
available  from  commercial  firms,  the  problem  is  to  determine 
the  value  of  commercial  satellite  images  to  the  CAF  and  to 
develop  a  system  use  strategy  for  the  purchase  of  images. 
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The  following  questions  are  addressed  in  this  research: 
Cl)  Ulhat  commercial  satellite  Imaging  systems  will  he 
potential  suppliers  for  surveillance  and  reconnaissance  data 


during  the  1390s? 

CS)  What  will  the  system’s  capabilities  be  for  these 
potential  suppliers?  Analysis  of  system  performance  will 
include  consideration  of  system  responsiveness  and  coverage, 
and  will  investigate  the  minimum  target  size  detectable  by 
the  system’s  sensors. 

C3)  How  can  commercial  satellite  images  be  evaluated? 

CHt)  Can  weather  satellites  detect  and  track  targets  of 
interest? 

CS)  Will  commercial  systems  require  cuing,  and  will 
they  provide  cuing  for  other  reconnaissance  resources? 

C6)  Will  a  combination  of  several  systems  outperform  a 
single  system? 

C7)  How  good  would  an  alternate  system  have  to  be 
before  it  outperforms  the  best  mix  of  available  commercial 
satellite  imaging  systems? 

CB)  How  much  information  on  Canadian  military  activity 
in  the  arctic  will  be  available  to  commercial  satellite 
surveillance  data  customers? 
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Scope 


The  complete  assessment  of  the  value  of  satellite 
Imagery  for  any  application  Involves  the  broad  areas 
discussed  in  the  next  chapter.  However,  analysis  of  all  of 
these  areas  for  detection  and  tracking  of  targets  in  the 
Canadian  North  is  beyond  the  scope  of  this  research. 

Instead,  this  research  concentrates  on  the  development  of 
algorithms  which  calculate  imaging  systems’  probabilities  of 
detection  and  tracking  targets  of  interest,  on  the  system 
wide  optimization  of  these  probabilities,  and  on  the 
validation  of  the  algorithms  which  are  developed. 

This  is  an  unclassified  thesis.  Representative,  rather 
than  actual,  targets  will  be  used  to  analyze  systems’ 
capabilities . 

Assumptions 

Simplifying  assumptions  may  be  necessary  to  complete 
analytical  calculations  to  determine  systems’  capabilities. 

Also,  although  information  could  be  made  available  to 
the  CAF  from  allied  military  reconnaissance  sources,  such 
sources  will  not  be  considered  as  alternative  collection 
options.  It  is  assumed  that  the  CAF  is  conducting 
unilateral  operations,  and  for  this  reason  foreign  military 
reconnaissance  is  either  not  requested  or  unavailable. 
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Research  Applications 

Although  the  analysis  presented  in  this  thesis  is 
specific  to  the  surveillance  of  Northern  Canada,  the 
algorithms  developed  in  this  research  could  be  applied  to 
determine  the  surveillance  value  of  commercial  satellite 


images  of  any  location. 


1 1 .  Current  Literature 


Organization 

Many  factors  are  involved  in  determining  the  military 
value  of  images  produced  by  commercial  satellite  systems. 

The  problem  is  that  many  of  the  factors  are  dependent  upon 
each  other.  For  example,  the  quality  of  an  image  and  its 
price  are  directly  related;  the  better  the  quality,  the 
higher  the  price. 

To  present  a  coherent  picture  of  these  factors,  it  is 
helpful  to  break  up  this  circle  of  dependencies  and  discuss 
segments  as  if  they  mere  independent.  In  this  way  the 
factors  contributing  to  the  military  value  of  images 
produced  by  commercial  satellite  systems  can  be  discussed 
under  broad  headings.  Four  such  headings  are: 

Cl)  the  quality  of  the  image; 

CS)  the  timeliness  of  the  image; 

C3)  the  cost  of  extracting  useful  information  from  the 
image;  and 

(4)  the  availability  of  the  image. 

IbjB— Qual itu  of .  the  Image 

The  first  major  Factor  to  consider  in  determining  image 
value  is  the  quality  of  the  image.  System  capabilities  are 
commonly  reported  in  terms  of  resolution.  The  problem  is 
that  there  are  several  types  of  resolution  which  can  be 
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Identified  For  satellite  systems,  and  some  of  those  have  a 
variety  of  measures  available.  As  a  result,  "Users  do  not 
understand  the  significance  of  resolution  figures  quoted, 
and  for  many  applications  significantly  overestimate 
capability"  CB9:S).  To  be  of  practical  value,  figures  given 
should  include  the  conditions  under  which  the  resolution  was 
achieved  and  can  be  expected  in  future.  For  example,  a  high 
resolution  system  that  is  pointing  at  the  ocean  may  not 
reveal  very  much,  yet  the  systems  capabilities  are  not  at 
fault.  In  contrast,  the  conditions  can  be  chosen  to  show 
the  system’s  best  capabilities. 

A  more  structured  approach  to  assessing  a  system’s 
capabilities  is  to  break  that  system  down  into  its 
constituent  parts  and  to  handle  the  parts  separately.  This 
sectioning  accomplishes  two  things:  it  reduces  the 
complexity  of  the  analysis  to  a  manageable  level,  and  it 
highlights  the  actual  conditions  under  which  the  system  is 
being  assessed.  In  this  regard  remote  sensing  models 
commonly  consist  of  three  elements:  the  scene  model,  or  that 
which  is  being  looked  at,  the  model  representing  the  medium 
through  which  eiectro-magnet ic  energy  travels  from  the  scene 
to  the  sensor,  and  the  sensor  model  (BB:1B3). 

Scene  Model i no .  The  scene  model  is  made  up  of  the 
description  of  the  features  within  the  field  of  view  of  a 
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sensor.  These  features  can  be  fully  described  by  their 
spatial,  temporal,  and  spectral  characteristics. 

Spatial  characteristics  describe  the  physical 
dimensions  of  features  and  uihere  they  are  in  the  scene. 
Coluell  uirites  that  the  spatial  composition  of  a  scene  can 
range  from  simple  to  complex,  and  has  vegetation,  soil, 
geological,  hydrological,  and  geographical  components 
CIS: 01).  Notwithstanding  this  complexity,  the  scene  should 
be  described  as  completely  as  possible  since  the  actual 
resolution  achieved  by  a  system  is  strongly  scene  dependent 
C50;4) . 

The  temporal  characteristics  of  a  scene  describe  how 
the  scene  changes  over  time  and  are  strongly  research 
dependent.  For  example,  if  remote  sensing  imagery  is  to  be 
used  to  map  stable  geological  formations,  the  scene  can  be 
characterized  as  displaying  little  change  through  time.  In 
contrast,  using  imagery  to  track  drift  ice  involves  a  scene 
that  changes  continuously  through  time.  Colwell  emphasizes 
this  point  since  the  temporal  characteristics  of  a  scene 
will  determine  whether  one  requires  a  single  look  at  an 
area,  or  whether  images  are  required  frequently  ClSiSS). 

The  third  set  of  characteristics  of  a  scene  model 
describes  the  spectral  properties  of  features  in  the  scene. 
Spectral  here  refers  to  the  electro-magnetic  spectrum.  In 
effect  the  features  within  the  field  of  view  of  the  sensor 
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reflect  and  emit  radiant  energy.  They  reflect  energy 
according  to  their  reflectivity,  and  emit  energy  according 
to  their  temperature  and  emissivity . 

In  addition  to  energy  reflected  and  emitted  by  objects 
in  the  scene,  a  third  factor  should  be  included  during 
spectral  analysis,  that  of  the  transmitting  medium  uihich 
lies  betujeen  the  scene  and  the  sensor.  Clearly  this  medium 
is  also  in  the  field  of  vieiu  of  the  sensor.  Therefore,  its 
reflectivity  and  emittance  must  be  taken  into  consideration 
(5:1264) . 

Atmospheric  Model inq.  The  second  element  of  a  remote 
sensing  model  represents  the  medium  between  the  scene  and 
the  sensor.  For  satellite  remote  sensors  that  medium  is  the 
Earth’s  atmosphere. 

The  atmosphere  interacts  with  electro-magnetic 
radiation  in  a  variety  of  ways.  However,  the  two  most 
important  of  these  interactions  are  scattering  and 
absorption.  These  interactions  are  important  since  they  do 
not  affect  all  wavelengths  of  radiation  equally.  For 
example,  "windows"  exist  for  some  wavelengths.  These  are 
specific  regions  of  the  spectrum  for  which  the  atmosphere  is 
essentially  transparent  (11:34).  On  the  other  hand,  there 
are  wavelengths  that  are  strongly  attenuated.  As  a  result 
the  distance  travelled  through  the  atmosphere  by  the 
electro-magnetic  energy  is  important. 


In  addition  to  accounting  For  scattering  and 
absorption,  an  atmospheric  model  should  consider  obscuring 
phenomena  such  as  clouds,  and  whether  the  obscuration  is 
rare  or  persistent  C13:06). 

Sensor  Modeling.  The  third  element  oF  a  remote  sensing 
model  is  the  sensor  model .  Although  imaging  sensors  on 
board  commercial  satellites  can  be  complex  in  design  and 
operation,  they  can  nonetheless  be  placed  into  tuo  general 
design  categories:  radiometers,  and  synthetic  aperture 
radars.  This  First  section  oF  sensor  modelling 
considerations  brieFly  explains  the  operating  principles  and 
how  resolutions  are  determined  For  each  sensor  category . 

The  second  section  shows  trends  in  sensor  systems.  For  a 
review  oF  current  and  planned  sensor  suites,  see  Appendix  A. 

Radiometers .  A  radiometer  Focuses  radiant  energy 
From  everything  within  its  Field  oF  view  onto  a  detector 
array.  This  "image"  is  then  preserved  directly  on  Film  or 
is  converted  into  a  digital  signal  For  subsequent 
transmission  to  a  receiving  site.  In  the  case  oF  digital 
signals,  an  image  is  produced  by  lining  up  the  individual 
picture  elements,  or  pixels,  a  process  which  in  effect 
reconstructs  the  original  scene,  square  by  square  C33:E5). 

If  the  detector  array  is  film,  the  problem  remains  as 
to  how  to  transport  the  information  contained  on  the  Film 
back  to  the  Earth.  To  do  this,  the  film  is  either  scanned 


by  a  second  sensor  and  the  information  converted  into  a 
digital  signal  as  described  above,  or  a  canister  containing 
the  film  is  Jettisoned  from  the  satellite  and  recovered  on 
Earth . 

Radiometer  Resolution.  Radiometer  resolution  is 
described  by  Jensen  as  being  characterized  under  four 
general  headings-,  spectral,  radiometric,  temporal,  and 
spatial  C41 : 4D  . 

The  spectral  resolution  of  a  radiometer  identifies  the 
portions  of  the  electro-magnetic  spectrum  to  uhich  the 
instrument  is  sensitive.  Commonly  imaging  radiometers  are 
sensitive  to  visible  and  near  infrared  radiation.  For 
example,  the  High  Resolution  Visible  CHRV)  radiometers 
operating  on  the  SPOT  satellite  collect  visible  wavelength 
radiation  in  four  channels  or  bands,  .5-. 59  >im,  .61-. 50  m"'i 
.79-. 89  ^m,  and  .51-. 73  CE0:495).  This  defines  the 

spectral  resolution  of  the  HRV . 

The  radiometric  resolution  of  a  sensor  is  its  ability 
to  record  many  levels  of  brightness  in  its  images  Cll:E24). 
Although  this  concept  can  be  applied  to  film  it  is  more 
important  in  systems  that  digitize  their  information  for 
transmission  to  Earth.  In  such  systems,  levels  of 
brightness  are  recorded  as  discrete  values  only.  Depending 
on  the  system  there  are  typically  64,  120,  or  256  values  for 
brightness  available  for  each  pixel.  While  64  levels  of 


brightness  may  appear  to  limit  the  system  only  slightly, 
this  limitation  nonetheless  introduces  errors  in  the  image. 

A  system's  temporal  resolution  is  a  measure  of  houi 
often  it  images  a  given  area  of  interest  (41:5).  Full 
discussion  of  temporal  resolution  is  presented  in  the  review 
of  image  timeliness. 

The  final  measure  of  resolution,  spatial  resolution,  is 
both  the  most  important  for  consideration  and  the  most 
complex  for  discussion.  The  complexity  arises  not  in  the 
concept,  since  spatial  resolution  is  basically  a  measure  of 
the  smallest  linear  separation  between  two  objects  that  can 
be  resolved  (41:4).  The  complexity  arises  since  there  are 
four  different  ways  to  determine  and  report  spatial 
resolution  (69:3). 

While  each  of  the  four  ways  carries  its  own  advantages 
and  disadvantages,  they  all  share  one  limitation.  That 
limitation  is  that  their  ability  to  resolve  objects  is 
calculated  for  objects  that  are  equally  intense  and  distinct 
from  the  background.  In  fact  a  system’s  resolution  is 
maximized  under  such  conditions.  Under  conditions  of 
moderate  contrast,  actual  resolutions  can  be  expected  to  be 
considerably  less  (69:8). 

One  way  to  calculate  spatial  resolution  is  to  consider 
only  the  geometry  of  the  system.  By  this  method,  the 
instantaneous  field  of  view  (IFOV)  is  determined  by  taking 
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the  product  of  the  range  to  the  target  and  the  detector 
size,  and  dividing  this  product  by  the  focal  length  of  the 
optics.  This  equation  is  shown  below. 

Rd 

IFQV  =  -  Cl) 

f 

Uhere  IFQV  =  instantaneous  field  of  view 

R  “  range  to  target  surface 

d  -  individual  detector  size 

f  *  focal  length 

This  calculates  the  nominal  ground  resolution  cell,  or 
the  area  seen  by  a  single  detector  in  the  array.  If  two 
objects  in  the  scene  fall  within  a  single  ground  resolution 
cell,  the  sensor  will  record  a  single  target,  so  the  objects 
will  not  be  resolved. 

The  problem  with  IFQV  calculations  is  that  there  are 
several  factors  which  degrade  this  ideal  capability.  Fusco 
and  Hsu  have  reported  that  actual  Landsat  Thematic  flapper 
CTfl)  data  spatial  resolution  is  between  40  and  50  meters,  in 
contrast  to  the  30  meter  IFQV  C30:1B1).  They  attribute  the 
difference  to  non-sensor  factors  including  geometric  re¬ 
sampling  and  atmospheric  effects.  UJatkins  and  Thormodsgard 
report  similar  effective  instantaneous  Fields  of  view 
(EIFQV)  and  attribute  the  degradation  to  the  sensor  optics. 
However,  they  argue  that  correction  For  the  optics  can 
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restore  the  actual  resolution  to  values  very  near  those 
predicted  by  the  IFOV  C71:S25). 

Again,  contrast  must  be  considered.  The  contrast  ratio 
is  defined  as  the  re  io  of  maximum  irradiance  received  in  an 
image  to  the  minimum  irradiance  received  CB4:206).  Jasani 
reports  that  in  a  scene  in  which  the  contrast  ratio  is  2:1, 
the  actual  spatial  resolution  achieved  is  approximately  2.0 
to  2.4  times  the  IFQV.  It  is  not  until  the  contrast  ratio 
reaches  1000:1  that  this  factor  is  reduced  to  1.4  times  the 
IFOV  C39: 10)  . 

Synthetic  Aperture  Radar.  In  contrast  to  a 
passive  radiometer,  which  relies  on  detection  of  electro¬ 
magnetic  radiation  originating  from  an  external  source, 
synthetic  aperture  radars  CSAR)  are  active  devices  which 
send  out  pulses  of  energy  and  build  an  image  based  on  the 
energy  reflected  from  the  scene.  The  advantage  of  a  SAR  is 
that  it  can  acquire  surface  imagery  anytime  of  the  day  or 
night  irrespective  of  the  surface  solar  irradiance  or  cloud 
cover  C26 : 641 )  . 

The  basic  principles  for  SARs  are  similar  to  those  of 
conventional  radars  and  directly  affect  SAR  resolution. 

SAR  Resolution.  Radars  determine  the  range  to  a 
target  by  timing  how  long  a  pulse  takes  to  return  to  the 
antenna.  It  is  important  to  note  that  the  pulse  of  energy 
that  is  being  sent  out  does  not  go  out  instantaneously,  but 
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requires  time  to  send.  As  a  result,  a  return  received  at  a 
time,  Tl,  may  be  energy  contained  in  the  initial  part  of  the 
pulse  that  uias  reflected  from  a  maximum  range,  energy 
contained  in  the  last  part  of  the  pulse  that  uias  reflected 
by  a  nearby  feature  in  the  scene,  or  some  combination  of 
time  and  distance  in  between  these  extremes.  Thus  the  range 
resolution  is  determined  by  the  time  it  takes  to  send  the 
pulse,  or  pulse  length  of  the  outgoing  signal.  Equation  E 
shows  the  relationship  between  pulse  length  and  range 
resolution  CE9;182D. 


c  T 

delta  R  -  -  C2) 

2  sinC0) 

Where  delta  R  *  range  resolution 

c  =  speed  of  light 

T  »•  pulse  length 

0  =•  depression  angle 

The  bearing  to  a  target,  or  azimuth,  is  determined  by 
recording  the  direction  in  which  the  antenna  is  pointing  at 
the  time  of  the  return  of  the  energy.  The  problem  is  that 
all  of  the  energy  that  is  sent  out  is  not  found  on  a  line 
drawn  normal  to  the  center  of  the  antenna,  hut  is 
distributed  around  that  line.  The  better  the  concentration 
of  the  energy,  the  finer  is  the  beam,  and  the  greater  is  the 
azimuthal  resolution  of  the  radar. 
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Two  things  control  a  radar’s  azimuthal  resolution:  it 
is  directly  proportional  to  the  wavelength  of  the  energy 
sent  out,  and  inversely  proportional  to  the  diameter  of  the 
radar  antenna  C33; 29). 


.7  S  X 

Ra  -  -  C3) 

D 

Where  Ra  *  width  of  the  beam 

S  *  slant  range  to  the  target 
X  -  wavelength  of  radar 
D  *  diameter  of  the  antenna 

To  improve  Ra,  two  options  exist:  to  increase  the  size 
of  the  antenna,  or  to  decrease  the  wavelength  of  the  radar. 
There  are  practical  limitations  to  each  of  these  options 
C33:30).  However,  it  is  possible  to  emulate  the  performance 
of  an  antenna  larger  than  the  actual  antenna.  This  is 
accomplished  by  illuminating  a  target  with  more  than  one 
pulse  of  radar  energy  as  the  aircraft  or  spacecraft  travels 
along  its  track,  and  storing  the  amplitude  and  phase  of  the 
returns  from  each  pulse.  These  returns  can  then  be  combined 
to  provide  azimuthal  resolution  equal  to  that  achievable  by 
an  antenna  as  long  as  the  full  path  flown  by  the  platform 
while  the  target  was  being  illuminated.  In  the  limit  this 
resolution  is  not  dependent  upon  range  to  the  target 
C33:30).  This  is  true  because  of  the  geometry  involved.  As 
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the  range  increases,  the  length  of  the  synthetic  aperture 
also  increases  since  the  target  uiili  be  illuminated  sooner 
and  for  a  longer  period  of  time.  At  its  best  the  finest 
azimuthal  resolution  of  a  SAR  is  half  the  length  of  the 
actual  aperture  of  the  radar  C33:30). 

L 

Ra  -  -  C4) 

2 

Where  Ra  ■  azimuthal  resolution 

L  =  length  of  the  actual  radar  aperture 

This  best  resolution  is  not  always  sought  since  it 
maximizes  the  processing  required  for  imaging,  which  in  turn 
decreases  the  number  of  images  that  can  be  acquired.  In 
practice  only  a  portion  of  the  available  path  is  used. 

Sensor  Trends.  The  first  important  trend  in  sensors  is 
that  there  will  be  more  commercial  sensors  flying  in  the 
future.  By  the  mid-1330s  images  may  be  available  from  eight 
or  more  systems.  Planned  imaging  satellites  are  described 
in  Appendix  A,  and  shown  in  Figure  1. 

A  second  important  trend  in  sensors  is  that  in  spite  of 
speculation  that  higher  resolution  sensors  will  be  employed 
on  platforms  in  the  1330s,  published  plans  are  to  orbit 
sensors  with  resolutions  similar  to  those  that  are  available 
now.  For  example,  Landsat  reports  a  spatial  resolution  of 
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Floure  1.  Commercial  Satellite  Imaging  Systems 

30  meters  for  its  thematic  mapper.  SPOT’s  radiometers  yield  j 

a  10  meter  spatial  resolution  when  operating  in  the 

panchromatic,  or  across  all  color,  mode.  Sojuzkarta,  the 

Soviet  Union’s  marketing  agency  for  satellite  acquired  ( 

images,  boasts  an  average  of  5  meter  spatial  resolution  for 

its  highest  resolution  photographs  CiBiBSB). 

A  third  important  trend  in  sensors  is  that  there  will  j 

be  more  SARs  in  the  future.  The  first  civilian  SAR  was 
SEASAT ,  launched  in  1378.  Since  then  several  SARs  have  been 
flown  on  the  shuttle.  By  the  mid-lSSOs  there  may  be  three  | 
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imaging  SARs,  one  each  on  the  Canadian  radar  satellite 
CRADARSAT),  the  European  Space  Agency’s  Earth  resources 
satellite  CERS-5),  and  the  folloui-on  to  the  Japanese  Earth 
resources  satellite  CJERS-1). 

Summaru  oF  Image  Qualitu 

In  summary,  image  quality  depends  on  the  scene,  the 
atmosphere,  and  the  sensor.  The  variety  of  components 
contributing  to  the  overall  quality  of  the  image  is  shown  in 
Figure  S. 
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Quality 


Floura  2.  Components  Contributing  to  Image  Quality 
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The  Timeliness  oF  the  Image 


The  second  major  Factor  to  be  considered  in  determining 
image  utility  is  the  timeliness  uiith  uihich  the  image  is 
received.  Not  surprisingly  timeliness  can  be  divided  into 
several  parts:  the  time  it  takes  a  system  to  acquire  the 
image,  the  time  it  takes  to  relay  that  image  to  the  Earth, 
and  the  time  it  takes  to  process  and  deliver  the  image. 

Image  Acquisition.  The  time  it  takes  a  system  to 
acquire  an  image  of  a  specific  target  depends  on  hou  quickly 
the  system  can  be  tasked  to  provide  the  image,  and  on  houi 
quickly  the  system  can  respond  to  the  tasking  and  position  a 
sensor  within  range  of  the  target. 

Tasking  a  system  can  take  time,  particularly  if  an 
exchange  of  documents  and  contracts  is  required  prior  to 
acceptance  of  the  tasking.  Once  the  tasking  is  accepted, 
additional  time  may  be  needed  to  re-program  the  sensor. 

System  responsiveness  to  a  tasking  is  determined  by  the 
sensors  on  board  and  the  orbit  of  the  satellite. 

Considering  the  sensors  First,  the  primary  feature  which 
governs  responsiveness  is  whether  the  sensor  is  fixed  or 
able  to  be  pointed.  For  example,  the  HRVs  on  SPOT  can 
collect  data  from  up  to  S7  degrees  either  side  of  the 
satellite  ground  trace.  This  reduces  the  maximum  time  to 
ru-visit  a  target  from  26  to  5  days  020:436). 
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The  satellite’s  orbit  is  also  important  in  assessing 
responsiveness.  Harris  points  out  that  there  is  a  tradeoff 
involved  between  the  spatial  and  temporal  resolutions  of  a 
system  C33:43).  This  occurs  since  higher  resolution  for  a 
fixed  detector  array  requires  a  lower  altitude.  Selection 
of  a  lower  altitude  means  that  the  sensor  will  see  a  smaller 
area,  or  swath,  on  each  pass.  If  a  smaller  swath  is 
observed  on  each  pass,  it  will  take  longer  to  image  the 
entire  surface  of  interest,  and  in  turn  it  will  take  longer 
to  re-visit  a  specific  area. 

Relay  Time.  The  time  it  takes  to  relay  the  image  to 
Earth  depends  on  the  system  strategy  for  that  relay.  Three 
relay  strategies  are  common.-  direct  transmission,  delayed 
transmission,  and  hard  copy  recovery. 

Direct  transmission  is  possible  if  the  sensor  produces 
a  digital  signal  representing  the  image  and  is  in  range  of 
either  a  relay  platform  or  a  ground  receiving  station.  The 
delay  introduced  by  this  strategy  is  measured  in  fractions 
of  a  second  and  is  not  significant  overall. 

Delayed  transmission  is  used  if  the  image  can  be  stored 
electronically  and  if  the  sensor  is  not  in  range  of  either  a 
ground  station  or  a  relay  satellite.  The  delay  introduced 
by  this  strategy  is  measured  in  minutes  since  typical  orbit 
periods  for  imaging  satellites  are  about  100  minutes. 
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To  recover  a  hard  copy  of  the  image  on  Earth  is  the 
most  time  consuming  strategy  of  the  three.  This  involves 
de-orbiting  a  film  canister  and  can  introduce  delays 
measured  in  months. 

Processing  Time.  The  time  required  to  process  data 
depends  on  the  information  needed  from  it.  Processing 
options  and  their  costs  are  presented  in  the  following 
discussion . 

Summaru  of  Image  Timeliness 

The  components  that  contribute  to  the  time  it  takes  to 
receive  an  image  are  summarized  in  Figure  3. 


Figure  3.  Components  Contributing  to  Image  Timeliness 
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The  Cost  of  Extracting  Useful  Information  From  the  Image 

The  third  major  factor  to  consider  in  determining  image 
utility  is  the  cost  of  extracting  useful  information  from 
the  image.  This  section  of  the  literature  review  defines 
useful  information,  the  data  processing  options  that  are 
available  to  assist  in  information  extraction,  and  the  cost 
in  time  and  money  of  those  data  processing  options. 

Useful  Information.  The  usefulness  of  information  can 
be  determined  by  comparing  the  information  available  to  that 
required.  As  noted  earlier,  the  information  required  is 
determined  by  the  specific  mission.  Here  that  mission  is 
surveillance . 

Surveillance  is  defined  in  general  terms  as  watch 
keeping.  For  the  military,  to  effect  surveillance  is  to 
detect,  identify,  and  track  objects  of  interest.  As  a 
consequence  of  this  definition,  for  information  to  have 
surveillance  value,  it  must  contribute  to  the  detection, 
identification,  or  tracking  of  objects  of  interest. 

Qblects  of  Interest.  Objects  of  interest  are  easily 
defined.  In  writing  on  the  requirements  For  real  time 
information,  Allen  points  out  that  if  information  is  needed 
in  real  time,  it  must  be  because  the  situation  is  changing. 
IF  the  situation  is  changing,  the  objects  of  interest  must 
be  moving.  Based  on  this  logic,  he  proffers  a  list  of 
targets  which  includes  tanks,  trucks,  missiles,  guns. 
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troops,  ships,  boats,  and  aircraft  C2:145).  A  similar 
definition  of  objects  of  interest  in  the  Canadian  North  has 
been  made  by  a  Canadian  Government  interdepartmental  task 
force.  Their  report  defines  fishing  boats,  survey  parties 
on  land  or  ice,  base  camps  on  land  or  ice,  drift  stations  on 
ice,  and  all  ships  as  objects  of  interest  C 14 : 112-1E3) . 

Detection  ■  The  first  task  in  surveillance  is  to  find, 
or  detect  an  object  of  interest.  Campbell  defines  detection 
as  determining  the  presence  or  absence  of  features  Cll:07). 
Houjever ,  this  definition  closely  resembles  the  definition  of 
identification.  To  prevent  confusion  in  this  analysis, 
detection  is  used  to  mean  finding  differences  in  the  scene. 
In  contrast,  identification  is  used  to  mean  determining  uhat 
caused  those  differences. 

Finding  differences  in  a  scene  can  be  approached  in  two 
uiays.  The  first  uiay  is  to  analyze  a  single  scene  to  find 
features  that  are  distinct  from  the  background.  This  is 
detection  based  on  the  spatial  characteristics  of  the  scene. 
The  second  uiay  is  to  look  for  changes  in  two  or  more  images 
of  the  same  scene  taken  at  different  times.  This  is  change 
detection  based  on  the  temporal  characteristics  of  the  scene 
C41:234) . 

Detection  may  not  be  a  straightforward  task.  Hord 
writes  that  resolution  and  detection  are  not  the  same  thing. 
For  example,  a  radar  may  detect  an  object  much  smaller  than 


its  ground  resolution  cell  if  the  object  reflects  sufficient 
energy  back  to  the  antenna  C35:3S).  Similarly,  a  small 
bright  source  may  send  out  sufficient  energy  to  raise  the 
overall  response  of  a  passive  detector,  and  thus  be 
detected.  In  each  case,  "sufficient  energy"  is  determined 
by  the  amount  of  signal,  or  energy  from  the  scene,  and  the 
amount  of  the  noise,  uihich  is  the  current  generated  in  the 
detector  by  everything  but  the  scene  being  imaged.  If  the 
signal-to-noise  ratio  is  high  enough,  approximately  3:1  or 
better,  the  target  can  be  considered  detectable  CSS: IE). 

Change  detection  involves  the  comparison  of  tuo  images 
of  a  scene.  This  comparison  can  be  made  before  or  after 
feature  identification  and  classification  have  been  done. 

If  a  comparison  is  made  before  classification,  two 
techniques  are  common:  image  differencing  and  image 
ratioing.  In  image  differencing,  the  first  image  is 
subtracted  from  the  second  and  a  constant  brightness  is 
added  back  in  to  accommodate  negative  values.  In  image 
ratioing  the  brightness  value  for  each  pixel  at  time  1  is 
divided  by  its  brightness  value  at  time  S.  For  both  of 
these  algorithms  the  interpreter  must  decide  what  the 
threshoid  level  is  above  which  a  pixel  is  considered  to  have 
changed  ( 41 : E34-S54 ) . 

Identification .  Although  detection  of  an  anomaly  is  an 
important  part  of  surveillance,  the  heart  of  the  mission  is 
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identifying  the  cause  of  the  anomaly.  In  civilian 
applications  the  task  of  identifying  features  is  referred  to 
as  classification,  and  necessarily  preceeds  the  counting, 
measuring,  and  delineation  of  features  C11:B7).  Very  often 
this  is  accomplished  by  investigating  small  portions  of  the 
scene  in  situ,  selecting  typical  objects  which  can  be  easily 
located  in  the  satellite  image,  and  then  comparing  elements 
in  the  image  to  the  reference  objects  whose  distribution  of 
characteristics  is  known  C60;i28).  The  other  option  is  to 
build  a  "Canopy"  model  based  on  radiative  transfer  theory. 
Such  models  are  usually  validated  by  comparing  predicted  and 
actual  values  of  reflected  and  emitted  radiation  C68:130j. 

The  issue  of  resolution  re-appears  at  this  point.  Some 
authors  argue  that  identification  requires  higher  resolution 
than  does  detection  C24:503.  Others  report  that  too  high  a 
resolution  can  cause  errors  in  classification  if  the  IFDV 
captures  individual  differences  within  a  population 
C19; 115).  When  this  occurs  some  members  of  the  population 
could  be  classified  incorrectly. 

Although  seemingly  otherwise,  these  two  approaches  are 
in  agreement.  This  is  true  since  authors  interested  in 
small  feature  identification  are  united  in  their  insistence 
on  the  requirement  for  high  resolution  systems,  while  low 
resolution  systems  are  supported  by  authors  interested  in 
large  scale  land  use  classification  studies.  Colwell 


confirms  that  selection  of  a  high  or  low  resolution  system 
depends  on  the  specific  study  being  conducted  C13:117). 

Specific  to  small  feature  identification,  it  is 
important  to  consider  the  interplay  of  resolution  and 
identification.  Jasanl  provides  a  very  good  discussion  and 
illustration  of  this  interplay  (39:11-16).  In  short,  the 
better  the  resolution  the  more  specific  the  identification. 

Jensen  states  that  for  identification  the  spatial 
resolution  of  a  system  should  be  less  than  half  the  smallest 
dimension  of  the  feature.  Even  at  that  he  cautions  that 
identification  is  not  guaranteed  if  the  feature  has  the  same 
spectral  response  as  does  its  background  (41:5). 
Representative  spatial  resolutions  required  to  detect  and 
identify  common  military  targets  are  provided  by  Jasani  and 
are  reproduced  in  Table  1.  However,  Doyle  cautions  that  the 
origin  of  the  numbers  is  unclear  and  the  contrast  conditions 
in  the  scene  unreported  CS4:50). 

Tracking .  To  track  an  object  is  to  maintain  a  record 
of  where  that  object  has  been.  Although  it  is  often  given 
as  the  third  task  in  surveillance,  tracking  does  not 
necessarily  follow  identification.  In  fact,  tracking  of 
unidentified  objects  is  a  common  occurrence  in  any 
surveillance  system,  and  often  occupies  the  time  between 
detection  and  identification  of  the  object. 
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Table  1.  Surveillance 

Resolution  Required 

(meters) 

□b lects 

Detection 

Identification 

Description 

Aircraft 

4.6 

1.5 

0.15 

Airports 

6.0 

4.6 

0.3 

Ports/Harbors 

30.5 

15.0 

3.0 

Bridges 

6.0 

4.6 

0.3 

Ships 

7.6 

4.6 

0.3 

Subs  on  Surface 

30.5 

6.0 

0.3 

Troops 

6.0 

2.0 

0.3 

Source: C39: 15) 


Whether  or  not  identification  has  occurred  is  important 
since  the  tracking  of  identified  and  unidentified  objects 
can  be  considerably  different.  Tracking  of  identified 
objects  can  be  through  periodic  contact  with  them,  with  the 
frequency  of  contact  depending  on  the  importance  of  the 
object  and  the  object’s  relative  velocity  in  the  scene.  In 
contrast,  tracking  of  unidentified  objects  requires  contact 
often  enough  to  ensure  that  the  same  object  is  being 
detected  each  time  a  contact  is  made.  Again,  the  object’s 
velocity  relative  to  the  fixed  elements  in  the  scene  mill 
affect  the  frequency  of  detections  required. 

Data  Processing.  Having  defined  useful  information, 
the  information  on  hand  can  be  evaluated  through  comparison. 
This  is  straightforward  except  for  one  thing;  the 
information  on  hand  is  not  static--it  can  be  processed  to 
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facilitate  information  extraction.  For  convenience, 
processing  is  divided  into  three  categories:  pre-processing, 
image  enhancement,  and  image  classification. 

Pre-Processing .  Pre-processing  is  also  referred 
to  as  image  rectification  and  restoration,  and  is  commonly 
undertaken  by  the  system  prior  to  delivery  of  the  image  to 
the  customer.  Pre-processing  seeks  to  accomplish  three 
things:  to  correct  for  geometric  distortions,  to  correct  for 
radiometric  distortions,  and  to  eliminate  noise  C40:611). 

Geometric  distortions  in  the  image  are  the  result  of 
both  systematic  and  non-systematic  effects.  Systematic 
effects  include  panoramic  distortions,  the  distortion 
introduced  by  the  motion  of  the  sensor  as  it  passes  over  the 
scene,  and  the  distortion  introduced  by  the  motion  of  the 
scene  itself  due  to  the  rotation  of  the  Earth.  Systematic 
distortions  are  uiell  understood  and  predictable,  and  can  be 
corrected  by  using  appropriate  algorithms.  In  contrast, 
non-systematic  distortions,  such  as  small  changes  in  the 
platform’s  altitude  or  attitude  C41:10S),  are  not 
predictable.  However,  non-systematic  distortions  can  be 
corrected  by  cross-referencing  the  image  to  ujell-known 
ground  paints  C48:614). 

Radiometric  distortions  arise  as  a  result  of 
differences  in  scene  irradiation  and  errors  in  the  sensor 
operation  C11:B45).  To  correct  for  differences  in  scene 


irradiation,  the  relative  position  of  the  sun  is  determined 
and  the  brightnesses  recorded  are  adjusted  or  normaiized  in 
each  image  acquired.  In  contrast,  sensor  performance  is 
modeiled  radiometricaiiy  both  before  flight  and  during 
flight  so  that  a  standard  correction  can  be  applied  to  all 
images.  Prior  to  launch,  sensor  performance  is  determined 
by  irradiating  the  detector  array  with  a  uniform  light 
source.  The  responses  of  individual  detectors  are  recorded 
and  used  to  correct  subsequent  images.  Similarly,  during 
flight,  absolute  calibration  of  sensors  can  be  determined  by 
presenting  a  standard  target  to  the  sensor.  For  example,  in 
flarch  13SB,  the  HRV  cameras  on  SPOT-1  were  aimed  at  a  test 
ground  site  at  White  Sands,  New  Mexico.  Comparison  of  the 
brightnesses  recorded  by  the  HRVs  to  those  recorded  by  the 
ground  site  and  an  helicopter  borne  radiometer  provided  data 
for  absolute  calibration  CB:66-7B). 

Noise,  as  noted  earlier,  is  false  information  that  is 
inserted  into  the  image  either  by  unwanted  photons  or  by  the 
sensor’s  electronics.  An  important  source  of  unwanted 
photons  is  atmospheric  scattering  and  absorption.  An 
approximation  of  these  effects,  and  therefore  correction  for 
them,  is  passible  since  near  infrared  radiation  is  not 
strongly  scattered  or  absorbed  by  the  atmosphere  C41:97). 

In  the  sensor’s  electronics,  although  noise  is  random, 
individual  detectors  may  report  characteristic  levels  of 
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noise.  A  comparison  of  the  responses  of  detectors  and  their 

I 

neighbors  can  identify  individual  differences  C11:S77). 

Image  Enhancement.  The  goal  of  image  enhancement 
is  to  increase  the  visual  distinction  betiueen  features  in  a 

i 

scene.  Therefore,  image  enhancement  carries  uith  it  the 
implication  that  information  extraction  is  to  be  completed 
manually,  using  the  human  eye  as  the  final  sensor.  There 

I 

are  four  general  techniques  used  to  effect  image 
enhancement;  contrast  manipulation,  spatial  feature 
manipulation,  edge  enhancement,  and  multi-imaging 

( 

manipulation . 

Contrast  manipulation  involves  the  selective  display  of 
only  those  pixels  that  have  specific  values  of  brightness. 

Tuio  examples  of  contrast  manipulation  are  grey  level 
thresholding  and  level  splicing.  In  grey  level  thresholding 
the  analyst  specifies  that  only  pixels  having  a  brightness 
value  higher  than  a  threshold  value  will  be  displayed.  In 
level  splicing  the  analyst  selects  brightness  groupings  and 
specifies  that  all  pixels  having  brightness  levels  uithin 
the  group’s  range  uiill  be  displayed  as  having  the  same 
brightness  C48;B27). 

The  spatial  frequency  of  an  image  is  defined  as  "the 
number  of  changes  in  brightness  value  per  unit  distance  for 
any  particular  part  of  the  image"  C41;13a).  In  spatial 
feature  manipulation  a  filter  is  used  to  emphasize  some 
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features  and  de-emphasize  others  based  on  their  frequency  of 
occurrence.  For  example,  one  low  pass  filter  algorithm 
averages  the  brightness  values  for  a  three-by-three  block  of 
pixels  and  replaces  the  center  pixel  brightness  value  with 
that  average.  In  contrast,  a  high  pass  filter  subtracts  the 
low  pass  filtered  brightness  value  from  the  original 
brightness  value.  The  effect  is  that  low  pass  filters 
preserve  details  which  occur  infrequently  in  the  image  but 
reduce  the  display  of  high  frequency  features,  while  high 
pass  filters  accomplish  the  apposite  (48:637). 

The  third  general  technique  used  in  image  enhancement 
is  edge  enhancement.  Edge  enhancement  operations  sharpen 
the  edges  of  some  features  by  comparing  neighboring  pixels 
and  adjusting  the  brightness  values  in  these  pixels  once  a 
threshold  of  difference  has  been  reached  (11:202).  The 
threshold  is  specified  by  the  analyst.  This  produces 
greater  contrast  for  the  eye  and  can  make  features  more 
prominent  and  easier  to  analyze  (41:144). 

The  last  general  technique  for  image  enhancement 
discussed  here  is  multi-image  manipulation.  Multi-image 
manipulation  does  not  involve  comparison  of  two  images  of 
the  same  scene  taken  at  different  times.  It  refers  to  the 
comparison  of  the  simultaneous  response  of  one  spectral  band 
of  the  radiometer  to  another  (48:650).  This  comparison  is 
made  by  dividing  the  brightness  recorded  in  one  band  by  that 
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recorded  in  the  other,  pixel  by  pixel,  throughout  the  image. 
Thus  the  ratio  of  one  band  response  to  the  other  is 
obtained.  Although  selection  of  the  bands  for  comparison  is 
largely  trial  and  error,  ratioing  has  been  useful  in 
reducing  differences  caused  by  topography,  shadoius,  and 
seasonal  variations  in  surface  irradiation  CHI: 135). 

Image  Classification.  Image  classification,  as 
defined  earlier,  is  the  identification  of  features  in  a 
scene.  Automated  image  classification  is  the  identification 
of  scene  features  using  pattern  recognition  rules  that 
compare  the  features  in  a  scene  to  a  reference  set  of 
spatial  and  spectral  characteristics  C40:BBB). 

Cost .  The  extraction  of  useful  information  from  an 
image  begins  with  receipt  of  the  image  and  may  require  use 
of  the  techniques  described  above.  Clearly,  once  an  image 
has  been  processed  sufficiently  to  allow  for  information 
extraction,  the  processing  need  not  be  continued.  At  that 
point  the  cost  of  the  information,  in  both  money  and  time, 
can  be  determined. 

The  information  cost  in  dollars  is  dependent  upon  both 
the  specific  system  tasked  to  provide  the  image,  and  the 
data  processing  techniques  that  may  have  been  used  to 
extract  the  information.  Taken  together  these  costs  should 
be  evaluated  using  normal  standards  of  cost  effectiveness 
and  in  full  consideration  of  the  alternatives  C15:H7). 
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The  information  cost  in  time  is  again  system  and 
processing  dependent.  Houiever,  it  is  taking  less  and  less 
time  tc  process  images.  In  the  extreme,  five  hours  of 
operation  of  a  SAR  on  a  shuttle  flight  required  one  year  of 
analysis  (“10: 11).  In  contrast,  the  ESA  ERS-1  SAR  plans  to 
be  able  to  provide  fast  delivery  data  uiithin  three  hours  of 
scene  irradiation  C25:537).  Moreover,  Bernstein  predicts 
that  near-real-time  distribution  of  information  will  be 
available  in  the  future  because  of  on-board  processing 
C7:74)  . 

Summaru 

In  summary,  useful  information  is  information  that  is 
needed  to  detect,  identify,  and  track  objects  of  interest. 
Such  information  may  be  embedded  in  an  image  and  require 
further  processing  to  be  extracted.  Overall,  the  processing 
of  data  contributes  to  its  value  but  also  increases  both  the 
time  required  to  extract  the  information  and  the  cost.  The 
elements  contributing  to  the  cost  of  extracting  information 
are  summarized  in  Figure  1. 
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Floura  4.  Components  Contributing  to  the  Cost  of  Extracting 
Useful  Information 
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The  Avallabilltu  qF  the  Image 


"The  passage  of  a  satellite  from  the  space  above  one 
country  to  that  above  another  requires  no  visa,  and 
technology  has  made  obsolete  the  concept  of  complete 
national  privacy"  C54:35D.  However,  there  are  many  things 
that  could  restrict  the  availability  of  images.  Thus  the 
fourth  major  factor  to  consider  in  determining  the  utility 
of  images  is  the  availability  of  those  images. 

First  is  the  issue  of  national  security.  McElroy 
writes  that  one  reason  for  the  early  support  of  weather 
satellite  data  was  that  their  relatively  coarse  resolution 
aroused  no  concerns  about  national  security.  However, 
concern  is  increasing  as  resolution  improves  and  systems 
become  good  enough  to  monitor  military  formations  and 
movements  C54;35,30). 

Second  is  the  issue  of  access.  Given  that  data  has 
been  gathered  and  information  extracted  from  it,  who  should 
have  access  to  this  information?  A  resolution  adopted  by 
the  United  Nations  General  Assembly  in  January,  1987, 
provides  guidelines  for  access  that  include  the  promotion  of 
co-operation  and  mandatary  passage  of  information  to 
affected  states  in  the  case  of  impending  natural  disasters. 
However,  the  resolution  does  not  address  two  situations  of 
concern  to  many  developing  nations:  that  information  of 
military  or  economic  value  could  receive  wide  distribution 


CH0;15),  or  that  system  access  could  be  restricted  due  to 
military  alliances  C54:35). 


Third  is  the  issue  of  cost.  Notuiithstanding  the 
provisions  of  the  United  Nations  resolution  mentioned  above, 
which  calls  for  availability  at  reasonable  cost,  it  is 
possible  that  only  affluent  nations  will  be  able  to  afford 
satellite-acquired  remote  sensing  information  CBIrSB). 

Fourth  is  the  issue  of  denial .  Denial  of  access  to  a 
system’s  data  could  occur  because  of  changes  made 
unilaterally  by  the  owner,  or  through  total  shut  down  of  a 
system  should  it  be  economically  inviable  C54:3B),  In 
either  situation,  an  investment  such  as  a  ground  receiving 
station  could  be  rendered  useless. 

Fifth  is  the  issue  of  risk  to  the  sensor.  Sloup  notes 
that  since  the  military  reconnaissance  satellites  us^d  by 
the  superpowers  are  considerably  more  powerful  than  are 
current  remote  sensing  satellites,  that  superpower  interest 
will  likely  be  in  keeping  track  of  what  others  learn  through 
such  data  CBB:80).  Given  the  inherent  vulnerability  of 
satellite  sensors  CB7:58D,  and  given  that  laser  interference 
with  satellites  is  possible  and  may  have  already  occurred 
C10;8B1),  intentional  dazzling  of  commercial  imaging 
satellites  is  possible. 

Sixth  is  the  issue  of  how  many  of  the  projected  systems 
will  actually  be  launched.  Jasani  and  Creasey  state  that 
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current  data  processing  and  telecommunications  markets  are 
barely  sufficient  to  support  the  cost  of  the  current  range 
of  products  CHO-.IS).  Future  systems  that  will  rely  on 
commercial  income  may  not  find  a  market. 

Seventh  is  the  issue  of  legalities.  Sloup  notes  that 
customary  laiu  has  established  that  space  reconnaissance  is 
legal  as  long  as  it  relates  to  self-defense  C66:79). 

However,  several  authors  suggest  that  while  this  is  true  of 
passive  sensors,  there  is  an  element  of  intrusion  involved 
in  active  sensors  and  that  this  issue  has  yet  to  be  resolved 
CB1:S5,  54:30). 

Eighth  is  the  issue  of  delay,  for  it  may  be  sufficient 
to  delay  provision  of  information  to  invalidate  it.  "For 
perishable  data,  delay  is  as  bad  as  denial"  C54:35). 

aymmam. 

The  issues  contributing  to  the  availability  of 
satellite  images  are  summarized  in  Figure  5. 
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Ill . 


The  review  of  current  literature  shows  the  many 
components  involved  in  the  evaluation  of  remote  sensing 
imagery.  Although  an  in-depth  analysis  of  all  of  these 
areas  is  necessary  to  fully  evaluate  satellite  imagery,  it 
is  not  possible  in  the  time  available  to  complete  this 
research.  For  this  reason,  this  analysis  is  limited  to  the 
evaluation  of  commercial  satellite  systems’  abilities  to 
detect  and  track  objects  of  interest. 

Sensor  Performance  Evaluation 

The  first  question  to  resolve  is  whether  a  sensor  can 
distinguish  an  object  of  interest  from  its  background. 

Thus,  Chapter  IV  develops  a  method  to  determine  the 
probability  that  a  sensor  will  detect  a  target. 

Algorithm  Development 

The  second  question  to  resolve  is  how  does  the  sensor 
compare  to  others  that  are  available.  Since  all  the  systems 
sell  their  information,  the  evaluation  of  their  relative 
worth  is  based  on  the  information  they  can  provide  and  on 
the  cost  of  that  information.  Thus,  Chapter  V  develops 
algorithms  to  compare  systems’  costs  and  information 
performances  and  identifies  techniques  to  select  the  best 
mix  of  sensors  from  those  that  are  available. 
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Analysis 

The  third  question  to  resolve  is  whether  the  algorithms 
that  are  developed  are  valid.  Although  the  quality  of  any 
analysis  depends  on  how  good  the  database  is,  developing 
actual  models  using  representative  targets  is  a  way  to 
investigate  the  face  validity  of  the  algorithms  proposed  in 
Chapter  V.  Therefore,  Chapter  VI  tests  the  optimization 
algorithms  by  building  example  analyses. 
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IV .  Sensor  Performance  Evaluation 


Sagacal 

Whether  a  sensor  distinguishes  a  target  from  its 
background  depends  on  tha  sensor,  the  target,  and  the 
background.  In  simple  terms,  it  is  all  a  matter  of  energy 
since  it  is  the  collection  of  the  energy  which  elicits  the 
response  Csignal)  from  the  sensor.  Current  satellite-borne 
sensors  distinguish  a  target  from  its  background,  or  detect 
a  target,  by  collecting  the  visible  wavelength  energy  that 
the  target  reflects,  or  by  collecting  the  infrared 
wavelength  energy  that  the  target  emits.  Since  reflection 
and  emission  are  different  phenomena,  separate  analysis  is 
required  to  determine  target  detection  probabilities  for 
visible  and  thermal  wavelength  sensors. 

Each  evaluation  of  sensor  performance  begins  with 
analysis  of  the  sensor  under  ideal  conditions.  Ideal 
conditions  means  that  the  target  is  fully  in  view  of  a 
single  detector  in  the  sensor’s  detector  array.  Additional 
complexity  is  added  to  the  analysis  once  the  ideal  case  has 
been  completed. 

Although  visible  and  thermal  wavelength  sensors  are 
different,  they  share  several  performance  degrading 
phenomena.  For  example,  both  are  prevented  from  imaging  a 
target  on  the  ground  if  there  are  clouds  in  the  target  area. 
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Since  this  is  true,  a  single  discussion  of  houj  cloud  cover 
is  accounted  for  in  an  analysis  is  enough.  For  this  reason 
the  discussion  of  common  degrading  phenomena  is  presented 
after  that  which  describes  the  evaluation  of  sensors  under 
ideal  conditions,  and  only  once. 

In  short,  the  topics  covered  in  this  chapter  are.- 
Cl)  Visible  Sensor  Analysis  —  Ideal  Conditions; 

CE)  Thermal  Sensor  Analysis  —  Ideal  Conditions; 

C3)  Non-Ideal  Considerations;  and 
C4)  Obscuring  Phenomena. 

Visible  Sensor  Analysis  -  Ideal  Conditions 

Visible  sensors  are  sensitive  to  radiation  in  the 
visible  wavelength  band,  which  is  approximately  .4  -.7  pm. 
Assuming  that  the  target  itself  is  not  a  source  of  visible 
light,  detection  of  a  target  depends  on  the  collection  of 
solar  energy  reflected  by  the  target.  Since  the  target  and 
its  background  receive  the  same  solar  irradiance,  they  will 
reflect  different  amounts  of  energy  to  the  sensor  if  their 
reflectances  are  different.  If  this  difference  is  large 
enough,  the  target  is  detected. 

An  implicit  definition  has  been  made  here  which  is 
important  for  the  remainder  of  this  analysis.  It  is  that 
"detection"  is  used  to  mean  the  distinction  of  the  target 
From  its  background  in  the  spatial  domain. 


One  measure  of  a  visible  wavelength  sensor’s 
sensitivity  is  its  noise  equivalent  reflectance  difference, 
or  NERD.  This  is  the  difference  between  the  target  and 
background  reflectances  which  is  Just  detectable  given  the 
target  entirely  fills  a  ground  resolution  cell,  the  target 
and  background  cells  are  receiving  equal  irradiance,  and 
there  is  no  atmospheric  attenuation  between  the  scene  and 
the  sensor.  The  NERD  corresponds  to  a  signal-to-noise  ratio 
of  1  under  conditions  of  minimum  scene  irradiance  C6S:405). 

Detection  depends  on  target  size  and  reflectance 
difference  with  respect  to  the  background.  A  target  which 
is  at  least  as  large  as  the  sensor’s  ground  resolution  cell, 
and  having  a  reflectance  difference  of  the  sensor’s  NERD  or 
more  will  be  detected  under  ideal  conditions.  Similarly,  a 
smaller  target  having  a  large  enough  reflectance  difference 
will  also  be  detected.  The  threshold  for  detection  is  that 
combination  of  target  size  and  reflectance  difference  which 
Just  corresponds  to  the  sensor’s  NERD. 

The  relationship  for  threshold  detection  is.  geometric 
and  is  derived  by  comparing  a  ground  cell  that  contains  a 
target  to  one  that  does  not. 

In  Figure  6  the  total  energy  per  second,  or  power, 
reflected  from  cell  A  is  the  product  of  the  solar 
irradiance,  the  background  reflectance,  and  the  cell  area. 
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Figure  6.  Ground  Cells  for  Comparison 


«a 


E_Ru,A_ 

s  o  c 


Where 


power  reflected  From  cell  A 
solar  irradiance 
background  reflectance 
area  of  ground  resolution  cell 


The  total  power  reflected  From  cell  B  is-. 


RtAt^ 


=  power  reflected  From  cell  B 
Aj^  =  area  of  background  not  obscured  by 
target 

Rj.  =  target  reflectance 
A*  =  target  area 


Where 


Detection  occurs  when  the  absolute  value  of 

a  b 

exceeds  the  threshold  as  specified  by  the  NERD,  or  uihen: 


-  Rv 


>  NERD 


C7) 


Where  NERD  •  noise  equivalent  reflectance  difference 


In  the  limiting  case,  threshold  detection  occurs  when 
the  left  side  of  equation  7  is  equal  to  the  NERD.  When  this 
is  true  the  minimum  target  area  which  results  in  detection 
can  be  expressed  as  a  function  of  the  reflectance  between 
the  target  and  the  background. 


A 


t 


NERD 


(B) 


This  general  relationship  still  assumes  a  clear 
atmosphere.  In  fact,  the  atmosphere  attenuates  some  of  the 
power  reflected  by  the  scene.  This  has  the  effect  of 
reducing  the  effective  sensitivity  of  the  sensor.  That  is, 
the  effective  NERD  to  be  used  is; 

NERD 

ENERD  =  -  C3;i 

^a 

Where  ENERD  =  Effective  NERD 

T_  =  atmospheric  transmission 

O 
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Combining  equations  0  and  9,  the  target  size  for 


threshold  detection  is  given  by: 


A 


t 


NERD 


CRt  -  Rfa) 


CIO) 


Using  this  relationship  threshold  detection  curves  can 
be  plotted.  A  typical  plot  is  shown  in  Figure  7.  Targets 
falling  above  the  threshold  detection  line  will  be  detected; 
those  falling  below  will  not. 
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Figure  7.  Generic  Threshold  Detection  Curve 
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This  analysis  makes  assumptions.  They  are: 

Cl)  Backgrounds  and  targets  are  Lambertian; 

C2)  Background  and  target  reflectances  are 
uniform  throughout  the  bandpass  of  the  sensor; 

C3)  Atmospheric  transmission  is  uniform  across 
the  bandpass  of  the  sensor; 

(4)  Sensor  response  is  linear  throughout  the 
bandpass  of  the  sensor; 

C5)  The  target  is  not  moving  uith  respect  to 
the  background;  and 

C6)  The  target  falls  in  a  single  ground 
resolution  cell. 


Thermal  Sensor  Analusis  -  Ideal  Conditions 

Thermal  detection  of  targets  depends  not  on  the  energy 
reflected  by  a  target  but  by  the  energy  emitted  by  the 
target.  Every  object  warmer  than  O’K  radiates  energy.  The 
amount  of  energy  radiated  is  not  uniform  across  all 
wavelengths  and  is  a  function  of  its  temperature.  For  a 
blackbody  the  energy  each  second,  or  power,  per  square  meter 
per  unit  wavelength  interval  radiated  is  given  by  the  Planck 
function  (53:7)  . 

3.74  X  10®  UJ 

n^CA.T)  =  -p - - -  — -  Cii; 

X^CexpCl.44  X  lO^/M)  -  1)  m^  ;im 

Where  \  =  wavelength  in  >im 
T  ”  temperature  C'K) 
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The  poujer  emitted  by  an  object  is  distributed  across 
all  wavelengths.  Figure  B  shows  the  spectral  blackbody 
power  per  square  meter,  or  exitance,  as  a  function  of 
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Figure  0.  Blackbody  Exitance  Curves  For  Objects 
at  300’K  Ctop)  and  250’K  Cbottom) 


wavelength  for  two  objects,  one  at  300’K,  and  one  at  E50*K. 

Since  exitance  is  wavelength  dependent,  and  since  a 
sensor  responds  only  to  a  portion  of  the  electro-magnetic 
spectrum,  the  power  output  in  a  specific  spectral  bandpass 
is  required  For  analysis.  This  is  calculated  by 


SE 


integrating  the  Planck  function  over  the  wavelength  interval 
for  which  the  sensor  is  sensitive. 

ncT)  -  f  x,T)  dx  ciS) 

XI 

The  total  power  emitted  by  an  object  in  a  specific 
bandpass  is  then; 


Bq  =  nCT)€Aj^  C13) 

Where  €  =  object  emissivity 

=  surface  area  of  the  object 

Given  that  the  sensor  response  is  linear,  a  ground  ceil 
containing  a  target  will  be  detected  Cthat  is,  register  a 
different  output  current)  if  the  sensor  receives  a 
sufficiently  different  amount  of  power  from  that  cell  as 
compared  to  a  cell  in  which  no  target  is  located.  One 
measure  of  how  much  power  is  sufficient  for  detection  refers 
to  the  source  of  the  power  difference,  temperature,  and  is 
the  noise  equivalent  temperature  difference,  or  NETD.  This 
is  the  minimum  temperature  difference  between  cells  which  is 
Just  detectable  by  a  sensor,  given  that  the  cells  are  at  a 
uniform  temperature  and  the  background  or  reference 
temperature  is  a  specific  value.  For  example,  the  Landsat  5 
Tn  thermal  band  is  reported  to  have  a  NETD  of  .5*K  at  300*K 
background  temperature  CB'ii'iSB). 
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One  problem  in  using  the  NETD  as  a  measure  of  sensor 
sensitivity  is  that  the  conversion  to  sensor  sensitivity  for 
other  background  temperatures  is  not  direct.  However,  the 
approximate  NETD  for  other  background  temperatures  can  be 
derived  based  on  the  elements  which  contribute  to  the 
signal-to-noise  ratio. 

In  general  the  signal  to  noise  ratio  for  a  sensor 
operating  in  the  linear  response  portion  of  its  response 
curve  is  determined  by  considering  the  arrival  rate  of 
photons  at  the  focal  plane,  and  on  the  sensor’s  electronics. 

Signal  output  from  the  sensor  is  dependent  upon  the 
arrival  of  photons.  The  relationship  is  CBS: EES) : 

ig  *  C#  photons/unit  time)K2^  CIH) 

Where  ig  “  signal  current 

#  photons  *  number  of  photons  incident 

*  constant  (charge  per  electron  times 
the  quantum  efficiency  of  the 
detector) 

Detector  noise  results  from  inherent  detector 
characteristics  such  as  sampling  rate  noise,  dark  current, 
and  shot  noise.  Dark  current  is  the  current  which  flows  in 
the  detector  when  its  field  of  view  is  dark.  Shot  noise  is 
associated  with  the  random  arrival  of  photons.  In 
simplified  form  detector  noise  is  given  by  equation  15. 


5^ 


*  CCshot  noise)^  +  other  sources ) CIS) 

Where  ”  noise  current 

Shot  noise  is  calculated  by: 

‘shot  '  CIS) 

Where  i^  =  average  signal 

e  *  charge  on  one  electron 
T  =  detector  integration  time 

Combining  these  equations: 

S  C#  photons* )K. 

-  *  - ;; - - - —  C17) 

N  CCC#  photons  JKj^e/r)'^  +  (other  sources)*^]  ‘ 

*  per  unit  time 

As  a  result  of  this  relationship  5/N  behaves 
differently  as  a  function  of  the  arrival  rate  of  photons. 

If  feui  photons  are  received,  then  the  total  noise  is 
dominated  by  other  sources  and  shot  noise  can  be  discarded 
from  analysis.  Conversely,  if  many  photons  are  received 
shot  noise  dominates  and  other  sources  of  noise  can  be 
discarded  from  analysis. 

This  gives  rise  to  three  distinct  portions  of  the  S/N 
verses  number  of  incident  photons  curve.  Beloiu  a  specific 
arrival  rate,  the  curve  is  linear.  This  is  the  portion  in 
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uhich  shot  noise  can  be  discarded  and  total  noise  is 
effectively  constant.  Above  a  specific  arrival  rate,  shot 
noise  dominates  and  S/N  increases  as  the  square  root  of  the 
arrival  rate.  Betuieen  these  tails,  shot  and  inherent  noise 
are  important.  This  relationship  is  shown  in  Figure  9.  In 
the  figure,  area  1  is  the  region  which  behaves  linearly. 
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Figure  9.  Representative  Signal  to  Noise  Curve 


As  noted  earlier,  NETDs  are  reported  For  a  specific 
temperature,  one  which  is  in  the  linear  response  portion  of 
the  curve.  For  the  thermal  sensor,  this  temperature 
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corresponds  to  a  specific  number  of  incident  photons.  As 
shown  in  Figure  9,  arrival  of  fewer  photons  than  the  number 
equivalent  to  the  NETD  should  result  in  a  linear  response  by 
the  sensor.  In  this  region  of  sensor  response  Carea  1) 
inherent  noise  dominates  and  can  be  assumed  to  be  reasonably 
constant.  Since  the  noise  is  approximately  constant,  and 
since  the  NETD  identifies  threshold  S/N  for  detection,  the 
signal  required  for  threshold  detection  remains 
approximately  constant  in  this  region. 

From  the  previous  discussion,  NETD  can  be  considered  in 
terms  of  numbers  of  photons.  It  can  also  represent  a 
specific  power  difference.  That  is,  an  NETD  of  S*K  0  300*K 
is  equivalent  to: 

NEPD  -  ®C302*K)  -  BC300*K)  CIB) 

Where  NEPD  =  noise  equivalent  power  difference 

fflC30S*Ki  =  power  received  in  sensor  bandpass 
from  object  filling  the  field  of 
view  at  30P*K 

fflC300*KD  =  power  received  in  sensor  bandpass 
from  object  filling  the  Field  of 
view  at  300*K 

In  the  linear  response  portion  of  the  curve,  when 
referencing  a  background  ceil  that  is  at  a  lower  temperature 
than  the  NETD,  if  the  power  difference  between  two  ceils 
differs  by  more  than  the  NEPD,  the  sensor  will  register  a 
different  response.  IF  that  difference  is  caused  by  the 
target,  then  the  sensor  has  detected  the  target. 


Using  the  NEPD  derived  above  it  is  possible  to  generate 


threshold  detection  curves  For  thermal  sensors.  As  in  the 
derivation  of  visible  wavelength  sensor  detection  curves, 
the  analysis  begins  by  considering  two  ground  cells,  one 
that  contains  a  target  and  one  that  does  not.  This  is  shown 
in  Figure  10. 


Figure  10.  Ground  Cells  For  Comparison 


The  power  emitted  in  the  sensor  bandpass  by  cell  A  is: 


•a  ■  ‘13) 

Where  ^  spectral  exitance  for  the  cell  at 

temperature  T^  in  the  sensor  bandpass 
A^  =  area  oF  the  cell 

=  emissivity  of  the  background 
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The  poiuer  emitted  in  the  sensor  bandpass  by  cell  B  is: 


•b  ■  "b''^b>''r«b  '20) 

Where  "  spectral  exitance  for  the  target  at 

temperature  in  the  sensor  bandpass 

*  area  of  the  target 

■  emissivity  of  the  background 

*  area  of  background  not  obscured  by  the 
target 


Detection  occurs  when  the  absolute  value  of  -  Bg  is 
greater  than  or  equal  to  the  NEPD,  or; 

I  "b'l'b>''c«b  -  -  "b'^b^^r^b  I  -  '"2PD  CEl) 

The  target  size  for  threshold  detection  is: 


A 


t 


NEPD 

"b'^b’^b  -  I 


CEE) 


Again  the  atmospheric  attenuation  must  be  considered. 


ENEPD 


NEPD 


CS3) 


Where  ENEPD  =  effective  NEPD 

T  =  atmospheric  transmission 

a 


Therefore  threshold  detection  is  calculated  using; 


NEPD/t. 


CE4) 
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Using  this  relationship,  threshold  detection  curves  as 
a  function  of  target  exitance  can  be  plotted.  Once  again  a 
combination  of  target  size  and  exitance  which  falls  below 
the  threshold  line  will  mean  that  the  sensor  will  not  see 
the  target  against  the  background.  However,  unlike  the 
curves  plotted  for  the  visible  wavelength  sensors,  thermal 
threshold  detection  curves  are  both  sensor  and  situation 
specific  and  cannot  be  generalized  for  all  targets  in  all 
situations.  For  example,  the  emissivity  of  the  target  will 
shift  the  curve  up  or  down  the  exitance  axis.  This  is  shown 
in  Figure  11  in  which  the  detection  curves  for  three  target 
emissivities  are  plotted.  Because  of  this  variability, 
thermal  detection  curves  apply  only  for  the  situation  for 
which  they  were  generated. 

Non-Ideal  Considerations 

The  detection  curves  developed  to  this  point  are  valid 
for  targets  which  are  completely  contained  in  a  single 
ground  resolution  cell.  However,  a  target  may  fall  in  the 
Junction  of  several  cells.  When  this  happens,  a  target 
which  is  larger  than  the  threshold  detection  size  may  not  be 
detected  since  it  could  present  a  small  enough  profile  to 
each  individual  detector  to  escape  detection. 

In  the  extreme  case,  a  target  could  be  equally  shared 
by  Four  ground  resolution  cells.  In  this  situation,  the 
minimum  actual  target  size  needed  to  guarantee  threshold 
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Figure  11.  Threshold  Detection  Curves  for  Several 

Target  Emissivities 


detection  is  four  times  the  size  of  the  target  under  ideal 
conditions . 

A  second  complication  is  encountered  when  there  are 
gaps  between  detectors  in  the  focal  plane  of  the  sensor  as 
shown  in  Figure  IE.  Again  the  effect  is  to  increase  the 
size  of  the  target  which  would  always  trigger  detection. 

The  actual  target  size  for  threshold  detection  in  all  cases 
is  calculated  using  equation  E5 . 
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Figure  IS.  Target  Overlapping  Gaps  in  Detector  Array 

^th  “  "^^ti  2XCA^.)-^  +  SYCA^^)-^  +  XY  CSS) 

Where  A^j^  =■  target  size  for  threshold  detection 

A^j^  =  target  size  for  threshold  detection  under 
ideal  conditions 
X,Y  =  gaps  between  detectors 

Finally  it  is  possible  that  a  sensor  could  detect  a 
target  under  ideal  conditions,  but  will  not  if  the  target 
overlaps  several  pixels.  Detection  of  the  target  now 
becomes  stochastic.  The  probability  that  the  target  will  be 
detected  can  be  approximated  by  considering  the  size  of  the 
target  and  the  size  of  the  threshold  detection  target. 

The  probability  of  detecting  the  target  is  one  if 
enough  of  the  target  falls  in  the  active  sensing  portion  of 
a  single  detector  in  the  detector  array.  The  probability 


that  enough  of  the  target  will  fall  in  the  active  portion  of 
the  detector  is  calculated  geometrically. 


Figure  13.  General  Detection  Probability  Geometry 


Consider  Figure  13.  For  the  target  to  be  detected  the 
center  of  the  target  must  fall  within  the  active  area  of  the 
detector.  If  the  target  Falls  on  the  edges  of  the  detector 
or  the  gaps  between  detectors  such  that  the  target  area 
remaining  on  the  active  detector  is  less  than  the  threshold 
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target  size  For  detection,  the  target  is  not  detected.  In 
that  the  center  of  the  target  is  equally  likely  to  fall 
anywhere  in  the  detector  array,  the  probability  that  the 
target  will  be  detected  is  the  area  that  would  result  in 
detection  divided  by  the  total  area. 


P 


det 


det 


tot 


CSS) 


Where  P^et  ""  Pr°tjability  of  detection 

^det  *  3rea  which  results  in  detection 
^tot  *  total 


The  total  area  in  which  the  center  of  a  target  could 
Fall  is  the  size  of  a  single  detector  plus  half  the  area  of 
the  gaps  which  surround  individual  detectors.  Considering 
square  detectors  with  equal  gaps  on  all  sides,  the  total 
area  is: 

Atot  “  ce?) 


Where  W  =  length  of  one  side  of  detector 
d  =  gap  between  detectors 


Again  the  geometry  is  used  to  solve  the  area  in  which 
the  target  will  be  detected.  In  Figure  14: 

Cl)  X  is  the  center  point  of  the  gaps  in  the  detector 
array ; 

C2)  Z  is  the  center  of  the  target  which  could  lie 
anywhere  along  the  diagonal  line;  and 

C3)  a  is  the  dimension  of  one  side  of  the  area  of  the 
target  which  lies  on  the  active  portion  of  the 
detector . 
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Figure  14.  Target  Gap  Overlap  and  Detection  Probability 


From  this  Cb+cD  is  equal  to  half  the  gap  between  detectors. 
For  detection,  the  area  remaining  in  view  of  the  detector 
Ca  )  must  equal  or  exceed  the  threshold  detection  size  under 
ideal  conditions  CA^^).  Also,  the  target  size  is  equal  to 

P 

four  times  the  square  of  size  Ca+b)  .  Solving  for  b: 


5 

b  *  -  -  CA^..)-=*  C2B) 

S 

The  area  in  which  the  center  of  the  target  may  fall  and 
still  result  in  detection  is  then: 

A.„^.  =  CU)+2b)^  C2S) 


Obscuring  Phenomena 


Having  considered  the  internal  complications  which 
occur  because  of  the  physical  layout  of  the  detector  array, 
it  is  time  to  consider  complicating  factors  which  are 
outside  the  sensor.  Such  factors  have  the  effect  of 
obscuring  the  sensor.  In  general  two  major  types  of 
obscuring  phenomena  are  found  in  the  Canadian  arctic: 
darkness  and  cloud  cover. 

Darkness .  Darkness  is  predictable  and  is  a  function  of 
latitude  and  time  of  year.  Of  interest  are  the  local 
sidereal  times  that  the  sun  rises  and  sets  at  a  given 
latitude.  Since  the  orbits  for  imaging  satellites  are  sun- 
synchronous,  they  overfly  the  same  latitude  at  the  same 
local  sidereal  time  each  pass  Coniy  the  longitude  changes) . 
As  a  result  the  time  when  the  satellite  is  overhead  a  given 
latitude  is  predictable.  If  the  satellite  is  overhead 
during  darkness,  and  if  the  sensors  on  board  require 
sunlight  for  target  irradiation,  the  obscuration  by  darkness 
is  complete  and  the  probability  of  detection  is  reduced  to 
zero.  On  the  other  hand,  if  the  satellite  flies  over  the 
target  while  the  target  is  being  irradiated  by  the  sun  the 
obscuration  by  darkness  is  zero. 

Although  there  are  times  when  the  sun  has  risen  but  is 
not  yet  high  enough  in  the  sky  to  provide  enough  irradiation 
of  the  ground  for  imaging,  this  analysis  assumes  there  is  no 
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middle  ground.  IF  the  sun  Is  up  the  target  can  be  Imaged. 


Therefore  the  Irradiation  factor,  If^^i  can  only  take  on  the 
values  zero  or  one . 

The  time  of  sunrise  and  sunset  Is  determined  by 
considering  the  right  ascension  and  declination  of  the  sun, 
and  the  latitude  of  the  observer  C'lrAlS).  The  formula  is; 


LST  -  .997S7Co(  ±  cos  ^C-tan0  tan?T)l  C301 

Where  LST  -  local  sidereal  time 

■:«:  =  right  ascension 

0  =  latitude  of  observation 

TT  =  declination 


Table  2.  Local  Sidereal  Times  for  Sunrise  and  Sunset 
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The  times  of  local  sidereal  sunrises  and  sunsets  are 
published  annually  in  the  Astronomical  Almanac.  Although 
there  Is  some  variation,  the  approximate  time  of  sunrise  for 
a  given  date  at  a  given  latitude  is  the  same  year  to  year. 
For  this  reason,  it  is  reasonable  to  use  the  almanac’s  times 
of  sunrise  and  sunset  throughout  the  analysis  period. 

However ,  the  almanac  only  provides  times  for  latitudes  up  to 
and  including  66*N.  To  determine  local  sidereal  times  for 
sunrise  and  sunset  at  latitudes  not  listed  in  the  almanac, 
the  formula  shown  in  equation  30  was  used.  The  times  for 
sunrise  and  sunset  are  given  in  Table  E. 

Determination  of  the  local  sidereal  time  of  satellite 
passage  over  a  given  latitude  requires  orbital  analysis. 
Since  an  imaging  satellite  Flies  in  a  circular  orbit,  its 
ground  speed  is  constant  throughout  its  orbit.  If  the  orbit 
is  inclined  at  SO*,  the  time  taken  to  travel  between  any  two 
latitudes  is  calculated  using  equation  31 . 

Tlat  “  CCBi-|3g)/360)Tp3^  C31) 

Where  ^lat  ~  time  to  travel  between  the  latitudes 

^2^  0  =  First,  second  latitude 

T„l^  =  orbital  period  of  the  satellite 
per 

A  satellite  in  an  orbit  at  any  other  inclination  has 
more  ground  to  cover  before  it  reaches  the  latitude  of 
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interest.  This  distance  is  calculated  by  considering  both 
the  latitude  of  interest  and  the  orbit  inclination.  Since 
the  satellite’s  ground  speed  is  constant,  the  time  to  travel 
betueen  the  equator  and  the  latitude  of  interest  becomes: 


flat  ■  ^ 


sinCLat ) 


sinC i ) 


per 

360* 


C3E) 


UJhere  Lat  ”  latitude  of  interest 


Having  ujorked  out  the  time  for  passage  between  the 
equator  and  a  given  latitude,  the  local  sidereal  time  the 
satellite  passes  over  this  latitude  is  its  equatorial 
crossing  plus  or  minus  the  transit  time.  The  plus  or  minus 
is  assigned  depending  on  whether  the  satellite  is  ascending 
or  descending  at  time  of  equatorial  crossing. 

Armed  with  the  local  sidereal  times  for  sunrise, 
sunset,  and  satellite  passage  overhead,  the  rule  for 
assigning  the  irradiation  factor  is: 


fac 


If  t 


sr 


<  t 


pass 


<  t 


ss 


=  0  Otherwise 


UJhere  t^^  =  time  of  sunrise 

t  =  time  of  satellite  passage 

tgg  =  time  of  sunset 


C33) 
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Cloud  Cover.  In  contrast  to  the  predicicability  of 
darkness,  analysis  of  the  second  obscuring  phenomenon, 


clouds,  is  stochastic.  For  those  sensors  which  will  suffer 
obscuration,  the  obscuration  is  taken  to  be  complete  when 
clouds  are  present.  Therefore,  the  information  needed  is 
the  probability  that  there  will  be  cloud  cover  in  the  target 
area.  Once  this  is  found,  the  probability  of  getting  an 
image  of  the  target  area  considering  the  cloud  cover  is 
given  in  equation  34. 


^obscl  ^  ^clcud^ 


C34) 


Where 


P 

obscl 

p 

cloud 


probability  of  observation  considering 
clouds 

probability  of  cloud  cover 


Conclusion 

In  summary  this  analysis  shows  that  the  overall 
probability  of  detecting  a  target,  given  that  it  is  in  the 
field  of  view,  is  the  sensor’s  detection  capability  under 
ideal  conditions  modified  to  include  the  complications  of 
array  geometry  and  external  obscurations. 

The  obscuring  phenomena  of  interest  for  visible 
wavelength  sensors  are  clouds  and  darkness.  Thus  for 
visible  wavelength  sensors,  the  expression  for  P(jet*’ 


0 


is  the  probability  of  detection  in  a  single  image 
considering  obscuring  phenomena,  is: 


^det*  "  ^det^fac*'^  ^cloud^ 

The  obscuring  phenomenon  of  interest  for  thermal 
sensors  is  cloud  cover .  Thus  for  thermal  sensors  the 
expression  for  given  by: 

^det*  ”  ^det ^ ^~^cloud ^ 
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General 

The  second  question  to  resolve  is  hoiu  a  sensor  compares 
to  the  other  sensors  which  are  available.  Although 
individual  comparisons  are  interesting,  it  is  more  important 
to  be  able  to  compare  all  the  sensors  in  terms  of  their 
contribution  to  the  mission  at  hand.  Conceptually  there  is 
a  system  of  commercial  sensors.  The  goal  of  this  chapter  is 
to  develop  the  rules  for  optimization  of  the  system’s 
performance  through  selection  of  the  sensor  or  mix  of 
sensors  which  provides  the  greatest  mission  accomplishment. 

Sensor  comparisons  are  based  on  the  individual  sensor 
detection  capabilities  as  discussed  in  Chapter  IV.  However, 
the  previous  chapter  only  shows  how  a  single  sensor’s 
probability  of  target  detection  can  be  calculated  given  that 
the  target  is  in  the  field  of  view  of  the  sensor.  How 
likely  is  the  target  to  be  in  the  field  of  view?  The  answer 
is  that  it  depends  on  the  mission  and  on  the  sensor.  For 
this  reason  this  chapter  builds  mission  specific  algorithms 
that  are  used  to  assess  the  system  nerf ormancu . 

Detect  j^on 

The  first  mission  in  surveillance  is  detection. 

However,  detection  itself  can  be  broken  into  two  distinct 


missions.  The  first  mission  is  cued  detection.  Can  the 
system  of  available  sensors  detect  a  target  given  the  system 
knows  where  to  look?  Information  from  a  source  suggests 
that  a  target  may  be  at  a  location.  The  measure  of 
performance  for  cued  detection  is  the  probability  that  the 
system  will  image  the  target  and  confirm  or  deny  its 
presence.  The  second  mission  is  non-cued  detection.  Can 
the  system  detect  a  target  if  it  does  not  know  where  to 
look?  In  non-cued  detection,  the  system  is  tested  as  a  wide 
area  surveillance  tool.  The  measure  of  performance  for  non- 
cued  detection  is  the  probability  that  the  system  will 
detect  a  target  given  that  the  target  is  in  the  area  of 
interest  but  the  location  is  unknown. 

Phenomenologically  cued  detection  is  a  subset  of  non- 
cued  detection.  Thus,  non-cued  detection  is  the  more 

demanding  of  the  two  detection  missions.  For  this  reason  ^ 

the  algorithm  to  calculate  the  probability  that  the  system 
will  detect  a  target  under  non-cued  conditions  is  developed 
first.  I 

Non-Cued  Detection  Algorithm  —  Stationari-i  Targets 

In  general,  non-cued  detection  is  pursued  over  an 
extended  period  of  time  during  which  many  images  of  the  ' 

total  area  of  interest  could  be  acquired.  Given  that  the 
images  that  are  acquired  do  not  overlap  and  that  the  target 

I 

is  not  moving,  the  probability  of  detecting  the  presence  of 
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a  target  in  the  area  of  interest  depends  on  the  probability 
that  the  target  will  he  detected  in  a  single  image,  the  size 
of  one  image  as  compared  to  the  total  area  in  which  the 
target  could  be,  and  on  the  number  of  images  acquired.  The 
relationship  for  a  given  sensor  is  given  in  equation  37. 


dot 


^one^im^det* 

^int 


C37) 


Where 


^dot 


'^int 


probability  of  detection  over  time 
area  of  a  single  image 
number  of  images  acquired 
probability  of  detection  in  single 
image 

total  area  of  interest 


The  number  of  images  that  are  acquired  will  depend  on 
the  time  the  target  remains  in  the  area  of  interest. 


^im 


^ut^tav 
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Where  =  #  images  taken  per  unit  time 

T^g^  =  time  target  is  in  the  area 


Also  the  maximum  number  of  images  taken  in  a  given  unit 
of  time  will  be  constrained  by  the  cost  of  the  images,  or: 


( 


I 


N 


ut 


C39) 


Where  =  budget  per  unit  time 

=  cost  of  a  single  image 
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Folding  equations  37,  38  and  39  together  yields  a 


general  expression  for  the  probability  of  detection  over 
time  of  stationary  targets  given  that  all  images  acquired 
are  non-overlapping. 


dot 


^one^S^tav^det* 

^int^si 
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Optimization  Non-Cued  Detection  —  Stationary  Targets 

Equation  37  gives  the  probability  of  non-cued  detection 
of  stationary  targets  for  individual  sensors.  To  identify 
the  best  sensor  or  sensors  to  use  in  a  given  situation 
requires  an  optimization  methodology. 

For  convenience,  equation  37  can  be  rewritten  in  the 
general  form-. 


^dot  *^im*^dotl 


Where 


^one^det* 


dot  1 


A, 


int 
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If  there  is  more  than  one  sensor  available 
expression  for  overall  probability  of  detection 
linear  combination  of  each  sensor's  individual 
the  number  of  images  acquired  From  that  sensor, 
combination  is  given  in  equation  42. 


,  the 
is  the 

•^dotl 

This 
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^dsys 

Where  f’dsus  ”  overall  probability  of  detection 
n  ^  “  number  of  sensors  available 

C(Slim)i  ”  number  of  images  acquired  from 
sensor  i 

CPdoti^i  ”  probability  of  detection  in  a  single 
image  from  sensor  i 

Equation  45  is  counter-intuitive.  However,  this  result 
is  important  and  must  be  accepted  for  the  subsequent 
analysis  to  have  any  merit.  By  way  of  example,  consider  the 
situation  of  a  checkerboard  with  a  single  checker  placed  on 
a  square  at  random.  A  player  is  blindfolded  and  given  the 
task  of  finding  the  checker.  He  searches  for  the  checker  by 
asking  that  a  third  party  look  at  a  number  of  specific 
squares  and  tell  him  if  the  checker  is  in  any  of  the 
squares.  The  restriction  on  the  number  of  squares  a  player 
may  choose  is  that  each  request  costs  him  money  and  he  has 
only  a  restricted  budget  available.  The  last  rule  is  that 
the  checker  is  not  moved  during  the  game  so  that  revisiting 
a  square  does  the  player  no  good.  The  game  begins. 

The  task  of  finding  a  stationary  target  using  satellite 
imagery  of  an  area  of  interest  is  the  checkerboard  game 
scaled  up  to  the  dimensions  of  the  area  of  interest.  The 
total  areas  of  interest  are  the  S4  squares  of  the 
checkerboard  and  the  7,000,000  square  kilometers  in  Northern 
Canada.  Picking  one  square  is  analagous  to  contracting  for 


one  image.  The  probability  of  finding  the  checker  is  one 
Cl)  given  the  player  request  that  the  square  holding  the 
checker  be  looked  at.  The  probability  that  the  player  will 
request  the  correct  square  is  N/64  where  N  is  the  number  of 
requests  made  at  the  start  of  the  game.  Thus  the 
probability  of  finding  the  checker  increases  linearly  with 
the  number  of  trials  taken.  Similarly  the  probability  of 
finding  a  stationary  target  in  Northern  Canada  increases 
linearly  with  the  number  of  images  contracted.  In  each  case 
the  probability  of  success  is  maximized  at  the  point  that 
ail  64  squares  are  chosen  or  all  the  territory  is  imaged. 
Beyond  this  point  both  the  analogy  and  the  analysis  break 
down.  Therefore,  for  the  linear  analysis  which  follows  to 
be  valid,  the  restriction  of  non-overlapping  images  must  be 
honored . 

The  first  constraint  on  is  that  there  is  an  upper 

limit  to  the  number  of  images  that  can  be  acquired  as 
constrained  by  the  available  budget.  That  is: 

Bj  C43) 

1-1 

Where  =  cost  of  one  image  from  sensor  i 

The  second  constraint  limits  the  number  of  images  that 
can  be  provided  by  the  sensors.  Since  the  sensors  are 
restricted  to  providing  non-overlapping  images,  once  the 
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entire  area  of  interest  has  been  imaged  further  acquisition 
is  terminated.  This  constraint  takes  the  form: 

Where  ^^one^i  “  image  taken  by  sensor  i 

The  third  constraint  is  actually  a  family  of 
constraints;  those  that  restrict  the  maximum  number  of 
images  that  can  be  acquired  from  a  specific  sensor  in  a' 
specific  period  of  time.  These  constraints  take  the  form: 

0  <  CN,„).  <  Ki  1:45 ) 

im  1  1 

Where  =  maximum  number  of  images  that  can  be 

acquired  from  sensor  i  in  the  period 

The  fourth  constraint  requires  explanation.  The  goal 
of  non-cued  detection  of  stationary  targets  is  to  achieve  a 
uniform  probability  of  detecting  targets  throughout  the 
entire  area  of  interest.  However,  the  area  0^"  interest  is 
not  uniform.  The  non-uniformity  occurs  since  there  are 
regions  within  the  area  of  interest  that  are  more  easily 
imaged  than  others  because  of  darkness  and  a  difference  in 
the  probability  of  cloud  cover.  As  the  model  stands  now  any 
optimization  methodology  would  select  areas  to  irnag-  that 
would  return  the  highest  individual  Pjjgt*  expense  of 
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coverage  of  lower  areas.  To  enforce  a  uniform 

probability  of  detection  throughout  the  area  of  interest, 
areas  of  common  obscuration  phenomena  must  be  identified, 
and  constraints  added  which  specify  that  each  unique  area 
must  achieve  the  same  P^sys" 

Having  built  the  expression  to  be  maximized,  and 
having  identified  the  series  of  constraints  which  must  be 
met,  the  task  of  optimizing  the  selection  of  sensors 
remains.  One  method  which  will  solve  the  combination  of 
equations  above  is  linear  programming  C49:7E),  This  can  be 
accomplished  using  equation  42  as  the  objective  function  to 
be  maximized  and  equations  43,  44,  and  45  plus  those 
specifying  equal  P^jgyg  far  each  unique  area  as  the 
constraints.  This  yields  the  maximum  probability  of 
detection  given  a  specified  budget.  A  second  option  is  to 
specify  a  minimum  acceptable  P^jgyg  for  equation  42  and  write 
the  linear  program  using  equation  43  as  the  objective 
function  to  be  minimized.  This  would  solve  for  the  minimum 
cost  to  achieve  a  constrained  probability  of  detection. 

A  problem  arises  in  that  the  value  for  P^jq^i  changes  as 
a  function  of  the  time  of  year,  type  of  target,  and  location 
of  the  target.  Because  of  this,  care  must  be  taken  to 
explicitly  define  what  is  to  be  maximized.  In  its  largest 
scale  the  analysis  could  produce  a  maximum  probability  of 
detection  of  all  target  types  located  in  all  possible 


locations  throughout  the  entire  year.  Analysis  at  this 
scale  is  impractical  if  the  numbers  of  sensors,  target 
types,  locations,  and  distinct  times  of  year  are  large.  In 
such  a  situation  either  simplifying  assumptions  are 
required,  or  subsets  of  the  large  linear  programme  can  be 
analyzed . 

In  using  these  algorithms,  one  final  caveat  is 
appropriate.  The  situation  could  arise  that  complete 
imaging  of  an  area  results  in  a  probability  of  target 
detection  belou  that  specified  by  operational  requirements. 
For  the  system  to  yield  a  higher  probability  of  detection, 
multiple  imaging  of  the  area  of  interest  ujould  be  required. 
That  is,  sensors  uiould  have  to  be  tasked  to  image  areas  more 
than  once.  This  would  have  the  effect  of  increasing  the 
probability  of  detection  for  those  areas  that  were  imaged 
more  than  once.  Hultipie  imaging  of  areas  is  discussed  in 
the  next  section  of  this  chapter. 

Non-Cued  Detection  Algorithm  —  Hoving  Targets 

If  a  target  is  moving  in  relation  to  the  background, 
acquiring  an  image  of  a  specific  ground  cell  at  time  tl  does 
not  absolve  the  detection  system  from  having  to  re-image 
that  area  at  any  other  time  in  order  to  see  what  is  there. 
The  probability  of  detecting  the  presence  of  a  target  in  the 
area  of  interest  again  depends  upon  the  probability  that  the 
target  will  be  detected  in  a  single  image  and  upon  the 
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number  of  images  acquired.  Houiever,  in  contrast  to  the  case 
of  stationary  targets,  the  relationship  between  the  number 
of  images  acquired  and  the  overall  probability  of  detection 
is  not  linear. 

The  first  case  to  consider  is  the  situation  where  a 
single  sensor  is  tasked  to  provide  multiple  images  of  the 
area  of  interest.  Such  a  situation  would  occur  if  the 
sensor  performance  in  detecting  a  specific  target  was 
dominant  in  relation  to  other  available  sensors.  Defining  a 
detection  to  be  a  success,  the  probability  of  success 
Cdetection)  in  a  single  image  is  P^jq^I  *  Then  the 
probability  of  X  successes  in  N  trials  C images  acquired)  is 
binomially  distributed  C21;3B).  For  example,  the 
probability  of  at  least  one  success  as  a  function  of  the 
number  of  trials  using  a  probability  of  success  of  .1  is 
shown  in  Figure  15. 
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Figure  15.  Probability  of  Detection  vs  Number  of  Trials 

for  PCSuccessi= . 1 


The  second  case  to  consider  is  the  situation  where  a 
variety  of  sensors  are  used  to  provide  an  image  of  the  area 
of  interest.  One  method  to  calculate  the  combined 
probability  of  detection  is  by  buiding  a  probability  tree 
based  on  the  individual  sensor's  probability  of  target 
detection  for  a  single  image,  ^dotl  ■ 

Two  hypothetical  sensors  can  illustrate  this  procedure. 
I  Given  their  respective  P^otl®  probability 

tree  which  results  from  acquiring  one  image  from  each  sensor 
is  shown  in  Figure  16.  Because  the  goal  of  multiple  imaging 
I  is  to  achieve  an  acceptable  probability  of  at  least  one 

detection,  the  outcome  of  primary  importance  is  the  one 
which  shows  "no  detection"  throughout.  The  probu'-'i  1  i ty  of 
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at  least  one  detection  is  equal  to  one  Cl)  minus  the 
probability  of  no  detections  throughout. 


^d>l  *  ^  ^nodet 


C45) 


Where  P(j>l  “  probability  of  at  least  one  detection 
^nodet  P’^Qtiabiiity  of  no  detection  throughout 


The  probability  of  the  no  detection  throughout  outcome 
is  the  product  of  the  individual  probabilities  of  no 
detection  for  each  image  acquired.  This  is  because,  as 
shown  in  the  figure,  there  is  always  one  all  negative 
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outcome  irrespective  of  the  number  of  images  acquired.  All 
other  combinations  contain  at  least  one  detection.  Each 


sensor’s  individual  probability  of  no  detection  is  one  Cl) 
minus  the  probability  of  detection  for  a  single  image, 
Pdoti’  Equation  47  provides  the  expression  for  the 
probability  of  at  least  one  detection  in  a  series  of  images 
acquired  from  a  variety  of  sensors. 
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Ulhere  ^dotll  ~  probability  of  detection  For  a  single 

image  From  sensor  1 

^dotlE  probability  of  detection  For  a  single 
image  from  sensor  2 

Pdotin  “  probability  of  detection  for  a  single 
image  from  sensor  n 

xl  =  number  of  images  acquired  From  sensor 

1 

xE  =  number  of  images  acquired  from  sensor 

2 

xn  =  number  of  images  acquired  from  sensor 

n 


This  relationship  holds  For  both  cases  discussed  above 
since  the  selection  of  multiple  images  from  a  single  sensor 
is  a  subset  of  the  general  case. 

□ptimization  Non-Cued  Detection  —  novino  Target 

As  is  the  case  for  stationary  target  detection 
optimization,  moving  target  detection  optimization  can  be 
approached  from  tujo  perspectives.  The  goal  is  either  to 


maximize  the  subject  to  the  constraints  of  budget  and 

image  availability,  or  co  minimize  cost  subject  to  the 
constraints  of  an  acceptable  and  image  availability. 

However,  since  equation  47  is  non-linear,  linear  programming 
cannot  be  used  directly  to  solve  this  model. 

The  probability  of  at  least  one  detection  can  be 
maximized  through  a  variety  of  non-linear  programming 
techniques.  In  its  general  form  equation  47  is  the 
abjective  function  to  be  maximized,  and  equations  43,  and  45 
make  up  the  constraints. 

To  minimize  the  cost  involved  in  achieving  an 
acceptable  again  non-linear  programming  techniques 

could  be  usea .  In  this  case  equation  43  becomes  the 
abjective  function  to  be  minimized,  and  equations  45  and  47 
became  the  constraints. 

In  either  case  the  non-uniformity  of  the  area  of 
interest  must  be  accounted  for.  This  will  generate 
additional  constraints  to  enforce  an  uniform  P^>^  for  each 
unique  region  within  the  area  of  interest. 

Cued  Detection  Algorithm  —  Stationary  Targets 

Cued  detection  of  stationary  targets  is  very  similar  to 
non-cued  detection  of  moving  targets.  Again  the  probability 
of  detecting  the  presence  of  a  target  depends  on  the 
probability  that  the  target  will  be  detected  in  a  single 


85 


image  and  the  number  of  images  acquired.  Again  the  logic  of 
multiple  imaging  of  an  area  applies.  However,  a  major 
difference  occurs  in  the  size  of  the  area  of  interest.  In 
cued  detection  of  a  stationary  target  a  single  location  is 
imaged.  This  effectively  reduces  the  area  of  interest  to 
the  area  of  a  single  image  from  the  sensor.  Thus  the 
probability  of  detection  by  an  individual  sensor  which 
applies  is  the  raw  probability  of  detection  amended  to 
include  appropriate  obscuring  factors,  or 
equation  For  calculating  the  probability  of  at  least  one 
detection  becomes; 


Where 


det*l , S ,  .  .  . , n 


probability  of  detection  in 
one  image  from  sensor 
1,2 . n 


Optimization  Cued  Detection  —  Stationary  Targets 

As  in  non-cued  detection  of  moving  targets,  the 
optimization  methodology  used  to  optimize  cued  detection  of 
stationary  targets  depends  on  whether  the  P(j>2^  is  to  be 
maximized  or  the  cost  is  to  be  minimized.  Again  a  variety 
of  techniques  is  available  to  accomplish  this  optimization. 
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Cued  Detection  Algorithm  —  noving  Targets 

Cued  detection  of  moving  targets  is  an  alias  for 
tracking.  Therefore,  the  algorithm  developed  in  this 
section  to  assess  a  sensor’s  capability  to  detect  moving 
targets  uihen  cued  is  actually  a  measure  of  the  sensor’s 
ability  to  track  a  target.  Similarly,  the  probability  that 
a  sensor  mill  detect  the  target  is  equivalent  to  its 
probability  of  tracking  the  target. 

There  are  tuio  broad  tracking  categories.  First,  a 
sensor  may  be  tasked  to  establish  a  second  contact  with  a 
target.  In  this  case,  the  first  contact  has  occurred  and 
constitutes  the  cueing  information.  Second,  a  sensor  may  be 
tasked  to  locate  a  target  for  which  at  least  two  contacts 
have  been  established.  In  this  case,  the  task  may  be 
simplified  since  some  prediction  of  the  target  location  is 
possible . 

Second  Contact .  Given  that  the  first  contact,  or  cue 
information,  is  accurate,  the  probability  that  a  system  will 
detect  a  target  depends  on  the  sensor’s  single  image 
orobability  of  detection,  the  size  of  the  image,  the 
position  of  the  first  contact,  the  velocity  of  the  target, 
and  upon  the  responsiveness  of  the  sensor.  Consistent  with 
an  approach  of  simple  first,  complex  second,  the  analysis 
begins  with  consideration  of  Fixed  sensors. 
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Sensors  Fixed  at  Nadir.  The  general  situation  of 
sensor  suiath  and  target  area  overlap  is  shown  in  Figure  17. 
In  the  figure,  X  is  the  position  of  the  target  at  time  tj^,  Z 
is  the  subsatellite  point  at  time  tg  when  the  satellite 
crosses  the  latitude  of  initial  target  detection  and  the 
sensor  is  u/ithin  range  of  the  possible  target  area. 


Figure  17.  Sensor  Swath  Overlap  of  Possible  Target 
Area 


The  possible  target  area  is  enclosed  in  the  circle  of  radius 

r,  where  r  is  equal  to  the  distance  that  the  target  could 

have  moved  if  it  were  travelling  at  top  speed  For  the  period  I 

tg  -  tj^.  The  diagonal  swath  is  the  territory  that  the 

sensor  could  image  on  the  current  pass. 

If  the  ground  swath  overlaps  the  possible  target  area,  i 

and  if  the  target  is  equally  likely  to  be  anywhere  in  the 
possible  target  area,  then  the  probability  that  the  sensor 
will  detect  the  target  on  this  pass  is  given  in  equation  43.  j 
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^track  ”  ^trkcon^det*  CH3) 

Where  f’track  “  Pr'atiability  sensor  mill  track  target 

on  this  pass 

Ptrkcon  ”  ratio  of  the  area  imaged  which 

could  contain  the  target  to  the  area 
which  could  contain  the  target 


In  the  limiting  case  the  sensor's  swath  just  touches 
the  area  in  which  the  target  could  be.  This  is  shown  in 


Figure  IB.  Sensor  Swath  Edge  Tangential  to  Possible 

Target  Area 

Figure, 18.  At  time  t^,  when  the  sensor  is  overhead  point 
Z3,  the  distance  between  Z3  and  X  is  given  in  equation  50. 


CZ3  -  X)  =  (t^-t,)s  +  5,,/E  C50) 

J  i  w 

Where  S^  =  width  of  the  sensor’s  ground 
resolution  cell 
s  =  target  speed 


as 


Angle  Ca)  is  equal  to  IBO*  minus  the  local  crossing 
angle  Cq) .  The  local  crossing  angle  is  a  Function  of  the 
orbit  inclination  and  the  latitude  of  interest.  The 
formulas  for  determining  the  local  crossing  angle  are 
contained  in  Appendix  E.  The  distance  between  Z*  and  X  is: 


CZ*  -  X) 


Ct3-4)s  ^  S^/S 


sin  Ca) 


C51) 


This  distance  can  be  converted  to  the  equivalent  number 
of  degrees  of  longitude.  Assuming  a  spherical  Earth,  the 
circumference  of  a  specific  parallel  of  latitude  is  given  by 
the  expression; 


Cirj_  =  STfRgCOsCL) 


CSS) 


Where  Cirj_  “  circumference  at  latitude  L 

R  =  radius  of  the  Earth 

L  -  latitude 

The  number  of  degrees  of  longitude  equivalent  to  the 
distance  CZ*  -  X)  is  then.- 


/  = 


Z*  -  X 


Cir 


L 


C3B0'  ) 


C53) 


The  equivalent  degrees  of  longitude  is  used  to 
determine  whether  the  sensor’s  swath  and  the  passible  target 
area  overlap.  IF  at  time  tg  the  subsatellite  longitude  is 
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ulithin  the  threshold  boundaries  C±  )  of  the  initial  target 
longitude,  the  sensor  suiath  and  possible  target  area  uiill 
overlap.  The  relationship  between  the  threshold  longitude 
for  overlap  and  is  linear  and  is  shown  for  60*N 


I 

sj 

5 

f 

o 

I 

s 


giqpid  Tim  Ov«r«) 


9  Ovts 


*0  Krota 

_i - 

20  tenets 


Figure  19.  Threshold  Longitude  vs  Elapsed  Time 


latitude  and  the  Landsat  5  swath  in  Figure  19.  Curves 

plotted  for  Landsat  5  and  SPOT  1  are  provided  in  Appendix  E.  ^ 

Having  determined  that  there  will  be  overlap  at  time 
tg,  the  question  remains  as  to  how  much  area  is  overlapped. 

This  overlap  will  depend  on  the  size  of  the  sensor  swath  and  ^ 

the  size  of  the  possible  target  area.  The  general  situation 
of  overlap  is  shown  in  Figure  20.  Since  the  longitude  of 

I 
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Figure  SO.  Sensor  Suath  Overlapping  Possible 
Target  Area 


the  subsatellite  point  2*  is  known,  the  distance  CZ*  -  Xj 
can  be  determined  using: 

Cl*-1^)  Cir, 

(2*  -  XI  =  - ^  C54) 

360* 

Where  1*  =  longitude  of  Z* 

1^  =  longitude  of  initial  contact 

Once  CZ*  -  X)  is  known,  the  distance  Y  is; 

Y  =  CZ*  -  X)  sinCa)  CSS) 

□nee  Y  is  known  the  area  of  the  swath  can  be  determined 
geometrically.  Two  situations  can  exist.  First,  Y  may  be 
less  than  half  the  swath  width  of  the  sensor,  in  which  case 


-  -4 

I 


J 


ss 


the  suiath  stradles  the  initial  target  location.  Second,  Y 
may  be  greater  than  or  equal  to  half  the  swath  width  of  the 
sensor,  in  which  case  the  swath  falls  in  only  one  half  of 
the  circle. 


Figure  El.  Geometry  For  Swath  Overlapping  Center  of  ( 

Possible  Target  Area 

Figure  El  shows  the  geometry  which  results  when  Y  is  ] 

less  than  half  the  sensor’s  swath  width.  In  this  case  the 
area  of  the  circle  which  is  overlapped  by  the  swath  is  given 
in  equation  56.  I 
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^SWC  “  ^cir  ‘^XAB^  ^^XCD  ^XCD^ 


C5B) 


Uhere 


area  of  swath  overlap  in  circle 

area  of  the  circle 

area  of  sector  XAB 

area  of  triangle  XAB 

area  of  sector  XCD 

area  of  triangle  XCD 


The  angle  C0)  is  calculated  using; 


casC0)  = 


Y  - 


Where  ^d  *  ~  ^1 


The  angle  FXD,  </3)  is; 


C57) 


casC/3) 


3^/2  -  V 


C50) 


The  area  of  sector  XAB  is: 


^XAB  *  CTrt^^s^) 

360 


C59) 


The  area  of  triangle  XAB  is; 


Axab  °  Y)'^tanC0) 


(BO) 
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The  area  of  sector  XCD  is: 

20 

SxcD  “  —  ^  ^ 

360 

The  area  of  triangle  XCO  is: 

^XCD  ”  “  Y)^tanC0)  C6E) 

The  analysis  for  determining  the  area  of  overlap  uihen 
the  suath  overlaps  the  initial  contact  point  is  not  quite 
complete.  There  are  two  special  situations  that  require 
additional  consideration. 

First,  the  situation  could  arise  that  the  subsatellite 
point  of  the  sensor’s  next  crossing  of  the  initial  contact 
latitude  exactly  corresponds  to  the  initial  detection  point 
This  is  shown  in  Figure  28.  In  this  case  the 
remains  at  unity  until  the  possible  target  area  exceeds 


Figure  82.  Geometry  for  Sub-Satellite  Point  Exactly 
Over  Initial  Target  Location 


35 


the  size  of  the  swath  width.  Therefore  when  equation  B3  is 
true  the  P^rkcon  assigned  a  value  of  1 . 


tdS  -  Sy/2  <  0 


CB3) 


Second,  in  the  situation  that  the  sensor’s  swath 
contains  the  initial  contact  point,  there  is  a  period  of 
time  when  the  possible  target  area  does  not  fully  extend 
through  the  sensor's  swath.  This  is  shown  in  Figure  23. 

For  any  t^j  less  than  the  time  required  to  travel  to  the  edge 


Figure  23.  Geometry  for  Swath  Overlapping  Initial 

Target  Location 


of  circle  1,  the  P 


trkcon 


must  be  1.  For  the  period  it  takes 


to  transit  from  circle  1  to  circle  2,  the  Formula  For 
overlap  is  given  in  equation  B4 . 


'^SUJC  "  '^cir  ‘-^XCD  '^XCD’’ 


C64) 


Figure  24  shoujs  the  geometry  uhich  results  when  Y  is 
greater  than  half  the  sensor’s  suiath  width.  In  this  case 
the  area  of  the  circle  which  is  overlapped  by  the  swath  is: 


'^SUC  "  ^XOC  “  ^XDC  '^^XAB  ~  '^XAB^ 


Where  *  area  of  sector  XDC 

^XDC  a^ea  of  triangle  XDC 
^XAB  ”  arsa  sector  XAB 
^XAB  *  area  of  triangle  XAB 


The  area  of  sector  XDC  is  then: 


^XDC 


aC30*  -  c) 


360* 


Cirtjj^s^) 


The  area  of  triangle  XDC  is: 

^XDC  "  CQsCc)t^sCY  -  Sy/S) 
The  angle  Cd3  is; 


cosCdD 


Y  ^  Sy/E 


tcjs 


The  area  of  sector  XAB  is: 


E  E 

Sxab  '  -  ‘'^^d  ®  ^ 

“  360*  ° 


The  area  of  triangle  XAB  is: 

Axab  *  Sy/S)t|^s  sinCd) 


C57) 


CGB) 


C63) 


C70) 


C71) 


Pointable  Sensors.  All  of  the  geometry  to  this 
point  in  the  analysis  is  used  to  calculate  the  overlap  of 
the  sensor  swath  and  target  area  for  a  sensor  which  looks 
only  directly  below  itself.  How  does  this  change  in  the 
case  of  satellites  that  can  point  their  sensors  at  targets 
which  are  offset  from  the  satellite  ground  trace'^ 


SB 


The  geometry  for  a  pointable  sensor’s  overlap  with  a 
target  area  is  shown  in  Figure  25.  The  first  difference  in 
this  analysis  as  compared  to  fixed  sensor  analysis  is  the 
value  of  the  distance  Y.  This  distance  is  again  calculated 
using  equations  55  and  56.  However  since  Y  is  equal  to  the 
distance  from  the  ground  trace  to  the  initial  target 
location,  in  the  threshold  case  it  is  the  distance  off  track 
the  sensor  is  pointing  C2-S*)  plus  half  the  swath  width  of 
the  sensor  plus  the  target  speed  times  the  elapsed  time  Cr). 


Figure  25.  Pointable  Sensor  Threshold  Overlap 


Again  the  swath  may  or  may  not  overlap  the  point  of 
initial  target  detection.  In  the  case  of  swath  overlap  with 
the  point  X,  the  distance  Y  will  be  less  than  the  distance 
C2-S*)  plus  half  of  the  swath  width  of  the  sensor.  This  is 


shown  In  Figure  SB.  In  the  figure  the  distance  dS  is  equal 
to  CZ-S*)  plus  half  the  swath  width  minus  Y. 

dS  -  CZ-S*)  +  Sy/a  -  Y  C7S) 

The  distance  dl  is  Y  minus  CZ-S*)  plus  half  the  swath  width. 
Using  these  distances  the  area  of  the  swath  overlap  with  the 


Figure  SB.  Geometry  for  Pointable  Sensor  Swath  Overlap 


If  the  sensor’s  swath  does  not  overlap  the  initial 
target  location,  the  swath  overlap  geometry  is  shown  in 
Figure  S4 ,  In  this  case  the  distance  XF  is  calculated  by 
subtracting  the  distance  CZ-S*)  and  half  the  swath  width 
from  the  distance  Y.  Thereafter  the  calculations  for 


overlap  follow  those  developed  in  equations  65  through  71. 


Conditional  Tracking  Probability  Calculation.  All 
of  the  geometry  shown  so  far  contributes  one  half  of  the 
information  needed  to  determine  the  conditional  tracking 
probability.  The  second  half  is  Just  the  total  possible 
target  area,  ^pQg •  It  is: 


pos 


C735 


Finally,  an  expression  for  the  probability  of  tracking, 
conditional  on  the  sensor  seeing  the  target,  can  be  written. 


^trkcon  (.7^:1 

The  Ptrkcon  sensor  specific.  Defining  the  sensor’s 
response  time,  t^,  to  be  the  elapsed  time  between  the 
initial  target  contact  and  the  time  of  sensor  crossing  of 
last  known  target  latitude,  the  P^rkcon  plotted 

against  sensor  response  time.  The  general  form  of  the  curve 
is  shown  in  Figure  E7,  which  is  made  using  the  Landsat  5 
swath,  an  initial  target  latitude  of  60*N,  and  target  speed 
of  5  Knots. 

Not  only  does  the  ^trkcon  ‘^screases  rapidly  as  the 
possible  target  area  increases,  but  also  more  images  are 
required  to  achieve  the  Ptrkcon'  relationship  between 

the  equivalent  number  of  images  required  to  achieve  the 
maximum  P^rkcon  response  time  is  shown  in  Figure  SQ . 
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Figure  S7 .  Conditianal  Tracking  Probability  verses 

Response  Time 

CTop  to  bottom  are  0,1,2,3,H,5*  Longitude  Displacement) 


Figure  27  shows  the  conditional  tracking  probability 
for  a  single  sensor,  a  single  target  speed,  and  a  single 
initial  target  detection  latitude.  There  are  many 
combinations  of  these  variables.  Nonetheless,  comparative 
plots  of  conditional  tracking  probabilities  can  be  used  as  a 
first  check  to  determine  if  a  given  sensor  could  provide  an 
acceptable  response.  For  this  reason  conditional  trncking 
curves  for  Landsat  5  and  SPOT  1  are  included  in  Appendix  F. 
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Figure  28.  Conditional  Tracking  Probability  and 
Equivalent  Number  of  Images  vs  Response  Time 


The  P^rkcon  calculated  above,  when  divided  by  the 
equivalent  number  of  images  required  to  achieve  it,  is  a 
measure  of  the  conditional  probability  of  tracking  in  a 
single  image,  or: 


P 


tim 


I 


p 

trkcon 

^qim 


C75) 


Where  ^tim  ”  conditional  tracking  probability  per 
equivalent  image 

E  im  ”  equivalent  number  of  images 


I 
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It  is  the  Ptim’  multiplied  by  the  Pj^g^*  that  gives 

the  probability  of  tracking  the  target  in  a  single  image,  or 

^trackl • 

^trackl  “  ^tim^det* 

Optimization  Tracking  Algorithm  —  Second  Contact 

The  algorithm  chosen  to  optimize  the  probability  of 
tracking  a  target  depends  on  whether  overlapping  or  non- 
overlapping  images  are  available. 

In  the  case  of  non-overlapping  images,  P^rack 
maximized  using  linear  programming.  The  abjective  function 
to  be  maximized  is  given  in  equation  77. 


^  ‘“'^im^i^tracki 
i-1 

Subject  to  the  constraint  set: 


C77) 


i  S  Bj  C7a) 

1*1 

0  <  x^  <  C73) 

Where  *  maximum  number  of  images  available  from 

sensor  i 

If  there  is  overlapping  coverage,  the  optimization 
becomes  untidy.  The  untidiness  occurs  since  the  objective 
function  must  account  for  areas  of  single  and  overlapping 
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coverage.  An  example  of  overlapping  swaths  is  shown  in 
Figure  S9 . 


Figure  29.  Overlapping  Swath  Coverage  Geometry 


In  the  figure,  Sensor  A  passes  overhead  at  time  tl, 
Sensor  B  passes  overhead  sometime  later  at  time  t2 .  The 
areas  of  single  and  double  coverage  depend  on  the  swath 
widths  of  the  sensors,  the  sensors’  inclinations,  and  the 
time  of  each  of  the  sensor’s  passages  overhead. 

The  areas  in  which  only  a  single  sensor  provides 
coverage  behave  linearly — that  is,  their  contribution  to 
overall  tracking  probability  provided  by  images  acquired  in 
these  areas  will  increase  linearly  as  a  function  of  the 
number  of  images.  Given  complete  overlap  of  multiple  imaged 


areas,  their  contribution  to 
linear.  However,  since  in  a 


the  overall  P^rack  non- 

given  situation  the  individual 


Ptrack  will  be  known,  and  since  only  one  image  From 


each  sensor  can  be  taken  of  an  area  on  the  ground  in  a 
single  pass,  the  combined  ^trackl+S  fixed  and  is: 


^trackl+2  “  ^  ^trackl  ^ ^"^trackS ^  ^ 


C80) 


Since  ^trackl+E  ®  fixed  value  it  also  contributes  to 
the  objective  function  linearly.  Therefore  optimization  for 
overlapping  coverage,  assuming  the  overlap  is  complete,  can 
be  accomplished  through  linear  programming.  The  form  of  the 
objective  function  to  be  maximized  is: 


^track  "  ^l^truckl  ^  ^2^track2  ^  ^l^-2^trackl+2 


cai) 


Where  g  “  number  of  non-overlaping  images  acquired 

’  from  sensors  1,2 
Xi^2  “  number  of  overlapping  images  taken 


The  constraint  set  is: 


XiCi 

+  XgCg  Xj^^gCCj^+Cg)  <  Bj 

C82) 

Xi  < 

^qiml'^^’^®®  Single 

Coverage 

1/Total 

Swath 

1)  C83) 

Xg  < 

^qima'-^^®®  Single 

Coverage 

2/Total 

Suiath 

2) 

^1+2  - 

E„;„iCArea  Double 
qiml 

Coverage/Total  Suath  Double) 

prior 

true 


In  practice  single  area  coverage  will  be  exhausted 
to  the  acquisition  of  overlapping  images.  This  is 
since  a  combined  image  yields  a  ^trackl+E  iJJhich  is 
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lower  than  the  sum  of  its  components  in  all  cases,  but  which 
costs  the  same  as  the  individual  images. 

Tracking  Algorithm  —  Follow-on  Contacts 

Once  a  target’s  position  is  plotted  for  two  points  in 
time,  its  location  at  any  subsequent  time  can  be  predicted. 
Assuming  the  target’s  velocity  remains  constant,  the  task  of 
establishing  follow-on  contacts  reduces  to  cued  detection 
and  the  algorithms  developed  for  cued  detection  of 
stationary  targets  can  be  applied.  The  only  glitch  that 
remains  is  determining  if  and  when  the  sensor  will  overfly 
the  target. 

The  question  of  simultaneous  target  location  and  sensor 
swath  overlap  can  be  broken  into  two  constituent  parts;  how 
soon  will  the  target  penetrate  a  specific  ground  swath,  and 
how  long  will  the  target  take  to  traverse  the  swath.  Uith 
these  two  pieces  of  information,  it  is  possible  to  determine 
whether  the  satellite  will  overfly  the  target  on  a  given 
pass.  That  is,  a  window  of  opportunity  for  imaging  during  a 
specific  pass  can  be  determined. 

The  first  piece  of  information  to  determine  is  how  soon 
a  target  will  penetrate  a  specific  ground  swath.  Figure  30 
illustrates  the  geometry.  The  shortest  path  to  the  swath  is 
perpendicular  to  the  inclination  of  the  sensor’s  orbit. 

This  distance  is  given  by  equation  04. 
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Figure  30.  Su/ath  Penetration  Geometry 

-  sinCa)j^*CF)  -  5^/2  (84) 

Where  Dg^  *  shortest  path 

F  -  #  km  per  degree  longitude 
at  latitude 

The  distance  to  the  suiath  increases  as  the  target 
heading  varies  From  the  perpendicular.  The  distance  to  the 
suiath  For  any  angle  other  than  the  perpendicular  is 
determined  by  dividing  the  shortest  distance  by  the  cosine 
oF  the  angle  between  the  target  heading  CH)  and  the 
perpendicular.  ThereFore,  the  distance  to  the  suiath  is; 


A  negative  result  For  D  indicates  that  the  target  is 
heading  auay  from  the  suiath.  This  would  occur  if  the  target 
were  heading  East  and  the  sensor’s  swath  under  consideration 
were  UJest  of  the  target’s  last  known  position. 

Once  the  distance  from  the  target’s  last  known  position 
to  the  swath  is  calculated,  the  time  it  takes  to  transit 
this  distance  is  Just  the  distance  divided  by  the  target’s 
speed.  This  provides  the  first  piece  of  information  needed 
for  the  imaging  window. 

The  second  piece  of  information  needed  for  the  imaging 
window  is  the  time  it  will  take  the  target  to  transit  a 
swath.  The  geometry  of  swath  transit  is  also  shown  in 
Figure  30.  The  difference  is  that  the  shortest  distance  in 
this  case  is  the  swath  width  itself.  Using  the  swath  width 
as  the  numerator  in  equation  B5,  swath  transit  times  for  a 
given  target  speed  can  be  worked  out. 

Optimization  Tracking  Algorithm  —  Follow-on  Contacts 

If  the  above  analysis  yields  several  imaging 
opportunities,  the  task  remains  to  optimize  the  selection  of 
sensors.  A  complexity  arises  since  the  value  of  the  target 
being  tracked  will  dictate  the  frequency  with  which  track 
updates  will  be  required.  At  one  extreme,  a  low  value 
target  may  not  require  tracking.  At  the  other  extreme  a 
high  value  target  may  require  continuous  track  updates. 


109 


Once  again  two  distinct  approaches  can  be  taken  to 


optimize  selection  of  sensors.  If  the  goal  is  to  maximize 
the  probability  of  at  least  one  detection  of  the  target,  the 
objective  function  becomes: 


d>l 


1  -  P 


nodet 


cas) 


Where  ^nodet  "  ^det*ii^  v, 

Pdet*i 1  “  single  image  detection  probability 
for  sensor  i  on  pass  J 


And  the  constraint  set  is; 

(B7) 

0  <  CBSJ 

Where  -  maximum  number  of  images  available 

from  sensor  i 

□nee  again  there  are  a  variety  of  non-linear 
programming  alternatives  which  can  optimize  the  objective 
function  subject  to  the  constraint  set. 

Summaru 

In  summary,  this  chapter  has  answered  the  second 
question.  That  is,  the  selection  of  sensors  from  those 
available  is  accomplished  by  optimizing  the  performance  of 
the  system  of  sensors.  The  system  performances  are  measured 
in  terms  of  its  probabilities  of  detecting  or  tracking 


targets.  Therefore,  the  selection  of  imagery  from  a 
specific  sensor  depends  on  houi  large  its  contribution  to  the 
overall  system’s  probabilities  of  detection  and  tracking  is 
in  comparison  to  the  options  which  are  available. 


Ill 


General 


The  third  question  to  resolve  is  how  valid  are  the 
algorithms  which  have  been  developed  in  Chapter  V.?  As  with 
all  models,  the  quality  of  the  output  will  depend  on  the 
quality  of  the  input.  Nonetheless  there  is  value  in 
developing  the  algorithms  for  specific  situations  since  this 
will  afford  face  validity  to  the  methodology  itself.  Thus, 
this  analysis  chapter  presents  an  example  formulation  in 
each  mission  area. 

Example  Analuses 

At  the  start  of  an  analysis  the  analyst  must  specify 
the  mission,  the  targetCs)  of  interest  and  the  time  Frame 
for  consideration.  Based  on  the  target  of  interest  and  the 
mission,  the  area  of  interest  can  be  defined.  Based  on  this 
area  and  the  time  frame  of  the  analysis,  all  possible 
backgrounds  can  be  described  in  spatial,  spectral,  and 
temporal  terms.  Similarly,  all  targets  can  be  described  in 
these  three  domains.  Next,  the  possible  sensors  can  be 
selected  based  on  the  time  frame  of  the  analysis. 

Using  this  database,  the  analysis  begins.  The  first 
question  to  be  answered  is  whether  the  sensor  will  detect 
the  target  against  the  background  under  ideal  conditions. 


HE 


This  requires  analysis  of  each  sensor  with  every  possible 
target  and  background  combination.  For  visible  wavelength 
sensors,  equation  10  is  used.  If  the  actual  target  area  is 
larger  than  the  threshold  target  area,  the  target  can  be 
detected  by  the  sensor.  For  quick  reference  Appendix  C 
contains  threshold  detection  curves  for  various  sensors. 
Similarly,  equation  S4  is  used  to  determine  if  a  thermal 
wavelength  sensor  can  detect  a  target  under  ideal 
conditions . 

Comparing  actual  target  areas  to  threshold  target  areas 
returns  a  of  one  Cl)  in  the  case  of  detection,  or  zero 

CO)  in  the  case  of  no  detection.  At  this  point  in  the 
analysis,  all  target,  background,  and  sensor  combinations 
which  return  a  of  0  are  discarded. 

To  this  point  the  steps  required  are  common  for  all  the 
missions.  They  form  a  checklist  of  actions  which  begin 
every  analysis.  In  contrast,  subsequent  analysis  depends  on 
the  mission  under  consideration.  To  continue  this  example 
it  is  assumed  that  the  mission  is  non-cued  detection  of 
stationary  targets. 

The  next  step  in  the  analysis  is  to  determine  the  Pj^g^ 
considering  detector  array  geometry.  This  is  accomplished 
using  equations  27,  20,  and  29.  It  is  this  value  of  Pdet 
which  is  carried  into  the  rest  of  the  analysis. 


The  analysis  continues  with  the  partitioning  of  the 
target  area  into  areas  of  common  obscuration.  Once  this  is 
done,  the  each  combination  of  sensor,  target, 

background,  and  location  can  be  determined  using  equations 
35  and  36.  If  the  irradiation  factor  for  a  combination  is 
zero  CO),  this  combination  is  discarded. 

Having  determined  the  P^jg^*  values  for  the  surviving 
combinations,  the  next  task  in  the  analysis  is  to  work  out 
the  corresponding  P^jq^I  '^^lue  for  each,  using  equation  Htl . 
The  area  of  interest  here  is  the  area  of  common  obscuration. 

At  this  point  the  linear  program  can  be  written.  In 
this  example,  the  objective  function  to  be  maximized  is 
equation  Hk2,  and  the  constraint  set  is  furnished  by 
equations  H3,  44,  and  45,  and  by  specifying  each  area  must 
return  the  same  P^jgyg-  Solving  the  linear  program  gives  the 
Pdgys  achievable  for  the  budget  specified. 

Appendix  H  provides  checklists  for  analysis  start  and 
for  each  mission. 

The  remainder  of  this  chapter  contains  example  analyses 
using  the  representative  target  set  described  in  Appendix  D. 
Each  mission  area  is  analyzed.  Sources  of  data  are  enclosed 
in  parentheses. 
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First  Analusis:  Non-Cued  Detection  oF  Stationaru  Targets 
START  —  Checklist  1 

1.  Mission:  Non-cued  detection  of  stationary  targets 
Target:  Survey  Parties 
Time  frame:  Month  of  January 
S.  Area  of  interest;  60*UI  -  I'lS’W,  B0*N1  -  B3*N 

3.  Possible  backgroundCs) :  Snow  CRBf:lB) 

4.  Background  characteristics; 

Spatial  -  Elevation  <  1000  meters  CRefilE) 

Spectral  -  Reflectance  0.3  ClB:7j 

-  Emissivity  (10-12  jjm)  .3B5  C1B:7j 

Temporal  -  No  change  in  time  frame 

5.  Target  characteristics;  (Appendix  D) 

Spatial  -  Vehicles  up  to  50  m^ 

Spectral  -  Reflectance  0.12 

-  Emissivity  .92,  1/4  surface  area  2*K 
warmer  than  ambient 

Temporal  -  No  change  in  time  Frame 
B.  Sensors  available;  NOAA  AVHRR 

Landsat  5  MSS,  TM 
SPOT  1  HRV 

* 

Sojuzkarta 

*  Not  considered  due  to  the  delay  in  receipt  of  the  image 
7.  Pjjg^  calculation; 

Visible  -  (R^  -  Rj^)  =  .78 
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Table  3.  Visible  Sensor  First  Analysis 


Sensor 

AtiCm^) 

^det 

AVHRR 

1000  CFig. 

C17) 

0 

nss 

67  CFig. 

C7) 

0 

TM 

10  CFig. 

C5) 

1 

HRV 

1  CFig. 

Cl) 

1 

Thermal  -  Target 

exitance  28. 

8  W/m^ 

Table  H. 

Thermal  Sensor  First 

Analysis 

Sensor 

A^iCm^) 

^det 

Th  Band  6 

1  CFig. 

GH) 

1 

Surviving  options 

:  Landsat  Ttl 

Band  1 , 

3  or  H 

Landsat  Tfl  Band  B 
SPOT  HRV 


S.  End  checklist  1 


Ctia£kiiat_£i._NQnrCuei_Qetect.iQD_Qf_StatiQnar!y_Iar:aets 
1.  Landsat  5  and  SPOT  1  detector  arrays  are  staggered,  and 
there  are  effectively  no  gaps  in  the  detector  array. 

Thus  if  remains  at  one  Cl).  For  all 
surviving  sensors  this  is  true. 


! 


I 


■ 
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2.  Partition  area  of  interest:  The  two  considerations  are 
darkness  and  cioud  cover  since  both  thermal  and  visible 
sensors  have  survived  to  this  point. 

Darkness  - 

A.  Period  of  satellites  -  Landsat  5  is  30.9  minutes 

-  SPOT  1  is  101.5  minutes 

B.  Using  equations  32,  33  and  the  local  sidereal  times  of 

sunrise  and  sunset,  the  time  each  sensor  passes 
overhead  a  given  latitude  can  be  worked  out  and  a 
value  for  determined.  Table  5  contains  these 

values . 

Table  5.  Determination 


Latitude 

Landsat  Tj^^^ 

^fac 

SPOT  Ti3, 

^fac 

0 

0945 

1030 

60 

0320 

1 

1013 

1 

62 

0320 

1 

1012 

1 

64 

0327 

1 

1012 

1 

■ 

66 

0327 

0 

1011 

1 

60 

0326 

0 

1011 

1 

70 

0926 

0 

1010 

0 

72 

0325 

0 

1010 

0 

• 

Cloud  cover  - 

Cloud  cover  data  is  taken  from  Reference  10. 
Superimposing  the  areas  of  darkness  and  cloud  cover  yields 
six  distinct  areas  of  common  obscuration.  These  are  shown 
in  F igure  31 . 
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Figure  31.  First  Analysis  Areas  of  Common  Obscuration 


3/4/5.  Pjiet*  >^dotl  calculation: 

A.  Approximate  areas  of  each  of  these  regions  are  shown 
in  Table  6. 


Table  6.  Areas  of  Regions  of  Common  Obscuration 

First  Analysis 


Region 


B 

C 

D 

E 

F 

G 


005450 

1078040 

3017460 

631300 

845500 

6E7450 


i 

B.  Table  7  contains  the  values  for  Pjjgt*  ^dotl 
each  region. 


I 
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Table  7. 

F  irst 

Analysis 

^det*  ^dotl 

Region 

P 

cloud 

Sensor 

^det*  ^dotl 

B 

.6 

LT 

.4  1.0BE-3 

HRV 

.4  l.SBE-4 

C 

.6 

LT 

.4  l.BBE-3 

LV 

.4  l.BBE-3 

HRV 

.4  1.9BE-4 

D 

.5 

LT 

.5  2.3BE-3 

E 

.4 

LT 

.6  2.B3E-3 

F 

.5 

LT 

.5  2.36E-3 

HRV 

.5  2.4BE-4 

G 

.5 

LT 

.5  2.36E-3 

LV 

.5  2.36E-3 

HRV 

.5  2.4BE-4 

Legend  :  LT  - 

Landsat  5  Band  B 

LV  - 

Landsat  5  Band  3 

HRV  ■ 

-  SPOT  1 

6.  The  goal  of  this  analysis  is  to  maximize  the  system  wide 
probability  of  detection  given  a  total  of  51,000,000.00  for 
the  time  frame  under  analysis. 

7.  A.  Since  Landsat  5  Band  6  and  Band  3  sensors  return  the 
same  probability  of  detection  and  cost  the  same,  where  both 
are  available  only  Band  3  is  carried  forward  for  analysis. 
Using  equation  42  as  the  template  the  objective  function  is 
given  in  equation  09. 
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CB9) 


-  1.B0E-3  XI  +  l.aaE-3  xa  +  E.36E-3  X3 
^  +  a.B3E-3  X4  +  a.3BE-3  X5  +  E.36E-3  X6 

+  1.9BE-4  X7  +  1.90E-4  XB  +  E.4aE-4  X9 
+  E.4B  E-4  XIO 


Where  XI 

- 

# 

images 

taken 

by 

Landsat 

in  region 

B 

xa 

- 

images 

taken 

by 

Landsat 

in  region 

C 

X3 

- 

# 

Images 

taken 

by 

Landsat 

in  region 

D 

X4 

# 

images 

taken 

by 

Landsat 

in  region 

E 

X5 

- 

# 

images 

taken 

by 

Landsat 

in  region 

F 

XB 

- 

# 

images 

taken 

by 

Landsat 

in  region 

G 

X7 

- 

images 

taken 

by 

SPOT  in 

region  B 

XB 

- 

# 

images 

taken 

by 

SPOT  in 

region  C 

X9 

- 

# 

images 

taken 

by 

SPOT  in 

region  F 

XIO 

# 

images 

taken 

by 

SPOT  in 

region  G 

7.  B.  The  cast  of  a  Landsat  Tfl  image  is  $4,600.00.  The  cast 
of  a  9P0T  image  is  $1000.00  (43:357).  Using  equation  43  as 
the  template  the  cost  constraint  is: 


4600  XI  4600  XS  4600  X3  +  4600  X4 
+  4600  X5  +  4600  X6  +  1000  X7  +  1000  XB 

+  1000  X9  +  1000  XIO  <  1000000  (90) 


7.  C.  Using  equation  44  as  the  template,  the  restriction  on 
overlapping  images  is: 


34aE5 

XI 

+ 

3600 

X7  < 

005450 

(91 ) 

34EE5 

XE 

+ 

3600 

VI 

CO 

X 

1E70O4O 

C9E) 

34EE5 

X3 

< 

3017460 

(93) 

34aE5 

X4 

< 

691300 

(94) 

34aS5 

X5 

3600 

X9  < 

045500 

(95) 

34EE5 

X6 

3600 

XIO 

<  6E7450 

(96) 

D.  The  final  set 

of 

constraints 

restrict  the 

probabilities  of  detection  within  each  region  to  be  within  a 


ISO 


specified  range  of  each  other.  This  prevents  domination  of 
one  region  of  clear  uieather  at  the  expense  of  coverage  in 


harder-to-image  regions. 

To 

do  this 

the  area 

of 

the 

region 

is  used  to  specify  a  local 

Pdotl 

this  figure 

is 

used  to 

balance  the 

regions . 

1 .7E-2 

XI 

+  1.0E-3 

X7 

- 

1 .07E-2 

-  1.13E- 

3  < 

■ 

1 

C37) 

1 .7E-e 

XI 

+  1.8E-3 

X7 

- 

1 .07E-2 

-  1.13E- 

3  > 

1 

C3B) 

1 .7E-2 

XI 

+  l.aE-3 

X7 

- 

5.67E-3 

X3 

< 

• 

1 

C33) 

1 .7E-E 

XI 

+  1.8E-3 

X7 

- 

5.e7E-3 

X3 

> 

1 

ClOO) 

1 .7E-2 

XI 

+  1.8E-3 

X7 

- 

2.37E-4 

X'i 

< 

• 

1 

ClOl) 

1 .7E-a 

XI 

1.8E-3 

X7 

- 

2.37E-4 

X4 

> 

1 

C102) 

1 .7E-2  XI  + 

1 . 

aE-3  X7  - 

2. 

02E-2  X5  - 

2. 13E-3 

X3  < 

• 

1 

C103) 

1.7E-2  XI  + 

1 . 

aE-3  xa  - 

2.02E-2  X5  - 

2. 13E-3 

XIO 

>  - 

• .  1 

C104) 

1.7E-B  XI  + 

1 . 

aE-3  X7  - 

2. 

73E-2  X6  - 

2.07E-3 

XIO 

< 

.1 

C105) 

1.7E-a  XI  + 

1  , 

aE-3  X7  - 

2. 

73E-2  X6  - 

2.B7E-3 

XIO 

>  - 

-.1 

ClOB) 

3.  This  linear  program  was  salved. 


Second  Analysis :  Non-Cued  Detection  of  [loving  Targets 


START  —  Checklist  1 

1.  Mission:  Non-cued  detection  of  moving  targets 
Target:  Ship 

Time  frame:  Month  of  July 

2.  Area  of  interest:  Arctic  uateriuays 

Approximately  3,000,000  km^  C1E:E3) 

3.  Possible  backgroundCs) :  Water 

4.  Background  characteristics: 

Spatial  -  Sea  level 

Spectral  -  Reflectance  0.03  C1B:7) 

-  Emissivity  CIO-IS  ym)  0.33  C1B:7) 

Temporal  -  No  change  in  time  frame 

5.  Target  characteristics:  (Appendix  D) 

Spatial  -  lEl  x  15  meters 
Spectral  -  Reflectance  0.85 

-  Emissivity  0.74,  1/10  surface  area  E’K 
warmer  than  ambient 

Temporal  -  In  area  of  interest  entire  month 
B.  Sensors  available;  NOAA  AVHRR 

Landsat  5  MSS,  TM 
SPOT  1  HRV 

7.  Pjjg^  calculation: 

Visible  -  CR^.  -  R^^l  =  0.7B 

lEE 


Table  8.  Visible  Sensor  Second  Analysis 


Sensor 

AtiCm'^) 

^det 

AVHRR 

1800  CFig.  C18) 

0 

MSS 

67  CFig.  C8) 

1 

Td 

8  CFig.  CB) 

1 

HRV 

1  CFig.  C2) 

1 

Thermal  -  Target 

P 

exitance  56.7  UJ/m 

Table  9. 

Thermal  Sensor  Second 

Analysis 

Sensor 

Ati Cm^ ) 

^det 

Td  Band  B 

1  CFig.  67) 

1 

8.  Surviving  options:  Landsat  flSS  Cany  band) 

Landsat  Td  Band  Cany  band) 

SPOT  HRV 

9.  End  checklist  1 

Checklist  3:  Non-Cued  Detection  oF  dovinp  Targets 

1.  As  discussed  in  the  first  analysis,  all  values 

remain  at  one  Cl). 

2.  Partition  area  of  interest:  The  only  consideration  is 

cloud  cover  since  the  value  is  one  Cl)  for  both  Landsat 

and  SPOT. 
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Cloud  cover 


Cloud  cover  data  is  taken  from  Reference  IB.  For  the 
month  of  July  the  probability  of  cloud  cover  is 
approximately  O.B  throughout  the  area.  As  a  result  the  area 
of  interest  is  not  partitioned. 

3/4/5.  .Pdotl  calculation: 

Table  10  contains  the  values  for  ^dotl  each 

sensor . 


Table  10.  Second  Analysis  ^dotl 


Region 

P 

cloud 

Sensor 

^det* 

^dot  1 

Entire 

.B 

TM 

.E 

E.EBE-3 

MSS 

.2 

E,BBE-3 

HRV 

.2 

E .40E-4 

6.  The  goal  of  this  analysis  is  to  maximize  the  system  wide 
probability  of  detection  given  5100,000.00  for  the  time 
frame  of  the  analysis. 

7.  A.  The  abjective  function  is: 

Pjj>l  -  1  -  {;Cl-B.EBE-3)^^Cl-e.2BE-3)^^Cl-S.4E-4)^^}  C107) 

Where  XI  =  number  of  images  taken  by  Landsat  TH 

X2  =  number  of  images  taken  by  Landsat  nSS 

X3  =  number  of  images  taken  by  SPOT 

7.  B.  The  cost  of  a  Landsat  tlSS  image  is  51,000.00  C  43: 357). 
Using  equation  43  as  the  template  the  cost  constraint  is-. 


4600  XI  +  1000  XE  +  1000  X3  <  100000 


ClOB) 


Since  the  area  of  interest  is  not  partitioned  in  this 
analysis,  additional  constraints  designed  to  return  even 
Pdj.1  values  From  each  unique  region  are  not  necessary. 

7.  C.  This  model  was  not  solved. 


Third  Analysis:  Cued  Detection  of  Stationary  Targets 


START  —  Checklist  1 

1.  flission:  Cued  detection  of  stationary  targets 
Target:  Survey  Parties 

Time  frame:  Months  of  April,  May,  and  June 

2.  Area  of  interest:  Target  location  5  Csee  Figure  3SD 

Banks  Island  1S3*U)  73*N 

Target  location  1 
Southampton  Island  BS'Ul  B5*N 

Target  location  4 

Axel  Heiberg  Island  S5*LU  BO“N 


Figure  3E .  Fixed  Target  Locations  and  Ship  Route 


lEB 


3.  Possible  backgroundCs) :  Snow  CRef:lB) 

Ground 

H.  Background  characteristics: 

Snoui  - 

Spatial  -  Elevation  <  1000  meters  CReF:12) 

Spectral  -  Reflectance  0.3  C16:7D 

-  Emissivity  CIO-IE  pm)  .3B5  C1B:7) 

Temporal  -  Change  in  time  frame  to  ground 

Ground  - 

Spatial  -  As  above 

Spectral  -  Reflectance  0.25  C1B:7) 

-  Emissivity  ClO-12  0.94  C1B:7) 

Temporal  -  No  further  change  in  time  frame 

S.  Target  characteristics:  (Appendix  D) 

2 

Spatial  -  Vehicles  up  to  50  m 
Spectral  -  Reflectance  0.12 

-  Emissivity  .92,  1/4  surface  area  2’K 
warmer  than  ambient 

Temporal  -  No  change  in  time  frame 
5.  Sensors  available:  NQAA  AVHRR 

Landsat  5  dSS,  Td 
SPOT  1  HRV 

7.  Pjgt  calculation: 

Visible  -  Background  of  snow,  CR^  -  Rj^ )  =  .70 

-  Background  of  ground,  CR^  -  Rj^)  =  .13 


Table 

11.  Visible  Sensor 

Third  Analysis 

Sensor 

A^iCm^^Snou  ^det 

2 

^i*-"*  ^Ground 

^det 

AVHRR 

1800  CFig.  C17)  0 

••  CFig.  CS3) 

0 

nss 

67  CFig.  C7)  0 

300  CFig.  CB) 

0 

Tn  Band  4 

10  CFig.  CS)  1 

S5  CFig.  C6) 

1 

HRV 

1  CFig.  Cl)  1 

4  CFig.  C3) 

1 

Thermal  -  Target  exitance  C0  S77 

•K)  4S.0  U/m^ 

Table 

IS.  Thermal  Sensor  Third  Analysis 

Sensor 

AtiCm  ^det 

^ti*'”’  ^Ground 

^det 

TM  Band  6 

1  CFig.  G5)  1 

1  CFig.  GS) 

1 

Surviving  options:  Landsat  TM  Band  4 

Landsat  Tfl  Band  B 

SPOT  HRV 

End  checklist  1 

Checklist 

4:  Cued  Detection  of 

Stationaru  Taraets 

1.  Table  13  shows  the  number  of  opportunities  each  satellite 
has  in  the  3  month  period  of  interest. 


Table  13.  Imaging  Opportunities 


Target 

Landsat  5 

SPOT  1 

Axel  Heiberg 

41 

174 

Banks 

SS 

110 

Southampton 

B 

73 

ISB 


calculation: 


^det  Band  4  is  O.Bl 

3.  PjjQ^*  calculation: 

Only  cloud  cover  uiill  obscure  the  target  locations  during 
the  time  Frame  of  the  analysis.  Table  14  shows  the  results 
of  Pjjat*  for  each  target  and  sensor  calculation. 


Table  14.  Third  Analysis  Pj^et* 


Target 

Month 

P 

cloud 

Sensor 

^det 

^det* 

Axel  Heiberg 

April 

0.5 

TM  Band 

4 

1 .0 

0.5 

HRV 

1  .0 

0.5, 

TM  Band 

5 

1.0 

0.5 

May 

0.7 

TM  Band 

4 

1.0 

0.3 

HRV 

1.0 

0.3, 

TM  Band 

6 

1  .0 

0.3 

June 

o.a 

TM  Band 

4 

0.01 

0.16 

HRV 

1  .0 

0.0 

TM  Band 

6 

1.0 

o.e 

Banks 

April 

in 

o 

TM  Band 

4 

1  .0 

0.5 

HRV 

1.0 

0.5, 

TM  Band 

6 

1.0 

0.5 

May 

0.75 

TM  Band 

4 

1  .0 

0.05 

HRV 

1.0 

0.05 

TM  Band 

5 

1.0 

0.05 

June 

0.8 

TM  Band 

4 

0.01 

0.16 

HRV 

1.0 

0.0 

TM  Band 

6 

1.0 

0.0 

Southampton 

April 

0.6 

TM  Band 

4 

1.0 

0.4 

HRV 

1.0 

0.4, 

TM  Band 

5 

1.0 

0.4 

May 

O.B 

TM  Band 

4 

1.0 

0.0 

HRV 

1  .0 

0-2* 

TM  Band 

5 

1  .0 

0.0 

June 

0.0 

TM  Band 

4 

0.01 

0.16 

HRV 

1  .0 

0.0 

TM  Band 

6 

1.0 

0.0 

*  Discarded 

From  Further  analysis  since  the  same 

*^det* 

is  available  From 

the  same  sensor  For  the  same 

cast . 

4.  The  goal  of  this  analysis  is  to  maximize  the  probability 
of  detection  of  targets  given  SIOOOO.OO  for  the  time  frame 
of  the  analysis. 

5.  A.  The  objective  function  is: 


-  1  -  CC  .043^^C  CIOS) 

C  .a)^^C  .5)^®C  .5)^^C  .75)^^°C  .75)^^^^ 

C  .04)^^^C  .0)^^^C  .a)^^‘*C  .6)^^^C  .6)^^® 

c .a)^^^c .a)^^®c .04)^^^c .0)^^°c .b)^^^) 

Where  xl  ”  #  Tfl  Band  4  images  taken  of  target  4  in  April 
xS  =  #  HRV  images  taken  of  target  4  in  April 
x3  ”  #  Tfl  Band  4  images  taken  of  target  4  in  nay 
x4  *  #  HRV  images  taken  of  target  4  in  nay 
x5  ■  #  Tn  Band  4  images  taken  of  target  4  in  June 
xB  ”  #  HRV  images  taken  of  target  4  in  nay 
x7  =  #  Tn  Band  B  images  taken  of  target  4  in  June 
x0  ■  Tn  Band  4  images  taken  of  target  5  in  April 
x3  ■  #  HRV  images  taken  of  target  5  in  April 


xlO 

- 

# 

Tn 

Band  4 

images 

taken  of 

target 

5  in 

nay 

xll 

- 

HRV 

images 

taken 

of  target 

5  in 

nay 

xl2 

Tn 

Band  4 

images 

taken  of 

target 

5  in 

June 

xl3 

- 

HRV 

images 

taken 

of  target 

5  in 

June 

xl4 

- 

# 

Tn 

Band  6 

images 

taken  of 

target 

5  in 

June 

xl5 

- 

# 

Tn 

Band  4 

images 

taken  of 

target 

1  in 

April 

xl6 

- 

# 

HRV 

images 

taken 

of  target 

1  in 

April 

xl7 

3B 

# 

Tn 

Band  4 

images 

taken  of 

target 

1  in 

nay 

Xl0 

3K 

# 

HRV 

images 

taken 

of  target 

1  in 

nay 

xlS 

K 

# 

Tn 

Band  4 

images 

taken  of 

target 

1  in 

June 

xEO 

- 

HRV 

images 

taken 

of  target 

1  in 

June 

xSl 

* 

# 

Tn 

Band  6 

images 

taken  of 

target 

1  in 

June 

5.  B.  The  cost  constraint  is: 


4B00  xl  +  1000  xS  +  4B00  x3  +  1000  x4  +  4B00  x5  CllO) 

+  1000  xB  +  4B00  x7  +  4600  xB  +  1000  xS 

+  4600  xlO  +  1000  xll  +  4600  xlS  +  1000  xl3 
+  4600  xl4  +  4600  xl5  +  1000  xlB  +  4600  xl7 
+  1000  xlB  +  4600  xlS  +  1000  x20  +  4600  xEl  <  10000 
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5.  C.  The  image  availability  constraints  are: 


xl  +  x3  +  x5  +  x7 

< 

41 

Clll) 

xE  +  xl  +  x6 

< 

174 

CUE) 

x8  +  xlO  +  xlS  +  xl4 

< 

as 

C113) 

x9  +  xll  +  xl3 

< 

110 

C114) 

xl5  +  xl7  +  xl3  +  xSl 

< 

a 

CHS) 

xl6  +  xia  +  xSO 

< 

73 

CUB) 

5.  □.  This  model  uas  not  solved. 


Note:  If  operations  dictate  the  requirement  that  the  system 
deliver  an  equal  probability  of  detection  for  each  target, 
additional  constraints  must  be  written.  These  constraints 
would  be  similar  in  both  form  and  effect  as  those  written  to 
ensure  areas  of  common  obscuration  returned  equal 
probabilities  of  detection  in  earlier  analysis. 
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Fourth  Analusis:  Second  Contact 
START  —  Checklist  1 
1.  Mission:  Second  Contact 
Target:  Ship 

Time  frame:  Month  of  August 

S.  Area  of  interest:  Arctic  uiateruiays 

Approximately  3,000,000  km^  C13:S3) 

3.  Possible  backgroundCs) :  Uater 

4.  Background  characteristics: 

Spatial  -  Sea  level 

Spectral  -  Reflectance  O.OS  CIS: 7) 

-  Emissivity  CIO-IE  ym)  0.3S  C1B:7) 

Temporal  -  No  change  in  time  frame 

5.  Target  characteristics:  CAppendix  D) 

Spatial  -  lEl  x  15  meters 
Spectral  -  Reflectance  0.B5 

-  Emissivity  0.74,  1/10  surface  area  S*K 
warmer  than  ambient 

Temporal  -  Ship’s  top  speed  10  Knots 

6.  Sensors  available:  NQAA  AVHRR 

Landsat  5  MSS,  TM 
SPOT  1  HRV 

7.  Pfjg^  calculation: 

Visible  -  CR^  -  R^^)  =  0.76 


13S 


Table  15.  Visible  Sensor  Fourth  Analysis 


Sensor 

^i 

Cm^^^ 

^det 

AVHRR 

1800 

CFig. 

CIB) 

0 

MSS 

67 

CFig. 

CS) 

1 

TM 

0 

CFig. 

C6) 

1 

HRV 

1 

CFig. 

CB) 

1 

Thermal  -  Target 

exitance  S6 

.7  W/m^ 

Table  16. 

Thermal  Sensor  Fourth 

Analysis 

Sensor 

^ti 

Cm*^) 

^det 

TM  Band  6 

1 

CFig. 

67) 

1 

0.  Surviving  options:  Landsat  MSS  Cany  band) 

Landsat  TM  Band  Cany  band) 

SPOT  HRV 

3.  End  checklist  1 

Checklist  5:  Second  Contact 
1.  The  First  step  in  this  analysis  would  normally  be 
supplied  From  the  sensor’s  owners.  To  show  the  model 
Formulation  it  is  assumed  that  the  initial  contact  with  the 
target  is  130*UJ  75*N,  and  that  the  position  oF  Landsat  5  is 
aO*W  75*N  and  the  position  oF  SPOT  1  is  70*UJ  75*N  Five  C5) 
hours  aFter  the  initial  contact  with  the  target. 

S.  The  Landsat  5  swath  overlaps  the  possible  target  area  in 
8.3  hours  as  it  crosses  75*N  latitude  at  123.5*UJ  longitude. 
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SPOT  has  tuio  opportunities  to  image  the  possible  target  area 
within  10  hours.  The  first  occurs  in  eight  CB)  hours  at 
75*N  1S1*UI.  The  second  occurs  in  10  hours  at  75*N  I'lB'U). 


3.  The  Pjjgt  remains  at  one  Cl)  for  all  imaging 
opportunities . 

4, 5,6,7, a.  Since  Landsat  MSS  imagery  is  less  expensive  than 
Th  imagery  but  equal  in  ^det  ’  only  nSS  imagery  is  considered 
in  the  rest  of  this  analysis.  Empirical  weather  data  is 
used  to  determine  ‘  ^trkcon  are  found  in 

Appendix  F.  The  P^i^  is  calculated  based  on  the  size  of  the 
possible  target  area  and  the  Ptrkcon ‘  Table  17  contains  the 
results  for  these  calculations. 


Table  17.  Second  Contact  Calculations  Summary 


Opportunity 

^det* 

p 

trkcon 

^tim 

^trackl 

Landsat  MSS 

O.S 

0.6B 

0.34 

B.BE-2 

SPOTClst) 

O.S 

0.5 

0.045 

3.0E-3 

SPOTCEnd) 

0.2 

0.12 

0.03 

B.OE-3 

3.  The  available  swaths  could  overlap.  However,  it  is 
assumed  that  the  sensors  will  be  restricted  to  imaging 
separate  areas  to  provide  the  broadest  coverage  possible. 
3.  A.  The  goal  of  this  analysis  is  to  maximize  Pt^ack 
a  budget  of  $10,000.00. 


9.  B.  The  objective  function  is: 


P^rack  “  B.0E-S  XI  +  9.0E-3  XE  +  B.OE-3  X3  C117) 

Where  XI  -  #  images  taken  by  Landsat 

XS  “  #  images  taken  by  SPOT  on  first  pass 

X3  ”  #  images  taken  by  SPOT  on  second  pass 


The  cost  constraint  is: 


4600  XI  +  1000  XS  +  1000  X3  <  10000  Cl 18) 


The  image  availability  constraints  are: 


0 

s 

XI 

<  E 

C119) 

0 

< 

xa 

<  11 

ClEO) 

0 

< 

X3 

<  4 

ClSl) 

9.  0.  This  linear  program  uias  not  solved. 


Fifth  Analysis;  Follouj-on  Contact 
START  —  Checklist  1 
1.  Mission :  Second  Contact 
Target:  Ship 

Time  frame:  June,  Duration  of  the  voyage 
S.  Area  of  interest:  Ship’s  Route  which  begins  at  BO’N  60*W 
and  proceeds  Northwest  until  the  ship  has  cleared  the  area. 
The  ship’s  route  is  shown  in  Figure  3S. 

3.  Passible  backgroundCs) :  Water 

4.  Background  characteristics: 

Spatial  -  Sea  level 

Spectral  -  Reflectance  0.03  CIS: 7) 

-  Emissivity  CIO-IS  um)  0.33  C16:7) 

Temporal  -  No  change  in  time  Frame 

5.  Target  characteristics;  (Appendix  D) 

Spatial  -  121  x  15  meters 
Spectral  -  Reflectance  0.85 

-  Emissivity  0.74,  1/10  surface  area  B*K 
warmer  than  ambient 

Temporal  -  Ship’s  top  speed  10  Knots 
B.  Sensors  available:  NOAA  AVHRR 

Landsat  5  MSS,  TM 
SPOT  1  HRV 

7.  Pjjg^  calculation; 

Visible  -  (R^.  -  Rj^)  =  0.7B 
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Table  IB.  Visible  Sensor  Fifth  Analysis 


Sensor  D 


AVHRR  1800  CFig.  CIB)  0 
nSS  67  CFig.  CB)  1 
Tn  B  CFig.  CB)  1 
HRV  1  CFig.  CB)  1 


1 


Thermal  -  Target  exitance  56.7  U/m 


Table  13. 

Thermal  Sensor  Fifth 

Analysis 

Sensor 

A^^Cm'^) 

^det 

TM  Band  6 

1  (Fig.  G7) 

1 

a.  Surviving  options:  Landsat  hSS  Cany  band) 

Landsat  Th  Band  Cany  band) 
SPOT  HRV 

S.  End  checklist  1 


3U)-Qn  Contact 

1,B,3  A.  The  total  time  for  the  voyage  is  approximately  B50 


hours . 


B.  Given  an  initial  Landsat  5  location  of  110*UJ  BO’N, 
and  a  SPOT  location  of  30*W  60*N  at  the  start  of  the  ship’s 
voyage  (arbitrary  for  demonstration  purposes),  there  uill  be 
five  C5)  imaging  opportunities  during  the  voyage.  These  are 
shown  in  Figure  33. 


! 
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Longitude 


Figure  33.  Scatter  Plot  of  Imaging  Oppartunities 


4.  The  imaging  opportunities  and  values  are  summarized 

in  Table  20. 


Table  20.  Imaging  Opportunities 


TimeChours) 

Sensor 

P 

cloud 

^det* 

47 

SPOT 

0.8 

0.2 

71 

SPOT 

0.8 

0.2 

94 

Landsat 

0.8 

0.2 

144 

SPOT 

0.8 

0.2 

218 

SPOT 

0.8 

0.2 

I 

) 

i 

f 

I 
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5.  Pjj>i  is  to  be  maximized  in  this  example.  The  available 
budget  is  52,000.00.  Therefore  the  objective  Function  is: 

Pd>i  -  1-CC  C122) 

Where  XI  “  number  of  images  taken  by  SPOT  on  its 
first  opportunity 

X2  “  number  of  images  taken  by  SPOT  on  its 
second  opportunity 

X3  -  number  of  images  taken  on  the  Landsat 
pass 

XH  -  number  of  images  taken  by  SPOT  on  its 
third  opportunity 

X5  =  number  of  images  taken  by  SPOT  on  its 
fourth  opportunity 


The  cost  constraint  is; 


1000  XI  -*■  1000  X2  +  1000  X3  C123) 

+  1000  X4  +  1000  X5  <  2000 


The  image  availability  constraint  is  that  each  sensor 
can  provide  only  one  image  of  the  expected  ship  location  on 
each  pass . 

B.  This  model  mas  not  solved. 


VII.  Conclusions  and  Recommendations  For  Further  Research 

General 

Satellite  remote  sensing  of  the  Earth  for  any  purpose 
is  a  complex  process  involving  the  scene,  the  atmosphere, 
and  the  sensor.  Equally  complex  is  the  task  of  evaluating 
the  images  that  are  available  from  remote  sensing 
satellites.  The  sensor  evaluation  procedures  and 
surveillance  mission  evaluation  algorithms  developed  in  this 
research  consider  only  a  small  part  of  the  entire  process. 

As  a  result  of  the  narrow  focus  taken  in  this  research 
some  questions  can  be  answered  while  others  remain. 

Therefore  this  concluding  section  of  this  report  separates 
conclusions  into  two  categories.  The  first  category  is  hard 
conclusions,  those  that  are  defendable.  These  conclusions 
are  presented  in  the  section  which  addresses  the  research 
questions  posed  in  Chapter  1 .  The  second  category  is  soft 
conclusions,  those  that  are  observations  on  the  interaction 
of  satellites  and  surveillance.  These  observations  are 
presented  in  their  own  section. 

The  final  section  of  this  chapter  presents  some  areas 
which  cQuid  be  investigated  by  follow-on  research. 
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Discussion  oF  Research  Questions 

Uihat  satellite  imaaino  sustems  mill  be  available? 

The  United  States,  Europe  CESA),  India,  Japan,  Canada, 
the  Soviet  Union,  and  France  plan  to  sell  satellite  acquired 
imagery  in  the  1390s.  Additional  sources  of  satellite 
acquired  imagery  may  include  Chinasat,  Hediasat,  Bressex, 
and  a  Soviet  digital  imaging  system. 

Uhat  mill  the  sustem’s  capabilities  be  for  these  suppliers? 

The  minimum  target  size  that  can  be  detected  by  a 
visible  wavelength  sensor  is  determined  by  the  sensor’s 
ground  resolution  ceii  CA^,),  the  sensor’s  noise  equivalent 
reflectance  difference  CNERO),  the  atmospheric  transmission 
CTg),  the  reflectance  of  the  target  CR^),  and  the 
reflectance  of  the  background  CR^) . 

A_  CNERD/t_) 

A  -  — E - ^  C124) 

I  Rt  -  «b  I 

The  minimum  target  that  can  be  detected  by  a  thermal 
wavelength  sensor  is  determined  by  the  sensor’s  noise 
equivalent  power  difference  CNEPD),  the  atmospheric 
transmission  ,  the  background  exitanc-  the 

emissivity  of  the  background  the  target  exitance 

ChtCTt^),  and  the  target  emissivity  . 


NEPD/t 

A  -  - 2 -  C125) 

I  I 


Analysis  of  visible  and  thermal  sensors  must  consider 
the  detector  array  geometry  and  its  affect  on  the 
probabiiity  of  target  detection.  For  both  types  of  sensor 
the  probability  of  target  detection  is  the  ratio  of  the  area 
in  which  the  target  could  fall  and  be  detected  CA^^g^)  to  the 
total  area  in  which  the  target  could  fall  f ^tot^ ' 


^det  ^det'^^tot 


C12B) 


How  can  commercial  satellite  images  be  evaluated? 

Commercial  satellite  images  can  be  evaluated  by 
comparing  the  information  they  provide  to  the  information 
required,  and  by  considering  the  cost  of  the  information 
they  provide.  The  specific  algorithms  used  For  sensor 
evaluation  against  a  mission  requirement  are  mission 
specific  and  can  be  optimized  using  linear  and  non-linear 
optimization  techniques. 

For  non-cued  detection  of  stationary  targets  the  system 
wide  probability  of  detecting  a  target  is  given  below  in 
equation  1E7 . 


j 
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dsys 


i-1 


C127) 


Uhera  ^dsus  "  system  wide  probability  of  detection 
CN^^Ji  "  #  of  images  taken  by  sensor  i 
CP^otl^i  “  probability  of  target  detection  in 
area  of  interest  in  a  single  image 
from  sensor  i 


For  non-cued  detection  of  moving  targets  the  system 
wide  probability  of  at  least  one  target  detection  is: 


Pd>i  =  1  -  c  n  ti-Pdotii^'''^ 

1*1 

Where  Pd>l  ”  system  wide  probability  of  at  least 

one  detection 

PdQtii  *  probability  of  target  detection  in 
area  of  interest  in  a  single  image 
from  sensor  i 

xi  *  #  of  images  taken  by  sensor  i 


For  cued  detection  of  stationary  targets  the  system 
wide  probability  of  at  least  one  target  detection  is: 


d>l 


1  -  C 


n 

n 

i*l 


Cl-P 


det*i 


^xi:^ 


C1B3) 


Where  ^det*!  *  probability  of  target  detection  in  a 

single  image  from  sensor  i 


For  second  contact  on  moving  targets  the  system  wide 


probability  of  tracking  is  given  in  equation  130. 


C130) 


I  ^track  "  ^trackli 

Where  *^track  “  system  wide  probability  of  tracking 
^trackli  ”  probability  of  tracking  the  target 
in  a  single  image  from  sensor  i 

I 

This  relationship  holds  for  overlapping  and  non¬ 
overlapping  images  by  considering  overlapping  images  to  be 
I  available  from  a  separate  sensor  having  a  P^rackl 

from  the  P^rackl  sensors  involved. 

For  folloui-on  contact  the  system  wide  probability  of  at 
i  least  one  detection  is; 


P 


d>l 


c  n 

i  =  l 


m 


n  ci-p 

j-i 


det*i  j 


)} 


C131) 


Where  ^det*i 1  ^  probability  of  detection  in  image  J 

of  sensor  i 


For  all  missions  the  constraint  set  is: 


£  ciae) 

1-1 

Where  =  cost  of  one  image  from  sensor  i 

Bj  -  budget  available  in  the  time  frame 

0  <  <  K.  C133J 

im  1  1 

Where  =  maximum  number  of  images  available  from 

sensor  i  in  the  time  frame 


For  non-cued  detection,  additional  constraints  are 
required  to  ensure  that  more  easily  imaged  areas  within  the 


area  of  interest  are  not  selected  over  areas  that  are  more 
difficult  to  image. 

Can  weather  satellites  detect. and  track  targets  of  interest? 

Current  weather  satellites  cannot  detect  or  track 
objects  of  interest  in  the  Canadian  North.  This  is  because 
geosynchronous  satellites  do  not  provide  images  of  Northern 
Canada,  and  the  minimum  target  size  for  detection  by  the 
NOAA  series  of  satellites  is  approximately  EOOO  square 
meters . 

Uiill  commercial  sustems  require  cuing? 

The  decision  on  whether  to  use  commercial  imaging 
sensors  as  cues  is  economic.  The  system  wide  probability  of 
detection  is  limited  to  a  maximum  which  corresponds  to  the 
probability  that  there  are  clear  skies  in  the  target  area. 

To  achieve  the  maximum  probability  of  detection  the  entire 
area  of  interest  would  have  to  be  imaged.  Assuming  that  the 
price  per  image  remains  constant  irrespective  of  the  number 
of  images  acquired,  the  cost  of  imaging  the  entire  area  is 
$212,000.00  if  Landsat  flSS  imagery  is  used.  The  cost  of 
imaging  the  entire  area  is  $376,000.00  if  Landsat  TM  imagery 
is  used.  The  cost  of  imaging  the  entire  area  is 
$2,000,000.00  if  SPOT  imagery  is  used.  If  this  imaging  were 
undertaken  in  August  the  system  wide  probability  of 
detecting  a  target  located  somewhere  in  the  target  area 
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uould  be  approximately  O.B  due  to  the  likelihood  of  clouds. 
Imaging  of  the  entire  area  mould  take  Landsat  16  days.  It 
mould  take  SPOT  5  days . 


If  cued,  Landsat  5  and  SPOT  1  sensors  mill  detect 
targets  of  interest  in  Northern  Canada  mith  a  probability  in 
a  single  image  equal  to  the  probability  of  clear  skies  in 
the  target  area.  Since  the  probability  of  clear  skies 
varies  both  spatially  and  temporally,  the  specific  location 
mill  influence  the  number  of  images  to  be  acquired  to  return 
an  acceptable  probability  of  at  least  one  detection. 

mil  a  combination  of  sensors  outperform  a  single  sensor'!* 

The  question  of  sensor  or  system  dominance  can  be 
decided  in  general  terms  for  each  mission.  This  is 
discussed  in  more  detail  in  the  observations  section  of  this 
chapter . 

Could  an  alternative  system  outperform  commercial  sensors? 

An  alternative  system  mould  have  to  outperform  the 
combination  of  sensors  available  from  commercial  sources. 
Therefore  it  mould  have  to  return  larger  system  mide 
probabilities  of  detection  and  tracking  for  the  same  cost, 
or  mould  have  to  return  equal  system  mide  probabilities  of 
detection  and  tracking  for  less  cost  than  the  commercial 
systems . 
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Can  commercial  sustems  detect  Canadian  militaru  activitu? 

Canadian  military  vehicles  and  ships,  if  they  have 
similar  characteristics  to  those  used  in  this  analysis,  can 
be  detected  and  tracked  by  current  ccmmercial  satellite 
sensors . 

Observations 

The  sixth  research  question  asks  whether  a  mix  of 
sensors  performs  better  than  does  a  single  sensor.  The 
short  answer  is  that  it  depends  on  the  mission.  A  brief 
discussion  of  each  mission  explains  this  point. 

In  non-cued  detection  the  goal  is  wide  area 
surveillance.  A  comparison  of  sensors  is  possible  through 
their  respective  costs.  The  P^latl  given  by 

equation  134. 


^dotl 


^one  ^det* 


'int 


C134) 


Where  ^one  ”  ^^'ea  of  a  single  image  from  the  sensor 
A^^^  -  area  of  interest 


Given  competing  sensors  have  equal  ^det* ’  and  since  the 
area  of  interest  is  the  same  for  each  sensor,  a  sensor 
producing  a  larger  image  will  be  favored  over  one  which 
produces  a  smaller  image.  The  solution  of  the  linear 
program  developed  in  the  first  analysis  shows  the  domination 
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□f  SPOT  by  Landsat  for  the  uiide  area  surveillance  mission. 
These  results  are  shown  in  Table  El . 


Table  El.  Summary  of  First  Analysis  Results 


Sensor/Area 

#  Images 

Landsat/B 

E3.a 

Landsat/C 

38.  S 

Landsat/O 

88. E 

Landsat/E 

17.0 

Landsat/F 

ES.O 

Landsat/G 

18.5 

SPOT  /B 

0.0 

SPOT  /C 

0.0 

SPOT  /F 

0.0 

SPOT  /G 

0.0 

The  explanation  for  this  domination  is  the  cost  of  the 
imagery.  Landsat  TM  imagery  covers  territory  at  a  cost  of 
50.13  each  square  kilometer.  SPOT  HRV  covers  territory  at  a 
cost  of  SO.ES  each  square  kilometer. 

In  cued  detection  of  stationary  targets  the  size  of  the 
image  is  not  important  in  the  selection  of  the  sensor.  This 
is  true  since  sensor  performance  is  compared  using 
individual  does  not  consider  the  relative  sizes 

of  the  swath  and  area  of  interest.  In  this  analysis  the 
^det*  established  for  Landsat  and  SPOT  as  the  probability 
of  there  being  clear  weather  in  the  target  area.  The  two 
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factors  which  are  important  in  sensor  selection  are  the  cost 
of  the  sensor  and  the  responsiveness  of  the  sensor.  For 
high  value  targets  it  may  be  necessary  to  confirm  or  deny 
its  presence  as  quickly  as  possible.  In  this  situation  the 
selection  of  a  sensor  will  be  driven  by  the  ability  of  the 
sensor  to  image  an  area  before  a  competitor. 

Selection  of  a  sensor  for  second  contact  with  a  target 
again  depends  on  the  situation.  The  preferred  option  will 
provide  the  highest  P^rkcon  those  available.  An 
interesting  outcome  of  the  analysis  is  that  pointable 
sensors  provide  slightly  higher  Ptrkcon  if  their 

ground  trace  does  not  overfly  the  initial  target  location 
directly.  This  is  shown  in  Figure  34.  The  reason  this  is 
true  is  that  the  sensor’s  swath,  and  thus  the  area  it  can 
cover,  is  slightly  larger  off  nadir. 

Finally,  attempting  to  acquire  Follow-on  contacts  with 
targets  Favors  pointable  sensors.  The  selection  of  one 
sensor  over  another  is  again  situation  dependent  since 
several  imaging  opportunities  may  be  present  For  each 
sensor.  The  final  selection  From  the  available 
opportunities  will  depend  on  the  value  of  the  target  and  the 
requirement  to  monitor  the  target  as  frequently  as  its  value 
demands . 


! 
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Figure  34.  SPOT  1  Conditional  Tracking  Probability 
60*N  Latitude,  5  Knots 

CTop  to  bottom  are  0,4,0,10,12*  Longitude  Displacement) 


Recommendations 

Detection  throughout  this  research  refers  to 
distinguishing  a  target  from  its  background  in  the  spatial 
domain.  It  is  also  passible  to  approach  detection  in  the 
temporal  domain.  This  requires  prior  information  on  the 
area  of  interest  and  the  development  of  algorithms  to 
separate  out  evidence  of  change  over  time.  Folloui-on 
research  could  address  the  problem  of  detection  in  the 
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temporal  domain  and  investigate  the  data  storage  and 
processing  u/hich  uiouid  be  required  to  accomplish  such 
detection . 

This  research  assumes  that  the  timeliness  and 
availability  of  images  is  ideal.  Given  that  these  factors 
are  not  necessarily  ideal,  houi  could  they  be  evaluated  in 
the  overall  determination  of  the  worth  of  satellite  imagery. 
Follow-on  research  could  develop  the  algorithms  to  use  in 
the  evaluation  of  timeliness  and  availability  and  could 
develop  an  overall  value  algorithm  for  satellite  acquired 
imagery  . 

The  third  mission  of  surveillance  is  identification. 
What  is  involved  in  this  mission?  Should  the  military  train 
photo-interpreters  or  develop  pattern  recognition  algorithms 
to  effect  identification?  Follow-on  research  could  identify 
the  infrastructure  and  costs  involved  in  these  options  and 
evaluate  the  options  against  the  mission  requirements. 

The  fourth  major  area  of  research  which  could  be 
pursued  is  the  addition  of  active  sensor  evaluation  criteria 
and  investigation  of  the  impact  of  active  sensors  in  the 
selection  of  imaging  sensors. 
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Appendix  A:  Sensor  Description 

This  appendix  contains  a  general  description  of  imaging 
sensors  that  are  either  current  or  are  planned  to  be 
launched  in  the  1390s.  Excluded  from  this  listing  are 
sensors  that  are  incapable  of  imaging  Northern  Canada, 
examples  of  uihich  are  the  GOES  geosynchronous  weather 
observation  satellites  and  the  Tropical  Earth  Resources 
Satellite  CTERS^  which  is  planned  to  fly  at  0*  inclination. 

There  are  four  notable  omissions  from  this  appendix. 

The  first  is  Bressex ,  a  Brazilian  shuttle-borne  system.  It 
is  omitted  since  little  information  is  available  on  this 
system.  The  other  omitted  systems  are  Chinasat,  Mediasat, 
and  a  possible  Soviet  digital  imaging  system.  Although 
these  systems  may  be  developed  and  become  available 
commercially  in  future  their  current  status  is  speculative 
at  best . 
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Earth  Observation  System  CEOS)  Polar  Platforms 


Country 

Launch 

Orbit 

Lifetime 

Main 

of  Origin 

Date 

Parameters 

Sensor 

USA  1S34  1  “  SB*  10  years  Active/ 

H  “  BE4  km  Passive 

EQX  “  0900  Descending 
H  -  BEH/54E/705  km 
EQX  ■  1330  Ascending 

Background:  The  EOS  is  sponsored  by  NASA  and  is  conceived 

as  the  next  major  operational  Earth  remote  sensing  system. 
The  two  polar  platforms  mill  carry  a  variety  of  sensors  and 
will  in  part  replace  the  NOAA  tueather  observation  satellite 


Imaging  Sensors ; 


Name 

Band 

BandpassC  ym) 

GRCCm) 

NERDC 

nODIS-N 

1 

.46- .48 

500 

.E- 

.5 

E 

.54-. 56 

500 

.E- 

.5 

3 

.66-. 68 

500 

.8- 

.5 

4 

.70- .7E 

500 

.E- 

.5 

5 

.87-. 89 

500 

.E- 

.5 

6 

.33-. 95 

500 

.E- 

.5 

7 

.43-. 44 

1000 

6 

.485-. 495 

1000 

9 

.515-. 5S5 

1000 

10 

.56- .57 

1000 

11 

.615-. 6E5 

1000 

IE 

.66-. 67 

1000 

13 

.68-. 69 

1000 

14 

.76- .77 

1000 

15 

.86- .87 

1000 

El 

1.07-1 .09 

500 

.E- 

.5 

EE 

1  .lE-1 .14 

500 

.E- 

.5 

E3 

1.54-1.56 

500 

.E- 

.5 

E4 

1.63-1 .65 

500 

.E- 

.5 

E5 

S.Ol-E.ll 

500 

.E- 

.5 

E6 

E.OB-S.ia 

500 

.s- 

.5 

nODIS-T 

? 

o 

1 

»-> 

o 

1000 

SujathCkm) 


1500 


500 


1  53 


AHRIR 


1-11 


VIS-NIR 


Name 


Band  BandpassC  jjim)  GRCCm)  NERDC/i)  SuathCkm) 


Tins 

7 

.40-12.5 

SAR 

1 

isnHz 

ITIR 

1 

.55-. 32 

2-5 

1.5-2.36 

7-11 

3.53-11.7 

HIRIS 

7 

.45-2.45 

30 

05  33-100 

15 

15 

50 

30  23 


nODIS-N 

nODIS-T 

AMRIR 

ITIR 

HIRIS 

Tins 

SAR 


noderate  Resolution  Imaging  Spectrometer  —  Nadir 
(This  nODIS  is  Fixed  at  nadir) 

T  is  for  tilt  along  track  ±  50* 

Advanced  Hedium  Resolution  Imaging  Radiometer 
CFolioui-on  to  NOAA  AVHRR) 

Intermediate  Thermal  Infrared,  proposed  by  Japan 
High  Resolution  Imaging  Spectrometer 
Thermal  Imaging  Spectrometer 
Synthetic  Aperture  Radar 


Other :  The  HIRIS  will  be  pointable  across  track  ±  20*  and 
along  track  ±  60*,  and  is  intended  to  be  a  targeted  sensor 
rather  than  in  continuous  operation. 


Sources :  C  23 ; 77 ) C  34 : 02 ) C 1 7 ; 105 ) C  BE : 352 ) C  3 1 : 37 1 ) ( 57 : 388 ) 
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European  Space  Agency  Resources  Satellite  CERS-1) 


Country 

Launch 

Orbit 

Lifetime 

of  Origin 

□ate 

Parameters 

ESA  1389  i  -  38.5*  3  years 

H  -  777  km 
C  *  3  days 

EQX  *  1030  Descending 


Main 

Sensor 


Active 


Background :  ERS-1  will  be  the  first  Earth  remote  sensing 

satellite  for  the  ESA.  It  is  intended  to  be  experimental 
and  will  be  the  forerunner  for  follow-on  satellites.  If 
ERS-1  is  a  success,  ERS-2  could  be  launched  in  1393. 


Imaging 

Sensors : 

Name 

Band 

Frequency 

GRCCm) 

SwathCkm) 

AMI 

C 

5.3  GHz 

30 

75 

AMI  -  Advanced  flicrowave  Instrument 


Other :  ESA  plan  to  make  fast  delivery  products  available 

within  three  hours.  A  fast  delivery  reception  facility  is 
planned  for  Gatineau. 


Sources :  C  53 : 70 ) C  25 : 537  H  32 : 375 ) 
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Indian  Remote  Sensing  Satellite  CIRS) 


Country 

Launch 

Orbit 

Lifetime 

of  Origin 

Date 

Parameters 

India 

? 

i  -  99* 

3  Years 

H  -  904  km 
C  -  82  days 


Background:  The  Indian  remote  sensing  satellites  are 

follouj-ons  to  the  Bhaskara  1  and  8  satellites  flown  in  the 
late  1970s  and  early  19B0s.  The  IRS  program  is  planned  as 
prB-operational  space-based  remote  sensing  system. 


Imaging  Sensors: 


Name 

Band 

BandpassCym) 

GRCCm) 

NERDC^) 

Swathe  km ) 

LISS-1 

1 

.45- .58 

78 

140 

8 

.58- .60 

3 

.63-. 69 

4 

.76-. 90 

LISS-8* 

1 

.45-. 58 

36 

.5 

m 

148 

8 

.58-. 59 

.5 

3 

.68-. 60 

.35 

4 

.77-. 06 

.3 

LISS  -  Linear  Imaging  Self-Scanning 

•  There  are  two  LISS-8  sensors,  each  swath  is  74  km. 


Other :  Three  levels  of  data  products  are  planned.  Level  1 

will  be  corrected  for  radiometric  errors  and  Earth  rotation 
and  will  be  delivered  in  three  days.  Level  8  will  be 
corrected  for  geometric  errors  and  will  be  delivered  in 
seven  days.  Level  3  are  referred  to  as  precision  products 
and  will  be  delivered  within  three  weeks  of  imaging. 


Sources:  C48 :  07K5B :  615)  C59 630)  C60 :  153) 


JANUS  Earth  Observation  Satellite  CJEOS) 


Country 

Launch 

Orbit 

Lifetime 

Main 

of  Origin 

□ate 

Parameters 

Sensor 

ESA 

i  -  SB. 5* 

4  years 

Passive 

H  -  7B0  km 
C  “  2S  days 
EQX  -  1000 


Background:  This  is  a  proposed  Earth  remote  sensing 

satellite  that  could  be  inserted  into  orbit  in  piggy-back 
Fashion  during  an  ESA  launch.  The  goal  is  to  produce  a  low 
cost  satellite. 


Name 

Band 

BandpassC  >im ) 

GRCCm) 

NERDC*/.)  5wathCkm) 

1 

.45- 

.52 

30 

.5  120* 

2 

.52- 

.53 

.5 

3 

.62- 

.68 

.35 

4 

.77- 

.86 

.3 

Source:  C55:B15) 


Japanese  Earth  Resources  Satellite  CJER5-1) 


Country 

Launch 

Orbit 

Lifetime 

of  Origin  . 

Date 

Parameters 

Japan  1S90  I  -  93.7* 

H  -  568  km 


Background :  In  1980  the  National  Space  Development  Agency 

of  Japan  began  a  research  and  development  program  for  remote 
sensing.  Both  the  JERS-1  and  HOS-l  are  products  of  this 
program . 


Imaging  Sensors: 


Name  Band  BandpassCyni)  GRCCm)  NERD  SujathCkm) 


VNR  1 

a 

3 

4 

SAR  Frequency 


.‘15-.5E  as 

.sa-.Bo 

.63-. 69 
.76-. 90 

i.a75  GHz  asxas 
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75 


VNR  -  Visible  and  Near-Infrared  Radiometer 
SAR  -  Synthetic  Aperture  Radar 


Sources:  C36 :  611 K53 :  71 ) 


Landsat  5 


Country 

Launch 

Orbit 

Lifetime 

of  Origin 

Date 

Parameters 

USA  March  i  -  98*  1981* 

1904  H  -  705  km 

C  ~  16  days 
EQX  ”  0945  Oescending 

Background :  Landsat  5  is  the  Fifth  land  satellite  to  be 

launched.  This  program  was  originally  funded  by  the  US 
government  but  became  a  commercial  venture  in  1965 .  Delays 
are  expected  for  the  launch  of  Landsat  B  so  Landsat  5  is 
being  operated  below  capacity  in  order  to  extend  the  life  of 
the  satellite. 

Imaging  Sensors : 


Name 

Band 

BandpassC  >im  ) 

GRCCm) 

NERDCV.) 

Swathe  km ) 

TM 

1 

.45-. 52 

30 

.65 

105 

2 

.52-. 60 

.57 

3 

.63-. 69 

.57 

4 

.76-. 90 

.33 

5 

1.55-1.75 

1.68 

7 

2 . 08-2 . 35 

2.00 

6 

10.4-12.5 

120 

NETD  .5 

•K  e  3oo*k: 

MSS 

1 

.50- .60 

00 

2 

.60- .70 

3 

.70- .00 

4 

.80-1.1 

TM  -  Thematic  Mapper 

MSS  -  Multispectral  Scanner 

Other :  Data  from  Landsat  satellites  can  be  relayed  through 

the  tracking  and  data  relay  satellites  CTDR5) .  Initial  MSS 
processing  is  accomplished  at  the  Goddard  Space  Flight 
Center  CGSFC)  .  Follouj-on  processing  and  storage  is  done  at 
the  Earth  resources  observation  system  data  center  in  Sioux 
Falls,  South  Dakota.  TM  data  can  be  Fully  processed  at 
GSFC. 

Sources :  C  56 : 356 ) ( 33 ; 63-73  j  C  64 : 438  5 


Landsat  B/7 


Country  Launch  Orbit  Lifetime 

of  Origin  Date  Parameters 

USA  1909/1991  i  -  90* 

H  -  705  km 
C  ~  16  days 
EQX  -  09*15  Descending 

Background :  The  future  of  Landsat  6  and  7  is  unclear. 

Landsat  5  was  originally  scheduled  for  launch  in  1900.  Now 
estimates  for  a  launch  date  range  to  1993.  Landsat  7  plans 
continue  however  funding  for  this  platform  is  questionable. 


Imaging  Sensors: 


Name  Band  BandpassC>im)  GRCCm)  NERDC?s)  SwathCkm) 


ETh  C7  spectral  as  found  in  Landgat  5) 


Pan 

.50-. 90 

15/5 

3.53-3.93 

lEO 

a* 

B.EO-0.75 

BO 

9* 

0.75-9.30 

BO 

lO.E-11 .0 

BO 

11 

11.0-11.3 

BO 

m 

from  VNIR 

10 

4 

from  SUIR 

EO 

ETM  -  Enhanced  Thematic  Mapper,  this  sensor  will  be  able  to 
emulate  the  MSS  found  on  Landsat  5 
ALS  -  Advanced  Landsat  Sensor,  possible  for  Landsat  7 
VNIR  -  Visible  and  near  infrared  wavelengths 
SUIIR  -  Short  wave  infrared  wavelengths 
*  -  Under  consideration 

Other :  EOSAT  has  conducted  a  market  survey,  which  included 

potential  customers  from  the  media,  to  investigate  the 
possible  market  for  5  meter  resolution  imagery.  Also  EOSAT 
has  declared  that  TORS  will  not  be  used  to  relay  imagery  of 
foreign  countries  in  near  real  time.  Instead  such  imagery 
will  be  recorded  for  subsequent  transmission. 

Sources:  C7E:74K11:  149)  C27:  73KE0:  005)  C33:  109) 
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Maritime  Observation  Satellite  CMOS-1) 


Country 

Launch 

Orbit 

Lifetime 

of  Origin 

Date 

Parameters 

Japan  130S  i  “  S3*  S  years 

H  *  308  km 

EQX  “  1100  Descending 


Background:  MOS-1  is  the  maritime  portion  of  a  Japanese 
research  and  development  program  to  establish  fundamental 
technologies  for  remote  sensing  satellites. 


Imaging  Sensors: 


Name 

Band 

Bandpass  C;im) 

GRCCm)  NERDC?i) 

SuathC 

MESSR 

1 

.51-. 53 

50 

100 

E 

.61-. 63 

3 

.7E-.ao 

4 

.80-1.1 

VTIR 

1 

.50-. 70 

300 

1500 

E 

6. 0-7.0 

E700 

3 

10.5-11 .5 

E700 

4 

11 .5-1E.5 

E700 

MESSR  -  Multispectral  electronic  self-scanning  radiometer, 
there  are  tu/o  on  MOS-1 

VTIR  -  Visible  and  thermal  infrared  radiometer 


Sources:  C33: 187) C53: 70) C44 ; 535) C43: 18) 
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NOAA  I/J/K/L/n 


Country 

Launch 

Orbit 

Lifetime 

of  Origin 

□ate 

Parameters 

USA  I /Oct  ISBB  i  -  3B.9* 

J/nar  IBBB  H  -  BBB  km 
K/Oct  1930 


Background;  NOAA  stands  for  the  National  Oceanic  and 
Atmospheric  Administration.  The  satellites  that  bear  this 
name  are  weather  observation  satellites.  NOAA  satellites 
are  folloui-ona  in  a  continuing  program  which  began  with  the 
launch  of  TIROS-1  on  April  1,  1960.  Typically  NOAA 
satellites  operate  two  at  a  time  and  provide  both  morning 
and  afternoon  imaging  of  weather  systems. 

Imaging  Sensors: 

Name  Band  BandpassC^im)  GRCCm)  NERDC^i)  SwathCkm) 


AVHRR  1 


3b 

4 

5 


.58-. 68  1000 

.73-1 .0 

I. 57-1.78 
3.53-3.93 
10.2-11.3 

II . 5-12.5 


3000 


§VHRR  -  Advanced  Very  High  Resolution  Radiometer 
-  Planned  for  NOAA  K/L/fl 


Other ;  NOAA  satellites  transmit  images  in  real  time 
directly  to  users.  The  resolution  at  which  those  images  are 
received  depends  on  the  ground  station  built.  Automatic 
picture  transmission  CAPT)  images  are  received  at  a  ground 
resolution  of  four  kilometers.  High  resolution  picture 
transmission  CHRPT)  images  are  received  at  a  ground 
resolution  of  one  kilometer. 


Sources :  C  65 : 44 1 ) C  63 : 52-67 ) C  37 : 22 1 K  58 : 3 ) 


i 
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RADARSAT 


Country 

Launch 

Orbit 

Lifetime 

of  Origin 

Date 

Parameters 

Canada  i  "  90.5*  5  years 

H  -  798  km 
C  “  16  days 
EQX  -  1030  Descending 

Background :  RADARSAT  is  Canada’s  first  Earth  observation 

satellite.  RADARSAT  is  intended  to  be  both  research  and 
operations  oriented.  The  main  sensor  uill  be  an  active 
microwave  imaging  radar.  Until  recently  the  United  Kingdom 
planned  to  develop  a  radiometer  to  fly  with  the  radar, 
however  the  UK  has  pulled  out  of  RADARSAT.  As  a  result,  the 
inclusion  of  a  radiometer  has  not  been  finalized. 


Imaging  Sensors : 

Name  Band  Frequency  GRCCm)  Swath C km) 

SAR  C  5.3  Ghz  E5CR)x8a(A2)  —  7  main  500* 

35CR)x2aCA2)  —  2  wide-swath  BOO 
10CR)x9CA2)  —  5  high  45 

resolution 

CR)  “  Range 
CA2)=  Azimuth 

*  -  The  seven  beams  each  have  a  swath  of  100  km  but  they 
overlap . 

Other:  Ground  stations  for  RADARSAT  include  Fort  Churchill, 

Shoe  Cove,  Gatineau,  Prince  Albert,  Fairbanks,  and  Stranraer 
in  Scotland.  Data  relay  will  be  possible  through  an  ANIK 
satellite . 


Sources :  C  53 : 70 ) C  45 : 33 ) C  46 : 56 1 ) C  47 : 409 ) Cl : 6 ) 
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Country 

Launch 

Orbit 

Lifetime 

Main 

of  Origin 

Date 

Parameters 

Sensor 

USSR  ?  ?  ?  ? 


Background:  Sojuzkarta  is  a  commercial  firm  established  by 

the  government  of  the  Soviet  Union  to  sell  satellite 
acquired  imagery  to  the  uorld.  Few  details  are  available  on 
platforms,  orbits,  and  sensors.  However,  Sojuzkarta  provide 
a  general  description  of  their  sensors  with  their  order  form 
and  price  list.  All  sensors  listed  below  are  cameras  and 
rely  on  physical  retrieval  of  film  canisters  to  transport 
the  film  to  the  ground. 


Imaging  Sensors: 


Name 

Altitude 

Ckm) 

Spectral 
Range  C>im) 

Resolution 

Cm) 

Swath 

Ckm  X  km) 

KATE-aOO 

250 

.50-.  60 
.60-. 70 
. 70- . 30 

15-30 

180  X  180 

MKF-B 

250 

.46-. 52 
.52-. 56 
.58-. 62 
.64- .60 
.70-. 74 
.7S-.90 

20 

140  X  200 

KATE-140 

250 

.50-. 70 

60 

270  X  270 

KFA-1000 

250 

.57-. 67 
.67-. 80 

5-10 

60  X  60 

Sources;  C52 : 1346) CBS ; 140) 
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SPOT  l/a/3 


Country 

Launch 

Orbit 

Lifetime 

of  Origin 

□ate 

Parameters 

France 

Feb  1366 

i  -  96. 5* 

2-3  years 

H  -  aas  km 
C  “  5  days 

EQX  ”  1030  Descending 


Background :  Systeme  Prcbatoire  d ’Observation  de  la  Terre, 

or  SPOT  began  as  an  ESA  project  but  changed  into  a  French 
lead  program  involving  only  Sweden,  Belgium,  and  private 
interests.  It  is  currently  overseen  by  the  Centre  National 
d’ Etudes  Spatial  CCNESD,  a  French  government  agency. 


Imaging  Sensors ; 


Name 

Band 

Bandpass  C;im) 

GRCCm) 

NERDCV.)  Swathe  km) 

HRV 

xsl 

.50-.  59 

20 

.5  117* 

xs2 

.61-. 68 

20 

.5 

xs3 

.73-. 09 

20 

.5 

Pa 

.51-. 73 

10 

.5 

yPV  -  High  Resolution  Visible 

-  There  are  two  HRVs  onboard,  each  with  a  swath  of  BO  km. 
The  overlap  of  these  HRVs  results  in  this  swath  width. 


Other :  SPOT’s  HRVs  can  be  pointed  up  to  27*  off  nadir. 

This  contributes  to  the  orbit  repeat  cycle  considerably 
lower  than  any  other  commercial  sensor.  SPOT  data  can  be 
received  directly  within  2600  km  of  the  satellite.  Ground 
stations  include  Toulouse,  Kiruna  CSweden),  Prince  Albert, 
and  Gatineau.  However  only  two  of  the  four  possible  data 
streams  can  be  sent  by  the  satellite.  SPOT  can  be 
reprogrammed  within  three  to  eight  hours  for  urgent 
requests.  Follow-ons  to  SPOT  1  are  planned.  SPOT  2  and  3 
are  planned  to  be  identical  to  SPOT  1 .  For  SPOT  4  the  goal 
is  for  complete  superimposabi lity  of  panchromatic  and 
spectral  channel  data.  Also  an  ocean  surveillance  sensor  is 
under  consideration. 


Sources:  CEO;  510)  C3:  3K  11  :  151  )  C8  :  1 1 15)  C  S  ;  1H3) 


Appendix  B:  Atmospheric  Transmission 


This  appendix  contains  plots  of  transmission  verses 
wavelength  generated  by  LOUTRAN  B  For  summer  and  winter 
sub-arctic  atmospheres  CFigures  B1  through  BB) . 

The  parameters  used  to  generate  the  SUnriER  SUB-ARCTIC 
and  WINTER  SUB-ARCTIC  transmission  profiles  are  given  in 
Table  B1 . 

Major  causes  of  attenuation  for  specific  wavelengths 
are  provided  in  Table  BE. 
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Table  B1 .  LDUTRAN 

B  Parameters 

Used 

Name 

Descriotian  Summer  Value 

Winter  Value 

nODEL 

4 

5 

I  TYPE 

Slant  Path 

3 

3 

lEnSCT 

Transmittance 

0 

0 

ni 

Temp/Pressure 

0 

0 

ns 

Uiater  Vapour 

0 

0 

n3 

□zone 

0 

0 

NOPRT 

Normal 

0 

0 

TBOUND 

Temperature 

S93.0 

S53.0 

SALE 

Albedo 

0 

0 

I  HAZE 

Haze 

0 

0 

ISEASN 

Season 

0 

0 

IVULCN 

Volcanic  Ash 

0 

0 

VIS 

Visibility 

0 

0 

ICIR 

Cirrus  Clouds 

0 

0 

IVSA 

Army  VSA 

0 

0 

RAIKRT 

Rain  Rate 

0 

0 

HI 

Initial  Altitude 

0 

0 

HS 

Tangent  Altitude 

0 

0 

ANGLE 

Zenith  Angle 

0/57 

0/57 

RANGE 

Path  Length 

1000 

1000 

BETA 

Centred  Angle 

0 

0 

RO 

Earth  Radius 

0 

0 

VI 

I  nit  Wavenumber 

BOO 

BOO 

VS 

Final  Wavenumber 

S5000 

55000 

OV 

Increment 

5/10/50 

5/10/50 

I 


i 


i 


i 


i 


I 


! 
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Bandpass  C;im) 

i  0.5 

0.71-0.73 

0.76 

o.ao-o.aa 

0 . 93-0 . 95 

1.10- 1.15 

1.30- 1.50 
1.90-1.90 

2.50- a.90 
3.00-3.10 

4.30- 4.40 
5.a0-7.50 

7.50- 9.00 

9 . 10- 10.0 
>  11.00 


Attenuation  Explanation 

Molecular  Scattering 
Water 

Carbon  Dioxide 

Water 

Water 

Water,  Carbon  Dioxide 
Water,  Water  Continuum 
Water,  Water  Continuum 
Water,  Water  Continuum 
Water 

Carbon  Dioxide 
Water,  Water  Continuum 
Water,  Water  Continuum 
Ozone 

Water  Continuum 


□  0 


.  □□ 


Figure  B2 .  Sub-Arctic  Summer  Atmospheric  Transmission 

.7-1.4  Micrometers 

CLower  line  for  Zenith  angle  27  degrees) 
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Wavelength  C 


Sub-ftrctic  Summer  Atmospheric  Tr= 

CLomer  line 

For  Zenith  engie  S7  degrees) 


□□ 


Wave  I ength 


.  DO 


Figure  B5 .  Sub-Arctic  Winter  Atmospheric  Transmission 

.H-.7  Micrometers 

(Lower  line  for  Zenith  angle  27  degrees) 


1  73 
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Figure  B7 .  Sub-Arctic  Winter  Atmospheric  Transmission 

1.4-7  riicrometers 

CLouer  line  for  Zenith  angle  27  degrees) 


1  75 


Wavelength  '  <^hometers;^ 


.  QO 


Figure  BB .  Sub-Arctic  Winter  Atmospheric  Transmission 

7-12.5  flicrometers 

CLou/er  line  for  Zenith  angle  27  degrees) 


1  76 
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Appendix  C:  Visible  Sensor  Threshold  Detection  Curves 


This  appendix  contains  threshold  detection  curves  for 
the  visible  wavelength  sensors  described  in  Appendix  A. 

The  plotted  line  is  generated  using  equation  10. 
Actual  NERD  values  were  used  for  Landsat,  SPOT,  and  IRS  ' 
sensors.  For  all  other  sensors  a  nominal  NERD  of  O.S\  was 
used.  Atmospheric  transmission  values  were  taken  from  the 
curves  presented  in  Appendix  B. 

The  threshold  line  platted  represents  the  boundary 
between  detection  and  non-detection  of  the  target  by  the 
sensor.  If  the  combination  of  target  size  and  reflectance 
difference  falls  on  or  above  the  threshold  line  the  target 
will  be  detected  by  the  sensor.  If  it  Falls  below  the 
threshold  line  the  target  will  not  be  detected. 
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Figure  C2 .  SPOT  4  Target  Detection,  Winter,  Zenith  27 
CTop  to  bottom  are  Xsl,  SUIIR,  Pa  Bands) 


1  79 


Reflectance  Difference 


35,00 


Figure  C3 .  SPOT  3  Target  Detection,  Summer,  Zenith  0 
CTop  to  bottom  are  Xsl ,  XsE ,  Xs3,  Pa  BandsJ 
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ectance  Difference 


35,00 


Figure  C4 .  SPOT  4  Target  Detection,  Summer,  Zenith  27 
CTop  to  bottom  are  Xsl,  SUJIR,  Pa  Bands) 
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Figure  C5.  Landsat  Thematic  Happer,  Winter 

(Top  to  bottom  are  Bands  7,  5,  1,  3,  4) 
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Figure  CB .  Landsat  Thematic  dapper ,  Summer 

CTop  to  bottom  are  Bands  7,  5,  1,  3,  4) 
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Figure  C7 .  Landsat  riul tispectral  Scanner,  Winter 
(Top  curve  is  Band  4,  bottom  is  Band  2) 
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Ref lectance  Difference 
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Figure  C8 .  Landsat  Mul tispectral  Scanner,  Summer 
(Top  curve  is  Band  4,  bottom  is  Band  B) 
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Figure  C9 .  Landsat  Enhanced  Thematic  flapper,  Winter 

CTop  to  bottom  are  Bands  7,  5,  1,  3,  4,  Pa) 
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Figure  CIO.  Landsat  Enhanced  Thematic  flapper,  Summer 

CTop  to  bottom  are  Bands  7,  5,  1,  3,  4,  Pa) 
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Figure  Cll .  LISS-1,  Winter 
CTop  to  bottom  are  Bands  1,  B,  3,  4) 
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Figure  C12.  LIS5-1,  Summer 

CTop  to  bottom  are  Bands  1,  E,  3,  4) 
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Figure  C13.  LISS-2,  Winter 
CTop  to  bottom  are  Bands  1,  E,  3, 
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Figure  C14 .  LI55-E,  Summer 

CTop  to  bottom  are  Bands  1,  E.  3,  4) 
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Figure  CIS.  Hypothetical  10,  5,  and  1  Heter 
Ground  Resolution  Cell  Sensors 
(Top  to  bottom  are  10,  5,  and  1  meter  threshold  curves) 
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gure  CIS.  SPOT  4  HRV  0.‘i3-0.47  Sensor 
Summer  and  Winter 
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Figure  C17 ,  NDAA  K/L/tl  AVHRR  Winter 
(Top  to  bottom  are  Bands  3b,  3a,  1,  and  E) 
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Figure  CIS.  NQAA  K/L/n  AVHRR  Summer 
(Top  to  bottom  are  Bands  3b,  3a,  1,  and  2) 
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Figure  C19.  nOS-l  HESSR  Winter 
(Top  curve  is  Band  1,  bottom  is  Band  2,  3,  or  4) 
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Figure  C20,  nOS-1  MESSR  Summer 
(Top  curve  is  Band  1,  bottom  is  Band  2,  3,  or 
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Figure  C21  .  riQS-l  VTIR  Summer 
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Figure  CBS.  JERS-1  VNR  Winter 
CTop  curve  is  Band  1,  bottom  is  Band  4) 
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Figure  CS3 .  JERS-1  VNR  Summer 
CTop  curve  is  Band  1,  bottom  is  Band  3) 
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Appendix  □;  Target  Descriotion 


A  Canadian  Government  report,  Satellites  and 
Sovereignty,  published  in  August  1977,  identifies  targets 
uhich  uould  be  subject  to  surveillance  in  the  Canadian 
arctic.  Tuio  target  types  are  of  interest:  survey  parties  on 
land  or  ice;  and  shipping.  This  appendix  contains  summaries 
of  the  target  descriptions  provided  in  the  government  report 
and  representative  coating  reflectances  taken  from  the  TPRC 
Data  Series,  Volume  9,  Thermal  Radiative  Properties  — 
Coatings . 
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General : 


Survey  parties  are  made  up  of  clusters  of  vehicles  and 
shelters.  The  vehicles  can  be  expected  to  be  metallic.  The 
shelters  can  be  made  of  a  variety  of  materials.  The  number 
of  vehicles  and  shelters  uill  vary. 

Spatial  Characteristics: 

Spatial  characteristics  uill  also  vary.  For  analysis  a 
representative  party  consisting  of  10  vehicles  and  Five 
shelters  uiill  be  used.  Vehicles  are  assumed  to  be  pairs  of 
10,  20,  30,  HO,  and  50  square  meters  plan  size  when  vieued 
from  nadir.  Shelters  for  analysis  are  tents  10  square 
meters  plan  size  when  vieued  from  nadir.  When  in  camp  all 
vehicles  and  shelters  are  assumed  to  be  clustered  uithin  a 
500  square  meter  area.  When  deployed  vehicles  are  assumed 
to  be  in  column  and  spaced  by  100  meters. 

These  choices  are  arbitrary  but  uithin  the  guidelines 
established  in  the  government  report  previously  cited. 

Spectral  Characteristics : 

Vehicles  and  shelters  are  assumed  to  be  painted  US  Army 
olive  drab  in  sufficient  quantity  to  provide  an  optically 
thick  coating.  The  reflectance  of  this  paint,  as  a  Function 
of  uavelength,  is  plotted  in  Figure  D1 . 

It  is  assumed  that  1/4  of  vehicle  plan  size  is  E*K 
uarmer  than  the  background  uhile  vehicles  are  operating. 
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Figure  D1 .  Reflectance  of  US  Army  Olive  Drab 
vs  Wavelength  CTPRC:5B9) 


Vehicles  are  assumed  to  be  in  constant  operation.  Shelters 
are  assumed  to  be  at  ambient  temperature. 

Temporal  Characteristics: 

Base  camps  are  assumed  to  be  stationary  but  can  be 
moved  monthly. 

Information  Requirements : 

Images  are  required  daily  and  must  be  available  within 
IE  hours  of  collection. 
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Targets  of  interest  are  single  ships  entering  and 


travelling  in  Canadian  arctic  uiaters. 

Spatial  Characteristics: 

Ships  may  be  10  to  200  meters  in  length.  For  this 
analysis  ships  mill  be  a  DDK  class  destroyer  uhich  is  115 
meters  long  and  15  meters  wide  C37;73). 

Spectral  Characteristics : 

Ships  will  be  made  of  steel  and  painted  wi  '■  aluminum 


0  2  ■«»  6  8  10  12  14 

wavelength  cmlcrometers^ 


Figure  DE .  Reflectance  of  Aluminum  Oxide  on  Steel 
vs  Wavelength  C70:735) 


[ 

I 

i 

i 


oxide  contact  coating.  The  reflectance  of  this  paint,  as  a 
function  of  wavelength,  is  plotted  in  Figure  DE . 

One  tenth  of  the  plan  size  of  the  ship  is  assumed  to  be 
B*K  warmer  than  ambient  temperature. 

Temporal  Characteristics: 

The  ship  is  assumed  to  travel  at  SO  knots  in  open 
water ,  and  at  3  knots  through  ice . 

Information  Requirements : 

Images  are  required  every  IE  hours  and  must  be 
available  within  IS  hours  of  collection. 
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lendix  E:  Threshold  Longitudes  for  Landsat  and  SPOT 


This  appendix  contains  threshold  longitude  curves  for 
overlap  of  the  possible  target  area  by  Landsat  and  SPOT 
sensor  siuaths.  If  the  sensor  passes  uiithin  the  threshold 
longitude,  the  sensor’s  siuath  and  the  possible  target  area 
uill  overlap. 

To  determine  the  local  crossing  angle,  the  follouing 
equations  are  used; 

cosC30-q)  =  casCLong)sin( i ) 

and 

sinCLong)  =•  tanCLat)/tanC  i  ) 

Where 

q  =  local  crossing  angle 

Long  =  longitude 

i  =  inclination 

Lat  -  latitude 
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igure  El.  Threshaid  Longitude  For  Overlap 
Landsat  5,  60’N  Latitude 
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Figure  EE.  Threshold  Longitude  for  Overlap 
Landsat  5,  70*N  Latitude 
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Figure  E3 .  Threshold  Longitude  For  Overlap 
Landsat  5,  0O*N  Latitude 
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igure  E4 .  Threshold  Longitude  for  Overlap 
SPOT  1,  BO*N  Latitude 
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This  appendix  contains  plots  of  conditional  tracking 
probabilities  For  Landsat  5  and  SPOT  sensors. 

One  anomaly  of  the  spreadsheet  used  to  generate  these 
curves  is  that  achieving  a  conditional  tracking  probability 
of  1.0  Forces  the  softujare  to  expand  the  y  axis  to  a  maximum 
score  of  l.S.  Since  a  probability  in  excess  of  1.0  is 
meaningless,  conditional  tracking  probabilities  exactly 
equal  to  1.0  are  restricted  to  a  value  of  0.93  in  this 
appendix  but  should  be  taken  to  be  1.0. 
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Figure  FI.  Landsat  5  Conditional  Tracking  Probability 
60*N  Latitude,  5  Knots 

CTop  to  bottom  are  0 , 1 , 2 , 3 , 4 , 5*  Longitude  Displacement) 
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Figure  F2.  Landsat  5  Conditional  Tracking  Probability 
B0*N  Latitude,  10  Knots 

(Top  to  bottom  are  0,1,2,3,H,5*  Longitude  Displacement) 
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Figure  F3 .  Landsat  5  Conditional  Tracking 
BO'N  Latitude,  20  Knots 
CTop  to  bottom  are  0,1, 2, 3, 4, 5*  Longitude  D 
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Figure  F4 .  Landsat  5  Conditional  Tracking  Probability 
70*N  Latitude,  5  Knots 

CTop  to  bottom  are  0,1,3,4,5,10*  Longitude  Displacement) 
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Figure  F5 .  Landsat  5  Conditional  Tracking  Probability 
70*N  Latitude,  10  Knots 

CTop  to  bottom  are  0,1,3,4,5,10*  Longitude  Displacement; 
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Figure  F6 .  Landsat  5  Conditional  Tracking  Probability 
70*N  Latitude,  20  Knots 

CTop  to  bottom  are  0,1,3,4,5,10'  Longitude  Displacement) 
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Figure  F7 .  SPOT  1  Conditional  Tracking  Probability 
60'N  Latitude,  5  Knots 

CTop  to  bottom  are  0,1, 2, 3, 4, 5*  Longitude  Displacement) 
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Figure  F0 .  SPOT  1  Conditional  Tracking  Probability 
60*  N  Latitude,  10  Knots 

Top  to  bottom  are  0,1, 2, 3, 4, 5*  Longitude  Displacement) 


Figure  F3 .  SPOT  1  Conditionai  Tracking  Probability 
BO*  N  Latitude,  80  Knots 

CTop  to  bottom  are  0,1,E,3,H,5*  Longitude  Displacement j 
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Figure  FIO.  SPOT  1  Conditional  Tracking  Probability 
70*  N  Latitude,  5  Knots 

CTop  to  bottom  are  0,1, 2, 3, 4, 5*  Longitude  Displacement) 
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Figure  Fll.  SPOT  1  Conditional  Tracking  Probabiiity 
70*  N  Latitude,  10  Knots 

CTop  to  bottom  are  0,1,S,3,4,5‘  Longitude  Displacement) 
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Figure  F12.  SPOT  1  Conditional  Tracking  Probability 
70*  N  Latitude,  20  Knots 

(Top  to  bottom  are  0,1, 2, 3, 4, 5*  Longitude  Dispiacement ) 
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This  appendix  contains  plots  of  threshold  target  sizes 
for  detection  as  a  function  of  target  exitance  for  specific 
combinations  of  targets  and  background.  Combinations  of 
target  exitance  and  area  which  fall  below  the  threshold 
curve  will  not  be  detected  by  the  sensor.  Those  that  fall 
on  or  above  the  line  will  be  detected  by  the  sensor. 

The  conditions  specified  below  each  graph  are  important 
since  the  graphs  only  apply  for  the  indicated  target  and 
background  conditions. 
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Figure  G1 .  Threshold  Target  Area  vs  Target  Exitance 
Target  Emissivity  .92,  Background  Emissivity  .94 
Landsat  5  Band  6  9ensor ,  Background  Temperature  253*K 
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Target  Exitance  11^23 


Figure  GE .  Threshold  Target  Area  vs  Target  Exitance 
Target  Emissivity  .3S,  Background  Emissivity  .34 
Landsat  5  Band  B  Sensor,  Background  Temperature  E75*K 
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Target  Exitance 


Figure  G3 .  Threshold  Target  Area  vs  Target  Exitance 
Target  Emissivity  .9E,  Background  Emissivity  .94 
Landsat  5  Band  6  Sensor,  Background  Temperature  E93*K 
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Figure  GH .  Threshold  Target  Area  vs  Target  Exitance 
Target  Emissivity  .92,  Background  Emissivity  .995 
Landsat  S  Band  6  Sensor,  Background  Temperature  253*K 
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Figure  G5 .  Threshold  Target  Area  vs  Target  Exitance 
Target  Emissivity  .92,  Background  Emissivity  .905 
Landsat  5  Band  B  Sensor,  Background  Temperature  275*K 
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Figure  GB .  Threshold  Target  Area  vs  Target  Exitance 
Target  Emissivity  .74,  Background  Emissivity  .99 
Landsat  5  Band  6  Sensor,  Background  Temperature  975*K 


232 


Figure  G7 .  Threshold  Target  Area  vs  Target  Exitance 
Target  Emissivity  .74,  Background  Emissivity  .99 
Landsat  5  Band  B  Sensor,  Background  Temperature  293*K 
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T arget  Exitance 


Appendix  Hr  Analysis  Checklists 


This  appendix  contains  the  checklists  used  to  calculate 
individual  remote  sensing  passive  sensors’  probabilities  of 
detection  and  tracking,  and  to  optimize  system  use  of  these 
sensors . 
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Checklist  1  —  Start 

1.  Specify  mission,  targetCs)  of  interest,  and  timeframe. 

2.  Define  area  of  interest. 

3.  Identify  backgrounds  possible  in  area  of  interest  during 

timeframe . 

4.  Describe  background’s  spatial,  spectral,  and  temporal 

characteristics . 

5.  Describe  target’s  spatial,  spectral,  and  temporal 

characteristics . 

6.  Identify  sensors  available  during  timeframe. 

7.  Calculate  Par  each  combination  of  sensor,  target, 

and  background  using  Equations  10  or  EH., 

CSee  Appendices  C  and  G) 

0.  Discard  all  combinations  with  P^g^  =  0. 

9.  Gn  to  mission  checklist. 
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Checklist  S  —  Non-Cued  Detection  oF  Stationaru  Targets 


1 . 


E. 


3. 

4. 

5. 
5. 

7  . 


8. 


8. 


Calculate  considering  array  geometry  using 

Equations  SB,  E7,  SB  and  S9. 

Partition  area  of  interest  into  areas  of  common 
obscuration . 

For  each  zone,  target,  background,  and  sensor  calculate 
the  Pjjgt*  using  Equations  35  or  3B . 

Discard  all  combinations  with  P^jg^*  “  0. 

Calculate  each  combination’s  P(jQtl  using  Equation  41. 
Specify  if  the  probability  of  detection  is  to  be 
maximized  or  if  the  cost  is  to  be  minimized. 

IF  goal  is  to  maximize  P^^gyg: 

A.  Equation  4S  used  as  objective  function,-  and 

B.  Equations  43,  44,  45  and  specification  that  each 
area  must  return  similar  Pjjgyg  make  up  the 
constraint  set. 

ELSE  goal  is  to  minimize  cost: 

A.  Equation  43  is  che  objective  function;  and 

B.  Equations  4S ,  44,  45  and  specification  that  each 
area  must  return  similar  P^jgyg  make  up  the 
constraint  set. 

Solve  the  linear  program. 


i 
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1.  Calculate  Pjjg^  considering  array  geometry  using 
Equations  SB,  57,  58  and  53. 

5.  Partition  area  of  interest  into  areas  of  common 
obscuration . 

3.  For  each  zone,  target,  background,  and  sensor  calculate 

the  '-‘Sing  Equations  35  or  36. 

4.  Discard  all  combinations  with  P^jg^*  *  0* 

5.  Calculate  each  combination’s  P^g^l  Equation  41. 

B.  Specify  if  the  probability  of  detection  is  to  be 

maximized  or  if  the  cost  is  to  be  minimized. 

7.  IF  goal  is  to  maximize 

A.  Equation  47  used  as  objective  function;  and 

B.  Equations  43,  45  and  specification  that  each 
area  must  return  similar  P^>i  make  up  the 
constraint  set. 

C.  Solve  using  non-linear  programming, 
a.  ELSE  goal  is  to  minimize  cost; 

A.  Equation  43  is  the  objective  function;  and 

B.  Equations  45,  47  and  specification  that  each 

area  must  return  similar  make  up  the 

constraint  set. 

C.  Solve  using  non-linear  programming. 
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Checklist  4  —  Cued  Detection  of  Stationary  Targets 

1.  Identify  sensors  which  will  have  imaging  opportunity  in 
timeframe . 

S.  Calculate  considering  array  geometry  using 

Equations  B6,  S7,  BS  and  B9. 

3.  Calculate  P(jgt»  using  Equations  35  or  36 . 

4.  Specify  if  the  probability  of  detection  is  to  be 

maximized  or  if  the  cost  is  to  be  minimized. 

5.  IF  goal  is  to  maximize  Pjj[>i: 

A.  Equation  48  used  as  objective  function;  and 

B.  Equations  43,  and  45  make  up  the  constraint  set, 

C.  Solve  using  non-linear  programming. 

B.  ELSE  goal  is  to  minimize  cost: 

A.  Equation  43  is  the  abjective  function;  and 

B.  Equations  45,  and  48  make  up  the  constraint  set. 

C.  Solve  using  non-linear  programming. 
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Checklist  5  —  Second  Contact 


1 .  Determine  longitude  and  time  of  crossing  of  last  known 
target  latitude  for  all  sensors  throughout  timeframe. 

S.  Determine  if  overlap  of  swath  and  possible  target  area 
overlap  using  Equation  53.  Cor  see  Appendix  E) 

3.  Calculate  Pjjg^  considering  array  geometry  using 

Equations  26,  27,  28  and  29. 

4.  Calculate  Pjjgt*  using  Equations  35  and  36. 

5.  Discard  options  with 

6.  IF  sensor  is  Landsat  5  or  SPOT  1; 

A.  See  Appendix  F  to  determine  Ptrkcon ' 

ELSE 

A.  Determine  value  of  Y  using  Equations  54  and  55. 

B.  IF  Y  <  half  the  swath  width  of  the  sensor: 

1.  Determine  using  Equations  57,  50,  59,  BD, 

61,  62,  63  and  64. 

2.  Determine  using  Equation  73. 

ELSE  Y  >  half  the  sensor’s  swath  width: 

1.  Determine  using  Equations  65,  66,  67,  BB , 

69,  70,  71  and  72. 

2.  Determine  Ap^^^  using  Equation  73. 

ALWAYS 

3.  Determine  Ptrkcon  using  Equation  74. 

7.  Determine  using  Equation  75. 

0.  Determine  P^rackl  using  Equation  76. 

9.  IF  available  swaths  do  not  cfverlap: 

A.  Specify  if  Ptrack  maximized  or  cost 

minimized . 

B.  IF  Ptrack  maximized: 

1.  Equation  77  is  the  objective  function. 

2.  Equations  70  and  79  are  the  constraint  set. 

C.  IF  cost  to  be  minimized; 

1.  Equation  70  is  the  objective  Function. 
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S.  Equations  77  and  73  are  the  constraint  set 
0.  Solve  using  linear  programming. 

10.  IF  available  suaths  do  overlap: 

A.  Specify  if  ^track  maximized  or  cost 

minimized . 

B.  IF  P^rack  maximized: 

1.  Equation  81  is  the  objective  function. 

S.  Equations  BE  and  83  are  the  constraint  set 
3.  Solve  using  non-linear  programming. 

C.  IF  cost  to  be  minimized: 

1.  Equation  88  is  the  abjective  function. 

S.  Equations  81  and  83  are  the  constraint  set 
3.  Solve  using  non-linear  programming. 


I 


! 
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1 .  Determine  longitude  and  time  of  crossing  of  last  knoun 
target  latitude  for  all  sensors  throughout  timeframe. 
B.  Determine  earliest  target  penetration  of  swath  and  suiath 
transit  time  using  Equations  S4  and  B5. 

3.  Identify  all  imaging  opportunities  in  timeframe. 

H.  Calculate  '-•sing  Equations  35  and  36. 

5.  Specify  i  f  Pd>  ^  to  be  maximized  or  cost  minimized. 

A.  IF  Pjj>2  to  be  maximized: 

1.  Equation  06  is  the  objective  function. 

2.  Equations  07  and  06  are  the  constraint  set. 

3.  Solve  using  non-linear  programming. 

B.  IF  cost  to  be  minimized; 

1.  Equation  07  is  the  objective  function. 

B.  Equations  06  and  SB  are  the  constraint  set. 

3.  Solve  using  non-linear  programming. 
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